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The  purpose  of  this  report  is  to  provide  a basis  for  assessing  the  merits  of  a proposal  by 
General  Electric  Corp.  to  detect  acoustic  signals  in  the  water  by  a laser  interferometer  in  which 
the  laser  beam  itself  is  totally  submerged.  The  basis  we  will  use  for  making  this  assessment  will 
be  to  compare  the  proposed  device  with  state-of-the-art  Navy  hydrophone  capability.  For  ready 
reference  we  take  our  model  of  (one  of)  the  best  Navy  hydrophones  available  to  daw  to  be  the 
USRD  H56,  whose  relevant  performance  we  summarize  as  follows: 

Henriquez  (J.  Acous.  So.  Am.  52,  p 1450  (1972))  discusses  the  design  of  the  H56  hydro- 
phone, which  he  constructed  to  measure  ambient  noise  in  10  Hz  to  60  kHz  band.  The  sensitive 
element  is  a tangentially  poled  PZT-5  cylinder,  of  which  two  are  incorporated  in  a single  unit. 

The  equivalent  noise  pressure  level  of  the  H56  is  38  dB  re  i |ifa.in  a 1 Hz  band  at  100  Hz.  The 

free-field  voltage  sensitivity  is  -185  dB  re  . Thus  the  internal  noise  of  the  H66  (in  a 1 Hz- 

-l.'SS  ' _ 

band)  appears  as  10  volt  at  the  terminals.  Knudsen  seastate  zero  is  62  dB  re  i Ufa.  in  a 1- 

IS  1 

Hz  band  at  100  Hz.  This  noise  in  the  water  appears  as  10  volt  at  the  terminals.  Thus  sea- 
state  zero  is  greater  than  the  internal  noise  of  the  H66  by  a factor  of  (approximately)  10  in  the 
frequency  region  less  than  1000  Hz,  making  the  H56  ideally  suited  to  the  study  of  ambient  sea 
noise. 


The  H56  active  PZT  element  lias  a capture  area  (i.e.,  the  area  which  “sees”  the  wave  front) 
of  approximately  200  It  mm.  ln  a typical  single  beam  gas  laser  experiment  of  the  doppler- 
velocimeter  type  the  capture  area  is  about  iSlt  MO  . The  ratio  of  capture  areas  is  thus 

about  13,000  in  favor  of  the  H56.  The  H56  is  a pressure-sensitive  device.  The  proposed  virtual 
acoustic  sensor  is  a particle  velocity-sensitive  device.  Thus  the  H56  measures  the  scalar  aspect  of 
the  acoustic  field  while  the  proposed  device  measures  (one  component  of)  the  vector  aspect  of. 
the  acoustic  field. 


The  chief  problem  of  the  virtual  acoustic  sensor  is  minuteness  of  the  threshold  signal  in 
comparison  to  the  dominating  noise  of  the  medium  and  of  the  detection  circuit  itself.  The 
following  report  is  devoted  to  an  analysis  of  this  signal-to-noise  ratio  problem. 
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Summary  and  Conclusions 

We  list  below  the  chief  conclusions  of  this  investigation. 

I.  The  method  of  detection  of  an  acoustic  signal  reviewed  in  this  report  is  that  of 
a.  laser  heterodyne  device  designed  to  measure  either  of  two  different  quantities  (1)  a sinus- 
oidal displacement  of  particles  in  a fluid  (2)  a fluid  velocity  (turbulent  or  laminar).  In  the 
first  case  the  success  in  detection  rests  on  the  capability  of  measuring  magnitude  of  the 
power  spectrum  of  a photodetector  current  in  the  presence  of  noise,  rather  than  in  the 
capability  of  measuring  a frequency  shift.  The  important  physical  quantity  in  the  signal 
processing  is  the  modulation  index,*  not  a Doppler  shift.  Hence  the  acoustic  sensor  in 
question  is  a true  displacement  device,  rattier  than  a velocity  device.  I.i  the  second  case 
the  acoustic  sensor  is  designed  to  measure  fluid  velocity.  Success  in  detev/don  rests  on  the 
capability  of  measuring  Doppler  shift  in  the  presence  of  noise.  It  is  a true  Laser  Doppler 
Velocimeter. 

Conclusions  on  the  two  applications  are  presented  below. 

II.  As  an  acoustic  sensor  of  displacement  we  estimate  the  capabilities  of  the  laser 
heterodyne  as  follows: 

a.  The  theoretical  magnitude  of  modulation  index  which  meets  the  Navy  threshold 
requirements  for  detection  of  a submarine  in  sea  state  1 at  100  Hz  is  of  the  order  of  8 X 10"® 
radian.  The  possibility  of  detecting  an  index  of  this  small  magnitude  is  the  core  of  the 
feasibility  study  in  the  accompanying  report. 

b.  Under  the  assumption  that  shot  noise  is  the  only  noise  in  the  circuitry  of  an 

acoustic  displacement  sensor  Massey  (1968,  Proc.  IEEE  56,  2157)  calculated  that  displace- 
-n  t 

ments  of  the  order  of  ICT  M could  be  measured  in  the  laboratory  for  a laser  wavelength 

-i° 

of  (o  *>"30»  10  fie  equivalent  to  a modulation  index  of  approximately  2 X 10"5  radian. 
Thus  the  Navy  threshold  to  be  achieved  is  (somewhat  smaller  than)  an  order  of  magnitude  less 
than  the  displacement  laser  heterodyne  capability  calculated  by  Massey  to  be  available  in  the 
laboratory  in  the  presence  of  shot  noise. 


•For  definition  tee  page  71. 
fM  ■ meter. 


c.  Power  spectral  broadening  can  be  due  to  causes  other  than  shot  noise.  When 

non-shot  noise  is  the  dominant  feature,  the  detectable  modulation  index  is  shown  by 
mathematical  modeling  to  be  a function  of  the  ratio  of  the  acoustic  frequency  to  the 
total  noise  bandwidth  4*  (in  units  of  radians/sec).  Orders  of  magnitude  calculations  show 
the  following:  when  the  bandwidth  of  noise  in  a power  spectrum  is  larger  than  the 
acoustic  frequency11,  the  modulation  index  is  of  the  order  of  unity.  The  smallest  displace- 
ment that  can  be  detected  is  then  about  4 X 10~8  meter.  This  is  five  orders  of  magnitude 
greater  than  the  Navy  threshold.  In  contrast  when  the  noise  bandwidth  in  a power 
spectrum  is  a fraction  of  the  acoustic  frequency  the  detection  capability  is  much  greater. 
Calculation  shows  that  a ratio  ***>  /ip  **  7 corresponds  to  a capability  of  detecting  a 
modulation  index  of  10*5,  which  is  neany  the  Navy  threshold.  This  means  that  at  an  acoustic 
frequency  of  100  Hz  the  noise  bandwidth  muft  not  exceed  14  Hz  in  order  to  measure 
acoustic  displacements  of  the  order  of  10"^  to  lO"^  meter.  The  possibilty  of  reducing 
noise  banc'width  to  the  limit  is  open  to  investigation. 

d.  In  all  calculations  of  c.  we  have  assumed  that  the  received  signal  from  the 
scattering  volume  is  large  enough  to  overcome  the  inherent  noise  in  the  photodetector 
circuit.  Under  certain  simplifying  assumptions  this  noise  is 


(e  ■ electron  charge,  “ laser  wavelength,  zj  “ detector  quantum  efficiency,*^  is  the 
received  power,  5 • detector  circuit  bandwidth  ^ h-  Planck’s  constant,  ■ speed  of 
light). 

The  received  power  collected  over  an  areata  in  direction  9 is  given  by 


(-a(V'Ki))  (*a)^ 


<?(*)  * scattering  function,"^  * laser  power,  0 i m attenuation  coefficient, 7?, distance 
to  / from  the  scattering  volume),  lm  depth  of  scattering  volume). 

Combining  the  two  equations  shows  that  the  minimum  detectable  displacement  (to  an 
order  of  magnitude)  is 


•This  i*  approximately  the  caae  of  a angle  laser  beam  disturbed  by  noisefat  100  Hz)due  to  Brownian  motion 
(aee  Part  IV  ).  Multiple  beams  (or  “diversity")  may  remove  this  limit  (see  Part  VI.). 


2 


*sniN 


£ 

4* 


( ) 


♦ 


At  a distance  of  30  meters  in  sea-water  the  backscattered  (i.e.  6 * ISO0)  laser  power  into 
a lens  of  diameter  0.1  meter  is  of  the  order  of  10*8  P , where  "T  is  the  laser  output  power. 
Hence  the  minimum  detectable  signal  (order  of  magnitude)  is 


Hence  we  conclude  that  to  achieve  the  Navy  threshold  at  a distance  of  30  meters  one 
requires  a laser  power  whose  magnitude  in  watts  is  about  equal  to  the  detector  circuit 
bandwidth  magnitude  in  herz,  provided  the  receiving  aperture  has  a diameter  of  0.1  meter. 
This  required  laser  power  can  be  reduced  by  increasing  the  receiving  aperture  and/or  the 
depth  of  the  scattering  volume.  The  possibility  of  increasing  receiving  aperture  is  very 
important  to  the  success  of  the  laser  heteodyne  detector,  and  is  open  to  investigation. 

e.  The  effect  of  Brownian  motion  is  serious  in  that  it  sets  an  irreducible  magnitude 
to  the  noise  bandwidth  A B noted  in  c.  above.  However  this  is  true  only  if  the  capture 
area  (i.e.  the  scattering  volume)  is  very  small.  It  is  well-known  that  increasing  capture  area 
reduces  the  effects  of  random  inputs  into  a detection  system  by  affording  an  opportunity 
for  increasing  the  ratio  (see  “diversity,”  Part  VI.). 

g.  The  effect  of  “platform  motion”  is  serious  only  if  this  motion  is  correlated  to,  or 
has  a Fourier  component  in  the  same  pass  band,  as  the  acoustic  sigrsl  to  be  processed. 

h.  The  effect  of  medium  inhomogeneities  including  gas  bubbles  is  serious  in  that  it 
degrades  the  laser  beam  coherence  over  long  path  length.  However  such  degradation  can 

be  overcome  by  increasing  the  receiver  capture  area.(sef;  Hodara  (1966)  Proc.  IEEE  54,  368) 

i.  The  acoustic  particle  displacement  (hence  pa.ticle  velocity )measured,  is  one  com- 
ponent of  a three-component  vector,  namely  the  component  in-line  between  receiver  and 
scattering  volume.  Two  component  LDV’  « have  fcwen  constructed  (Created  1971,  J.  Phys. 

E.  4,  585;  Blake  1972,  -J.  Phys.  E.  5,  623;  Grant  Mid  Orioff  (1973),  Appl.  Optics  12,  2913). 
The  vector  nature  of  particle  displacement  does  not  appear  to  be  a ragjor  problem. 

j.  A space  array  of  laser  beams  (•  virtual  Array)  is  essentially  an  assembly  of  electro- 
optic hydrophones  of  dipole  type.  The  signal  processing  of  the  returns  from  such  an  array 


*For  derivation  aee  page  69. 


is  conventional.  The  virtual  array  has  the  same  size  as  a real  array  at  the  same  frequency. 

k.  Since  a virtual  array  of  laser  beams  measures  the  local  particle  displacement  the 
effects  of  reflecting,  diffracting  or  scattering  bodies  must  be  taken  into  account.  1116 
virtual  array  must  be  located  far  enough  away  from  such  bodies  in  order  to  give  a true 
statement  of  the  incoming  acoustic  particle  motion. 

III.  As  a sensor  of  fluid  velocity  we  estimate  the  capabilities  of  the  laser  heterodyne  as 
follows: 

a.  Yeh  and  Cummins  (1964  Appl.  Phys.  Letters  4,  176)  concluded  that  they  could 
detect  (in  the  laboratory)  constant  (i.e.  laminar)  velocities  as  low  as  4 X 10"®  meter/sec 
at  a scattering  angle  of  30°. 

b.  Edwards  et.  al  (1971,  J.  Appl.  Phys.  4^  837)  estimated  that  under  conditions  (in 
the  laboratory)  where  thermodynamic  diffusion  of  molecules  was  the  limiting  factor  they 
could  detect  constant  velocities  as  low  as  10*^  meter/sec. 

Both  of  these  values  are  five  orders  of  magnitude  greater  than  the  Navy  threshold. 

c.  The  basic  limit  in  the  use  of  the  laser  Doppler  velocimeter  for  measurement  of 
turbulent  velocity  fluctuations  is  the  Doppler  ambiguity  (or  uncertainty  in  measuring  a 
frequency  shift)  due  to  extraneous  time-varying  modulation  of  the  laser  beam.  These 
modulations  are  introduced  by  finite  transit  time  of  particles  through  the  scattering  volume, 
turbulent  fluctuations  across  the  scattering  volume,  mean  velocity  gradients,  and  circuit 
noise.  Doppler  ambiguity  limits  spatial  and  temporal  resolution. 

d.  The  only  measurable  velocity  inturbulent  flow  is  the  Eulerian  random  velocity 
1^(4)  averaged  over  the  scattering  volume.  This  is  the  sum  o i a mean  velocity  t to  (*) 
and  a fluctuating  velocity  U«(0*  Th*  power  spectrum  of  turbulence  consists  of  a mean 
(Doppler)  frequency  shift  broadened  by  the  spectrum  of  the  fluctuating  components.  As 
noted  in  c.  the  resolution  of  the  power  spectrum  of  turbulence  (that  is,  its  separation 
out  of  the  noise)  is  limited  by  the  Doppler  ambiguity  (DA).  If  the  frequency  broadening 
of  the  turbulent  velocity  fluctuations  (namely  the  quantity  we  wish  to  measure)  is  of  the 
same  order  as  the  broadening  due  to  Doppler  ambiguity  (which  is  the  noise  we  wish  to 
avoid)  then  there  is  no  way  of  telling  them  apart  If  *)«  is  the  mean  Doppler  shift  due 
to  ifftft)  then  the  condition  of  resolution  is 
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(This  Doppler  ambiguity  poses  a severe  limit  in  the  determination  of  the  turbulence  spectrum.) 
It  is  fundamental  to  recognize  that  the  measurement  of  laminar  flow  (Z  which  is  non-random 
differs  from  the  measurement  of  t’irbulent  flow  which  is  random.  In  the  latter  case 
there  is  a largest  wavenumber  (or  highest  cut-off  frequency)  that  is  measurable  for  a fixed 
Reynolds  number  and  fixed  scattering  angle.  Thus  the  entire  power  spectrum  of  velocity 
turbulence  is  unattainable.  A simple  estimate  of  the  largest  measurable  wavenumber  is 
in  which  L is  the  largest  dimension  of  the  scattering  volume.  Thus  if  L is  a number 
fixed  by  the  LDV  the  largest  turbulence  wavenumber  measurable  is  L , and  the  rest  of  the 
spectrum  is  unresolvable.  Hence  if  the  presence  of  submarine  turbulence  is  to  be  determined 
by  examining  turbulence  scale  sizes  less  than  L meters,  the  LDV  method  fails.  It  can  be 
revived  by  reducing  the  scattering  volume.  However,  such  reduction  is  accompanied  by 
increase  in  Doppler  ambiguity  since  space  is  sampled  over  a shorter  time  interval.  Mathemat- 
ical modeling  shows  that  there  is  an  optimum  size  of  scattering  volume,  , given  by 


( € •*  rate  of  dissipation  of  turbulent  energy  per  unit  man,  V - kinematic  viscosity, 
Reynolds  number  based  or.  the  mean  velocity  * ),  • * angle  of  scattering). 

Wavenumbers  greater  than  «*  not  resolvable  because  of  Doppler  ambiguity.  The 

symbol  is  the  inner  scale  (meters)  of  turbulence.  When  6 ■ 180°, 


It* 

so  that 


Vo  (S>/q7 


0 Z5tl. 


L 


•pr 


L (o*s7  )(o-s)  Uo4  » 


If  the  turbulent  velocity  is  1 meter/sec,  the  optimum  scattering  volume  is  1/4  of  the  inner  scale  of 
turbulence.  This  is  a very  severe  restriction.  Any  attempt  to  decrease  the  scattering  volume  only 
increases  the  Doppler  ambiguity. 
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e.  George  and  Lumiey  (1973,  J.  Fluid  Mech.  60,  p.  321)  state  that  “ . . . estimates 
show  that  the  possibility  of  measuring  dissipation  spectra  in  high-speed  or  in  geophysical 
flows  using  Doppler  velocimeters  is  quite  remote.” 

Of  course  the  detection  of  submarine-induced  velocity  turbulence  may  not  depend  on 
determination  of  dissipation  spectra.  It  will  however  depend  on  determining  some  portion 
of  the  velocity  spectrum  (versus  wavenumber).  Hie  portion  that  can  possibly  be  used  is 
a subject  of  investigation. 

f.  The  basic  analysis  (presented  in  Part  II  of  this  report)  emphasizes  temporal  corre- 
lation of  the  velocity  spectrum  at  a single  point  in  space.  By  use  of  two  velocimeters  to 
sample  different  scattering  volumes  one  may  make  two-point  velocity  correlations.  In  this 
way  additional  statistical  moments  of  the  velocity  turbulence  can  be  obtained.  These  may 
be  essential  in  detecting  the  presence  of  a wake  of  a submarine. 

Two-point  velocity  correlations  in  a turbulent  fluid  have  been  successfully  measured 
by  CIark(1970  Ph.  D.  Thesis,  U.  of  Virginia,  Charlottesville). 

It  is  to  be  noted  that  single-point  statistics  are  insufficient  to  characterize  turbulent 

flow. 

g.  C.  J.  Bates  (July  1974,  DI8A  INFO.  No.  16,  p.  6-10,  779  Susquehanna  Ave., 

Franklin  Lakes,  New  Jersey  07417)  has  studied  Doppler  ambiguity  bandwidths  in  (laboratory) 
pipe  flow  of  water  using  an  LDV.  In  a 10  inch  pipe  at  flow  rate  of  1.4  ft/sec  (Reynolds 
number  2.066  x 10®)  he  found  the  following  ratios  of  spectral  broadening: 


Core 

Wall 

(1) 

transit  time  spectral  broadening 
mean  Doppler  frequency  shift 

0.636 

0.636 

(2) 

turbulence  spectral  bioadening 
transit  time  spectral  broadening 

1.074 

38 

(3) 

mean  velocity  gradient  spectral  broadening 
turbulence  spectral  broadening 

0.000719 

0.76 

Thus,  in  accordance  with  lie.  above  the  smallest  unambiguous  displacement  that  can 
be  measured  in  Bates'  experiment  is  estimated  to  be  10*®  meter. 
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IV.  The  key  problem  of  the  laser  heterodyne  detector  of  low  level  acoustic  signals  is  its  very 
small  S/N  ratio.  To  improve  this  ratio  General  Electric  Co.  proposes  to  use  the  diversity  tech- 
nique. Diversity  is  discussed  in  Part  VI  of  this  report.  Successful  application  of  diversity  is 
based  on  these  assumptions:  (1)  a multiplicity  of  “channels”  (or  copies)  of  the  transmitted 
signal,  (2)  Rayleigh  fading  in  each  channel,  (3)  the  received  signal  plus  noise  are  statistically 
independent. 

In  the  laser  heterodyne  detector  the  incoming  signal  plus  noise  is  generally  unknown,  al- 
though pattern  recognition  leads  to  positive  identification.  The  noise  of  the  system  is  primarily 
due  to  the  medium  itself,  and  is  correlated  over  significant  distances  in  the  medium  for  many 
frequencies  of  interest.  There  is  also  some  question  as  to  whether  the  amplitudes  of  the  received 
signals  are  Rayleigh  distributed. 

The  application  of  diversity  thus  will  require  more  insight  into  the  fading  properties  of  the 
transmission  channel  from  laser  to  field  points  and  back.  If  the  spatial  sampling  is  at  least  a 
half-wavelength  of  the  acoustic  frequency  to  be  detected,  and  if  the  fading  is  statistically  inde- 
pendent between  channels  it  is  anticipated  that  diversity  transmission  will  increase  the  S/N  ratio 
significartly.  The  nominal  increase  is  10  logioN,  where  N is  the  number  of  channels. 

V.  Rough  estimates  of  the  detection  limits  of  laser  heterodyne  systems  may  be  made  with  knowl- 
edge of  times,  distances,  absorptions,  etc.,  which  are  of  significance  to  the  generation  and  trans- 
mission of  laser  light  in  water.  These  estimates  appear  below. 


Significant  Distances  and  Times 


RMS  acoustic  particle  displacement,  planewave  equivalent  to 


seastate  1 at  100  Hz 

3.4  x 10"13M 

Radius  of  inner  most  electron  orbit  of  the  hydrogen  atom 

0,5  x 10-10M 

Laser  frequency 

5.8  x 1014Hz 

Laser  wavelength 

5.145  x 10-?M 

Sound  Speed  in  water 

1.5  x 103M 

One  period  of  acoustic  wave  at  100  Hz 

1 x 10"2s 

Speed  of  light  in  water 

2.26  x 103Ms-1 

One  period  of  lasr  r light 

1.72  x 10"15s 

2-way  travel  time  over  30  M range  in  water 

2.65  x 10~7s 

Laser  wavelength  in  water 

3.86  x 10"7M 

Time  to  traverse  one  quarter  laser  wavelength  in  water 
Depth  of  Volume  Interrogation  *hJf  laser  wavelength)  in 

0.426  x 10"15s 

water(this  depth  can  be  made  longer). 

1.93  x 10-7M 

Ratio  of  laser  power  backscattered  ' to  laser  power  transmitted 


in  an  "average  ocean'* 


"5 

7.5  meter  (at  100  Hz) 
0.1  meter: 

R (meter) 


at  range  R for  a scattering  length  of 
and  a circular  receiving  aperture  of  diameter 
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Brownian  Motion  and  Brillouin  Effect 

a.  The  mean  distance  travelled  in  time  t by  a colloidal  particle  in  water  at  temperature 
(T)  = 300°K,  dynamic  viscosity  tj  - 10“3M“2s,  diameter  a of  particle  = 3 x 10~®M,  is 


- *10  /\ft  me+cr 


Thus  the  mean  distances  travelled  in  significant  times  due  to  Brownian  motion  are: 

(a)  in  one  period  of  acoustic  wave  at  100  Hz  0.7  x 10“^M 

(b)  in  one  period  of  laser  light  2.75  x 10“1^M 

(c)  in  time  to  interrogate  two  quarter  waves  of  laser 

wavelength  (2-way  travel  over  distance  equal  to 

a quarter  wavelength)  1.9  x lO^M 

(d)  in  time  for  acoustic  amplitude  to  go  from  zero  to 

maximum  (i.  e.  one  quarter  period)  0.33  x 10-^M 

(e)  the  time  to  cover  the  two-way  travel  between 

laser  and  particle  (2  x 30M  ■ 60M)  0.341  x 10“®M 

b.  When  observed  with  visible  light  the  Brillouin  eflect  in  water  is  equivalent  to  an 
acoustic  wave  at  approximately  10^®  Hz  modulating  the  laser  beam.  The  laser  wavelength  is 
then  associated  with  two  satellites  at  wavelength  separation  of  ~ 0.05A  (see  Appendix  E). 
Currently  Available  Laser  Systems 

Many  types  of  lasers  useful  for  heterodyne  experiments  are  available.  The  following  set 
of  specifications  can  be  considered  representative  of  more  recent  design  achievement. 
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SPECIFICATIONS* 


Brightness:  Model  K 1500  is  normally  furnished  with 

selected  ruby  laser  crystals  to  provide  a 
brightness  up  to  7.5  x 10^  watts/cm^/ 
steradian.  Through  special  selection 
techniques,  the  K 1500  can  be  supplied 
with  brightness  of  3 x 10*^  watts/cm 2/ 
steradian. 

Peak  Power:  1.1  Gigawatt  (Max.) 

Beam  Divergence:  Available  from  2.4  mr  FAHE  down  to 

1.2  millitadians 

Line  Width:  Less  than  0.06  A per  single  line  achieved 

spectral  component  with  optional 
accessory  (KLMS). 

pulse  Energy:  10-15  Joules 

Ruby  Life:  Depends  on  power  levels;  typically  300-400 

shots  at  1.1  gigawatt;  much  longer  at  lower 
power  levels. 

Beam  Size:  Approximately  1.9  x 1.76  cm  (elliptical). 

Can  be  corrected  to  round  shape  with 
optional  sapphire  prism. 

Repeatability:  .*10%  for  10  shot  series 

Wave  Length:  6943  A (Ruby) 

Q-switch:  Pockels  Cell 

Pulse  Width:  10-15  Nanoseconds 

Jitter:  ±10  Nanoseconds 

Repetition  Rate:  2 PPM 

Power  Supply:  10  KV 

Laser  Head:  4"  x 9/16  Ruby  Oscillator  and  9”  x 3/4" 

Ruby  Amplifier.  Water  cooled;  helical 
flash  lamp. 


♦KORAD  K1500  (2520  Colorado  Ave.,  Santa  Monica,  Calif.  90404). 
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Comment:  The  Summary  and  Conclusions  noted  above,  together  with  auxiliary  data 
on  Brownian  motion,  Brillouin  effect,  available  laser  systems,  etc.  were  presented  to  afford 
the  reader  a condensed  background  of  pertinent  facts  for  judging  the  merit  of  the  proposed 
use  of  laser  heterodyne  systems  to  measure  acoustic  signals.  A deeper  appreciation  of  the 
optical-acoustical  interaction  requires  a more  detailed  mathematical  model.  This  modeling 
is  presented  in  the  remainder  of  the  report.  The  reader  should  bear  in  mind  however  that 
while  the  results  of  this  investigation  are  based  on  an  intensive  effort  to  reach  the  key  issues 
involved  the  points  of  view  presented  do  not  exhaust  the  range  of  possibilities  in  regard  to 
both  analysis  and  experiment. 
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Introduction 


Present  day  techniques  for  the  detection  of  underwater  sound  are  based  on  the  local 
interaction  of  acoustic  pressure  with  a piezoactive  device  (vz.  hydrophone)  which  converts 
alternating  pressure  into  alternating  electric  current.  The  limit  of  detectability  is  circuit 
noise  (Johnson  noise).  Sound  is  also  made  visible  by  shadow  photography,  or  schlieren, 
in  which  the  local  change  in  illumination  depends  on  the  change  in  the  gradient  of  den- 
sity of  the  fluid  as  the  pressure  wave  passes  through.  Since  acoustic  pressure  (or  density) 
waves  are  scalars  the  determination  of  the  direction  of  the  wave  is  accomplished  by  arrays 
of  hydrophones  in  which  phase  differences  between  elements  are  used  to  maximize  sen- 
sitivity in  the  direction  of  the  approaching  wave.  These  arrays  must  be  several  wave- 
lengths long  to  be  effective,  but  length  can  be  reduced  by  use  of  acoustic  multipoles. 

Hydrophone  arrays  in  current  use  suffer  the  following  deficiencies.  (1)  at  low  fre- 
quencies an  array  several  wavelengths  in  size  becomes  very  long,  making  for  costly  struc- 
tures to  hold  them,  or  making  for  errors  due  to  lack  of  stability  in  arrays  which  are 
trailed  behind  ships  (towed  arrays).  (2)  all  sound  detection  is  local.  Hence  if  the  array 
is  placed  on  the  ship  hull  it  is  disturbed  very  strongly  by  presence  of  the  interfering  hull, 
by  flow  noise,  by  bubbles,  by  hull  motion,  etc.  If  towed,  the  long  array  suffers  from 
local  flow  noise,  from  bearing  ambiguities,  from  catenary  curvature  in  the  steady  drag 
condition  etc.  Other  deficiencies  in  hydrophone  arrays  axe  sensitivity  to  depth  of  opera- 
tion, lack  of  precision  calibration,  and  phase  and  amplitude  errors  due  to  non-uniform 
elements,  and/or  non-uniform  spacing. 

It  is  clearly  advantageous  to  devise  a system  of  sounu  detection  which  is  not  local. 

A promising  technique  is  to  use  an  optical  heterodyne  to  detect  the  phase  difference  be- 
tween a reference  laser  beam  and  a laser  beam  scattered  from  suspended  particles  in  mo- 
tion in  afluid.  The  feasibility  of  this  technique  is  investigated  in  the  body  of  this  report. 
Magnitude  of  Threshold  Underwater  Acoustic  Signal  to  Be  Detected 

In  anti-submarine  warfare  the  threshold  of  detection  of  a submarine  by  passive  lis- 
tening is  a somewhat  arbitrary  number.  It  will  be  taken  here  to  be  the  magnitude  of  a 
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signal  of  the  level  of  noise  in  the  ocean  at  sea  state  1 at  specified  frequency  in  a one  herz 
band.  For  purposes  of  this  analyses  we  choose  a frequency  of  100  Hz.  The  average  deep- 
water ambient  spectrum  level  of  plane  wave  equivalent  of  rms  noise  in  a 1 Hz  band  is 
found  in  several  references  (Horton,  Fundamentals  of  Sonar,  p.61,  Urick,  Principles  of 
Underwater  Sound  for  Engineers,  p.  168).  Two  values  of  noise  are  quoted,  one  without 
shipping  noise,  and  the  second  with  shipping  noise.  We  choose  here  the  smallest  magni- 
tude of  the  threshold  by  excluding  shipping  noise.  This  we  state  the  threshold  of  detec- 
tion in  terms  of  the  following  selected  plane  wave  equivalents: 

100  HZ  Plane  Wave  Equiv.  in  1 Hz  Band 

Spectrum  pressure  level  >70  dB  re  N/M2 

(or)  Particle  velocity  2.13  x 10**0  Ms** 

(or)  Power  spectrum  level  -132  dB  re  1 watt  M*2 

(or)  Particle  displacement  3,  4 x 10**0  ^ (.  3.4Xio*0 

Angstrom) 

The  addition  of  shipping  noise  increases  the  threshold  amplitude  by  approximately  20  dB. 
As  the  frequency  rises  the  level  of  ambient-noise  spectra  first  increases  (i.e.  the  ocean  be- 
comes more  noisy)  to  about  - 65  at  500  Hz,  then  declines  at  a rate  of  about  6 dB  per 
octave  until  it  reaches  the  thermal  noise  limit. 

Proposed  Method  of  Detection 

In  the  publication  “Virtual  Operture  Sonar”  R.M.  Ameigh  et  al,  (General  Electric  Co. 
973-SH-347-973-SO-1 42-01,  March  1970)  it  is  proposed  to  use  an  optical  Doppler  radar 
to  operate  underwater  for  the  measurement  of  the  motion  of  natural  colloidal-type  particles 
.suspended  in  sea  water.  Devices  of  this  type  have  the  generc  name  of  Laser  Doppler  Veloc- 
imeter  (LDV).  Since  LDV's  come  in  different  arrangements  it  will  be  useful  here  to  de- 
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In  this  system  a local  oscillator  beam  (or  "reference  beam")  couples  with  the  scattered 
light  from  a second  beam  at  the  focal  point  of  the  principal  lens  and  is  then  viewed  di- 
rectly by  the  photodetector  which  detects  the  Doppler  (Goldstein  and  Kried,  (1967), 

J.  Appl.  Mech  21*  813).  The  second  type  or  differential  heterodyne  arrangement  is  shown 
in  Fig.  2. 


AFERMte 


In  this  arrangement  the  light  scattered  from  the  common  focal  region  is  focused  by  a lens 
on  to  the  photodetector  which  detects  the  Doppler  shift  due  to  motion  in  the  focal  volume 
(Rudd,  1969  J.  Phys.  E.  2,  55).  The  third  type  is  the  symmetric  heterodyne  arrangement, 
Fig.  3 
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In  this  arrangement  the  scattered  light  bundle  is  stopped  by  an  aperture  with  two  slit 
openings,  which  are  combined  by  a lens  and  focused  on  the  photodetector. 

It  is  noted  that  in  all  three  arrangements  of  LDV  it  is  necessary  to  use  the  heterodyne 
technique  to  retain  phase  information.  This  means  a reference  signal  (of  coherent  light) 
must  be  mixed  with  the  scattered  light.  The  method  of  providing  the  reference  signal 
distinguishes  the  varieties  of  LDV. 

Unified  Mathematical  Modeling 

Two  beam  interferometry  is  briefly  reviewed  in  Appendix  A.  More  complicated 
devices  have  been  reviewed  by  Wang  (1972,  J.  Phys.  E.  5,  763).  He  has  proposed  unified 
mathematical  models  for  the  varieties  of  LDV  shown  above.  These  are  listed  here  for 
easy  reference. 

Type  I.  Local  oscillator  heterodyne 

A reference  beam,  amplitude  , phase  4*.  > frequency  W,  is  assumed  to  have 

the  form, 

EJt)  - a„  (-  iuxt  -t  i W) 

The  scattered  beam  is  a collection  of  light  scattered  by  N parti ces,  velocity  phase  , 
frequency  ^ , amplitude  , i.e. 
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. »jt  i i -jfc. 


f (*)« Z w|»  [-1^  + 4 1 ^fi 

?•• 

in  which  |(  = ‘(j  - being  the  scattered  wave  vector.  In  brief:  there  is  one  incident 
and  one  (generally  different)  scattered  direction.  The  power  spectral  density  of  the  photo- 
detector current  is, 

S (uj) - L ct,  (cl*  + 2 S ) 4 eV  ( al  4 1 oj ) &») 

in  • • 

+ eV'a  „*  1 4 j (<u  ± K-yf) 

l 

+ **  ZZL<4^^[co+K*  (yP- 

?*• 

in  which  G , <*  are  photodetector  parameters. 


Type  II  Differential  heterodyne 

Fig.  2 shows  that  there  are  two  incident  beams  with  wave  vectors  and  jfejs 
The  scattered  light  has  the  mathematical  form, 

»L  , 

Ev<*>=  L af  «>(»  (-iM,<){  *p  [i  ((ts-SaV*t 

In  brief:  there  are  two  incident  directions  and  one  (generally  different)  scattered  direction. 
The  power  spectral  density  of  the  photodetector  current  is, 

H.  ,yr  . 


SM 


A’f  a * ♦ <eV‘(M 

L — r i • 

7C  | ' 
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+ eV’  1 1 Mi  | £ [<ui (K-kri)-  (&-$)] 

r* ' 1 K 

* i 4(«  -*  fe.-^  i fc  - *»)•  & -Ji )] 


Type  III  Symmetric  Heterodyne 


It>  the  symmetric  heterodyne  the  scattered  light  Held  is  given  by 

N. 

*»f  [i  +iU) 


In  brief:  there  is  qq£  incident  direction  and  two  scattered  directions.  The  power  spectral 
density  is  quite  complex,  but  it  can  be  retrieved  in  the  cats  of  common  particle  velocity  by 
replacing!  ^zZ  - (jj^"  ^«i)in  the  above  formula  for  Type  II. 

Concluding  Remarks:  The  power  spectra  formulas  derived  by  Wang  (loc.  cit)  are  spectral 
lines,  the  line  widths  not  being  considered.  These  spectra  serve  to  identify  th*  prominent 
features  of  the  Doppler  - induced  frequency  shifts.  However,  by  failing  to  present  line- 
widths  the  spectra  do  not  ixsist  us  in  estimating  the  detectability  of  these  lines  in  the 
presence  of  noise  originating  both  in  the  fluid  and  in  the  LDV  itself.  Since  the  inherent 
noises  in  the  systems  will  te  our  chief  concern  we  will  be  compelled  to  find  and  develop 
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more  detailed  models.  A list  of  the  important  noises  in  LDV  systems  is  given  here. 

Noises  in  LDV  Systems 

The  linewidths  of  the  Doppler-shifted  spectrum  originate  from  the  following 

causes: 

1.  laser  light  non-chromaticity. 

2.  Brownian  motion  of  the  colloidal  particles. 

3.  velocity  gradients  across  the  scattering  volume. 

4.  fluctuations  of  velocity  in  the  scattering  volume  (i.e.  turbulence). 

5.  angular  uncertainties  due  to  the  divergence  of  the  incident  beam  and  detector 
angular  aperture. 

6.  finite  passage  of  scattering  particles  through  the  laser  beam. 

7.  temporal  jitter  of  the  electronic  LDV  apparatus. 

8.  temporal  and  spatial  changes  in  the  index  of  refraction  of  the  scattering  fluid. 

In  order  to  account  for  such  a profusion  of  quantities  which  broaden  the  Doppler-spectrum 
we  must  find  their  order  of  appearance  in  the  more  detailed  model  cited  above.  Since 
the  model  to  be  developed  is  very  complex  it  will  be  useful  to  first  present  a physical 
picture  of  a generic  LDV,  following  which  we  will  develop  the  math  model  with  greater 
understanding. 
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Simplified  Model  of  the  Laser  Heterodyne 
Detector  of  Sound 


Simple  Model  of  a LDV 

The  optical  heterodyne  (or  Laser  Doppler  Velodmeter  (LDV))  can  be  modeled  in  a 
simple  way  (Rudd,  J.  Sd.  Inst.  Vol.  2 Series  2 1969,  p.55),  Fig.  A 


MA'XA 
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Here,  a laser  beam  is  first  expc  nded  (Station  1)  then  focused  by  a transmission  lens  in 
the  fluid  (2)  then  collected  by  a receiving  lens  (3)  on  to  a photodetector.  In  order  to 
produce  “beats"  (i.e.  a fringe  system)  the  beam  from  (2)  is  masked  (MASKA)  by  a 
screen  with  two  splits  before  it  is  brought  to  a focus  in  the  fluid.  Alternatively  tho 
beam  is  masked  at  position  3 by  MASK  B,  after  it  leaves  the  focus  in  the  fluid.  Mask 
A generates  a real  set  of  fringes.  Mask  B generates  a virtual  set  of  fringes. 

The  two  positions  of  the  masks  correspond  to  two  systems  of  LDV  in  current  use. 

A system  with  Mask  A describes  the  Goldstein-Kreid  experiment  (J.  Appl.  Mech.  1967 
34  813-8).  In  this  experiment  a real  set  of  fringes  is  generated  at  the  focal  volume  on 
the  Quid  by  beam  splitting  the  laser  light,  then  bringing  the  two  beams  to  a focus  at  an 
arrival  angle  between  them.  Scattering  then  occurs,  which  is  then  detected  (as  in 
Fig.  1). 

A system  with  Mask  3 describes  the  Yeh-Cummins  experiment  (Appl.  Phys.  Letters 
4 176-8  1964).  Here  scattering  occurs  Irst.  The  scattered  light  is  then  brought  together 
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with  the  unscattered  (or  reference)  beam  (at  angle  24)  to  form  a virtual  fringe  system 
which  is  then  detected  by  the  photodetector. 

Simple  Description  of  Photocurrents 

When  MASK  A is  the  position  of  the  effective  screen,  and  a set  of  real  fringes  is 
available,  a scattering  particle  entering  the  focal  region  and  crossing  the  fringes  sampies  the 
succession  of  light  and  dark  bands.  The  scattered  light  is  alternately  bright  and  dark,  with 
maximum  brightness  occurring  in  the  center  of  the  fringe  system.  The  two  scattered  laser 
beams  , E are  modulated  by  the  fluid  flow  (at  velocity  ),  and  are  multiplied  by 
the  photomultiplier  to  give  an  electric  current  proportional  to  their  product,  i.e. 

T OG  + h*  + Es,  Esv  CBS  ) 
i *• 

in  which  Ys  « the  Doppler  shift  in  the  laser  frequency  due  to  a (constant)  velocity  of  the 
scattering  particle.  A sketch  of  this  current  is  shown  in  Pig.  S 


When  MASK  l is  in  position  a set  of  virtual  fringes  is  found  in  the  focal  volume.  A 
particle  crossing  this  volume  at  constant  velocity  V alternately  scatters  and  does  not  scatter 
light  as  it  is  illuminated  with  the  light  and  dark  bands  of  the  fringe).  The  scattered  light 
is  combined  with  a reference  beam  Eg(i.e.  the  second  beam  generated  by  MASK  B) 
which  is  detected  v a current  X where 

IoC  V Ik  * c*v(m  ) 

A sketch  of  this  current  is  shown  here, 
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The  Doppler  Shift 

To  understand  the  Doppler  shift  (which  is  the  same  in  both  experiments)  we  first 
consider  the  focal  volume  to  be  infinite,  i.e.  we  first  consider  an  infinite  set  of  fringes  to 
be  available,  each  of  equal  intensity.  A scattering  parti  i crossing  there  fringes  at  right 
angles  (to  the  fringe  lines)  samples  the  alternating  light  and  darkness  (as  noted  earlier). 

This  has  the  same  effect  as  making  the  particle  stationary  and  allowing  the  set  of  fringes 
to  slide  over  the  particle.  Let  the  speed  of  the  fringe  motion  be  . If  the  wavenumber 
of  the  fringes  is  **^^then  the  radian  frequency  of  appearance  of  (say)  the  brightest 
scattered  light  is  Kefy  (radians  per  sec),  that  is,  the  time  required  to  ride  over  one  wave- 
length is  ^e/^sT  ( in  units  of  seconds/cycle),  so  that  the  frequency  of  alternating 
darkness  and  light  is  T = which  is  a "Doppler  Shift."  Since  light  frequency  -Q. 
is  given  by  we  also  see  that  fcJ,*  -Q«  {^/cg  ) , Or  $«  (*A) 

We  next  take  the  focal  volume  to  be  finite,  that  is,  we  consider  the  set  of  fringes  to 
be  bounded  in  space.  The  time  history  of  the  scattering  of  light  is  then  as  follows:  zero 
before  the  (finite)  set  of  fringes  rides  over  the  particle;  abruptly  oscillatory  when  the 
fringes  cross  the  particle;  abruptly  zero  again  at  the  termination  of  the  fringe  set  The 
frequency  of  appearance  of  the  light  and  dark  scattered  light  must  be  measured  from  the 
entire  time  record  of  -o^ta-MP.  This  record  can  be  visualized  as  a infinite  sinusoid 
screened  by  a black  screen  with  a single  rectangular  window.  The  time  history  of 
scattered  light  is  thus  a product  of  two  functions  (sinusoid  times  window).  The  Fourier 
spectrum  of  this  time  history  is  the  convolution  of  the  spectrum  of  each  function  separately, 
i.e.  the  convolution  of  a delta  function  5 (SI-  Ht)  and  a siny/y.  The  resultant  spectrum 
(in  the  range  (o*o  is  a sinX/^  centered  at .Qg. Thus  the  spectral  representa- 

tion of  the  time  history  of  the  crossing  of  a particle  through  a finite  set  of  fringes  has  a 
spectral  width  of  the  first  lobe  of  the  sinX/(  spectrum.  If  is  the  (finite)  time  of  ^ 
transit  sf  the  particle  then  the  bandwidth  of  the  first  lobe  of  sinX^f  t*  (A^) 

(Brace  well  “The  Fourier  Transform"  p.  128,  and  p.  368).  Thus  the  spread  in  the  Doppler, 
A^, , due  to  finite  transit  time  At,  is  A ■Jjs  (fit,)  This  spread  is  essentially  a repre- 
sentation of  uncertainty  in  the  Doppler  (see  Figs.5 , 6). 

A second  approach  to  the  calculation  of  Doppler  spread  due  to  finite  transit  time  is 
to  use  the  model  given  by  Fig.  1,  and  consider  the  masks  to  be  two  slits  of  the 
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equivalent  interferometer.  The  aperture  of  the  equivalent  slits  is  the  spatial  function 

'(«•  n(^Vn(!*) 

in  which  (1  is  a rectangle  function  of  unit  height,  centered  at  X'-tV  , and  Vf  units  wide. 
The  Fourier  transform  of  this  aperture  is, 


■PU')-  «(«!*) 


(1.4) 


(Bracewell,  p.  283).  The  power  spectrum  is  (approximately)  the  square  of  this,  or 


s(rt-  •(*)'( 


ft 


VSt/\ 


coi 


(1.6) 


If  we  set  the  variable  >/^r  "t}  (^-  “ interpreted  as  time)  then  the  “width  of  the  inter- 
ferometer” is  calculated  in  seconds,  A't  tec.  and  the  first  zero  of  S(f)occun  at 

A*  $ *1)  $'(4^)  (“  before).  We  thus  model  the  spectrum  of  the  photodetector 

current  by  S(-f)* 

Basic  Coherence  Requirements 

In  the  model  of  Fig.  4 the  basic  optical  system  is  that  of  a two  beam  interferometer. 
A review  of  this  system  will  help  define  coherence  requirements.  Let  E,  (-  4^  ylbe 
the  comple^amplitude  of  the  first  beam  at  distance  5,  from  the  two-slit  screen,  and  let 
E (=  jbt  the  second  complex  r iplitude  observed  at  the  same  observation  point 

but  traveling  distance  iS*  The  total  intensity  of  light  «Tu 

1=  EE*  = V-A,*  ->  2A,A,cm (it ) 

in  which  Jl*  is  the  optical  path  difference.  Since  cos  («)  is  periodic  we  define 

an  order  number  (K.D.  Mielenz,  p.  89  "Electron  Beam  and  Laser  Technology”  ed.  L. 

Marton,  1968),  m - 41  s 0$€<l  ^integral. 

Ill  „ * 

k 

This  means  that  the  phase  kl  hu  an  uncertainty  in  order  number  Am)  due  to  an 
uncertainty  of  path  difference  caused  by  finite  extension  of  the  source  (i.e.  of  the  slit) 
and  to  an  uncertainty  in  wavenumber  dfe  due  to  lack  of  monochromaticity.  Assume 


Now,  to  quantities  of  small  A It 
A 'Wi » 
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next  that  there  is  an  incremental  uncertainty  in  ^ i.e. 

U * 2Tt<W  + ft  7E  lAml 

(I,8) 

if  (4wi=  % , it  is  seen  that  the  change  in  intensity  aI  due  to  uncertainty  in/w  is 

- A -t-  Ax  - 2 Aj 

Thus  the  uncertainty  in  intensity  cancels  the  cosine  term  of  the  intensity  itself  leading  to 
the  conclusion  that  the  condition  (ftW  |=  V2  results  is  incoherent  superposition  of  the 
two  light  beams,  with  the  disappearance  of  observable  interference  effects.  Thus  the  limit 
of  coherence  is  given  by  * 

(M>)  =±1t  as. 

One  can  arbitrarily  define  a coherence  condition  if  (A  m | 1/4,  or  if, 


(4  if  )*+ 


The  two  parameters,  4A,  A It  , require  separate  investigation. 

Case  I.  Assume  source  (i.e.  slit)  is  infinitesimaly  small  so  that  /is  exactly  specifiable, 

[i't.  A$-o)  and  assume  that  the  source  emits  a spectrum  of  finite  width  i(nlt)cf. 

The  uncertainty  in  the  order  number  is  den  where 

t»*±l  c ii  (*,-<■)  ,n„ 

U JCt  *cs 

The  temporal  coherence  length  ol  the  source  is  thus  found  from  the  condition  £ A It  = i 1C 

c«K 

to  be 


Thus  the  greatest  path  difference  over  which  fringes  can  be  observed  is  controlled  by  the 
frequency  spread  of  the  source.  When  tl  t source  is  white  light,  <*)*-&>,  is  very  large, 
and  the  temporal  coherence  length  is  small,  i.e.  fringes  are  observed  only  over  very  small 
path  distances.  In  contrast  laser  light  has  a value  ~ j > IQ**"*  Hence 

~ ill ^ 3y  IQ*1  cm  or  30  million  kilometers. 

In  a multimode  laser  the  Doppler  width  is  much  greater,  e.g.  a neon  line  is 
)(q  ^ 0.3  rad/cm.  Since  the  mode  separation  in  a laser  of  length  L 1 4 
> if  i» ««n  that  the  number  of  possible  modes  is 


M (“‘3 1 _ L 
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For  a laser  length  L s 100  Cwv  } 0.03  rad/cm.  Thus  10  modes  may  be  in  operation. 

The  coherence  length  with  10  modes  is  thus  reduced  to 

^ 20  orrw 

o i « 

Case  II.  Spatial  Coherence 

Assume  one  (slit)  source  has  a finite  (width)  dimension  (say  QQ).  and  emits  perfect 

monochromatic  light  ( <&*  x (i)0  ~ ^ ).  Let  ot  be  the  angle  at  observation  point  ? which 

encloses  QQ  (Fig.  7*).  From  Q there  is  emitted  a wave  which  arrives  at?  simultaneously  with 

a wave  from  the  corresponding  $ of  the  second  slit,  with  a path  difference  X The  geometric 

picture  of  two  wavefronts  at?  can  be  duplicated  at  Q by  assuming  P to  issue  two  wave- 

fronts  which  intersect  at  (J(of  one  slit)  with  the  same  path  difference  /.  The  wavefronts 

w,  > ^ also  intersect  in  the  plane  of  the  figure  along  a line.  Assume  point ^ is  on  this 

line  of  intersection.  A is  chosen  such  that  one  wavefront  passes  through  Q of  the  slit. 

I 

The^  second  wavefront  is  at  a path  c .fference  (to  Q ) given  by  QS.  If  the  slit  dimension 
QQ  is  projected  on  the  line  of  intersection  to  give  dimension^  , then  the  path  difference 
os  is  (approximately)  . Thus  the  effect  of  finite  slit  size  is  to  introduce  a path  length 
uncertainty,  from  zero  at  Q then  along  all  intermediate  points  to  -igd  at<?  . Using  the 
center  of  the  slit  as  reference  the  uncertainty  is  zero  to  - $<*(l  which  is  of  the  above 
formula  Eq.  1.7.  The  source  is  spatially  coherent  if 


« £ 


or 


(1.14) 


The  largest  source  (width)  with  which  interference  fringes  can  be  observed  is  called 
which  is  twice  the  above  number,  or 

W. « ^ k (U6> 

Spatial  coherence  is  thus  controlled  by  slit  width  and  wavelength  of  the  laser  light. 
Statistical  Nature  of  Coherence 

Interference  by  n single  pair  of  wavetrains  is  too  rapid  for  any  detector  to  follow. 
Under  actual  experiment  the  set  of  fringes  generated  by  the  equivalent  mask  is  the 
statistical  average  of  a large  number  of  wavetrains.  The  coherence  time  of  a typical  gas 
laser  is  of  the  order  of  10“ 2 sec.  Although  this  allows  transient  tracking  by  a detector  of 


F»C,  rl 
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single  pair  interference  fringes,  the  time  required  for  the  remainder  of  the  experiment  and 
signal  processing  only  allows  statistical  averages. 

For  ideal  two  beam  interference  the  intensity  at  the  observation  point  is  (as  noted 
above) 

I--  A*  ■>  L*  xA,K  *4*) 

(1.16) 

and  the  fringe  visibility  is 


A A ^ . 

V*  A 


For  a general  two  beam  interference  experiment  this  is  replaced  by  a statistical  statement, 


I = I , + Iz  + 2 J I,  L I Vu  (r)  | [<*»  ft)]  (1.18) 

in  which  is  the  complex  degree  of  mutual  coherence  of  two  points  St>  in  the 

wavefield,  |'9,I  | ( Of  (li  are  its  magnitude  and  angle,  and  'IT  is  a time  shift  between  the 
signals  fromS|,an4  5^,  The  angle  It)  is  obtained  from  the  phase  condition  for  the 
appearance  of  bright  fringes, 

<yu  (t) » a.7i  N,  N«  in-teg&r  (L19) 

Hie  magnitude  (aln  | is  obtained  from  observation  of  the  two  beam  interferometer 


measurement  of  visibility  at"^  called 


t 


(t)  1 * Mtnea*  a 

\t  I was  ~ Jrtm 

It4A*  +Ihh 


yrxet , ” 


(1.20) 


(1.21) 


Degradation  of  coherence  is  caused  by  several  factors:  (1)  increase  in  path  length  (2) 
change  in  particle  numbers  and  particle  density  in  the  Quid  (3)  change  in  index  of  refrac- 
tion of  the  Quid.  A typical  experiment  to  test  changes  in  coherence  is  to  generate  a set 
of  interference  fringes  in  space,  and  observe  the  distortion  of  these  fringes  due  to  the 
above  factors.  Hie  degradation  of  coherence  in  a two-beam  interferometer  experiment 
using  an  argon  laser,  <5.  rAM.  beam,  5145  A0  has  been  studied  by  R.  E.  Lee  et  al  (Opto- 
electronics) 5.(1973)  41-51.  The*e  authors  used  distilled  water,  coastal  water  and 
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fabricated  salt  water  as  transmitting  media.  All  water  samples  showed  approximately  the 
same  value  of  normalized  degree  of  coherence  (about  0.97),  independent  of  their  nature 
and  independent  of  path  length  (for  paths  up  to  3 meters).  Agitating  the  water  mechanically 
by  moving  a stick  back  and  forth  along  the  water  trough  showed  no  change  in  fringe  pattern, 
except  in  the  fabricated  salt  water  where  a transient  reaction  (i.e.  spotty  unsteady  motion 
of  the  fringes)  was  observed  due  to  a transient  non-homogeneous  condition  of  the  water- 
salt  mixture.  This  transient  disappeared  after  a few  moments.  The  interference  patterns 
generated  in  space  were  recorded  by  a traversing  photomultiplier.  These  records  show 
(a)  relatively  smooth  oscillations  of  light  intensity  due  to  fringes  for  tap  water  and  distilled 
water  (b)  ragged  or  spiky  oscillation  of  light  intensity  for  coastal  water  and  fabricated 
salt  water,  both  containing  suspended  particles.  The  spikes  appear  to  correspond  to  the 
observed  low  frequency  local  wandering  of  the  fringes.  It  is  concluded  that  any  inhomo- 
geneity of  the  medium  which  disturbs  the  phase  uniformity  of  the  laser  beams  will  decrease 
the  visibility  of  the  fringes  (i.e.  decrease  coherence).  As  noted  earlier,  these  disturbers  are 
(1)  suspended  particles  which  introduce  (moving)  delta  function-type  density  changes  (2) 
temperature  changes  which  also  introduce  mass  density  changes  in  the  medium.  When 
these  inhomogenities  are  functions  of  time  they  cause  tl.e  fringe  system  to  undergo  wriggly, 
sloshing  motion,  with  transient  brightening  and  fading.  The  frequency  spectrum  of  the 
photomultiplier  current  shows  considerable  broadening  ("noise”,  "Doppler  noise”). 

Coherence  Requirements  for  LDV 

The  temporal  coherence  requirements  for  a two-beam  interferometer  noted  above  are 
somewhat  modified  when  fluid  motion  is  specifically  to  be  measured.  If  the  spread  of 
wavelengths  in  the  source  is  X*  and  the  spread  of  the  Doppler  is  -J,  - fi  then  one 
requires  that 

— > X©  , So  - ^o-Str 

ve«dengflv,^«y  (1.22) 

(Rudd,  loc.  at  eq.  (13)).  Similarly,  the  spatial  coherence  requirement  for  an  LDV  is  given 
by  specifying  the  width  d of  the  equivalent  slit  in  the  two-beam  interferometer, 

i-r-A 

IS** cc  \ I (1.23) 

(Rudd,  eq.  (16)).  The  definition  of  of  is  the  same  as  in  section  on  Basic  Coherence 
Requirements. 
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A Phasor  Description  of  the  Received  Signal  in  an  LDV 

A single  particle  crossing  the  (finite)  fringe  pattern  at  constant  velocity  samples  the 
pattern  and  generates  the  (scattered)  wave  packet  shown  in  Fig.  (6).  If  its  motion  is 
random  it  will  generate  pieces  of  wavepackets  in  a random  manner,  thus  introducing 
random  noise  whose  spectral  bandwidth  is  inversely  proportional  to  the  time  duration  of 
the  packet.  A time  record  of  random  particle  motion  in  an  LDV  will  be  a series  of  spikes. 

If  many  particles  are  simultaneously  crossing  the  fringe  pattern,  each  scattered  wave- 
packet  will  have  different  phase.  The  net  amplitude  of  the  wavefront  scattered  by  all  the 
particles  will  be  proportional  to  the  square  root  of  their  number. 

If  the  particle  samples  the  fringe  pattern  periodically  (frequency  6^ ) it  will  generate  a fre- 
quency modulation  of  the  original  carrier  wave  (frequency Xl£)  that  generated  the  fringe 
pattern.  The  net  time  function  of  the  scattered  wave  (for  particle  velocity  v,  acoustic  modu- 
lation mFm)  is-  Es  * A a*  t 1 { Su* + n s,n  + --i 


Of 


E*  = A, 


^ Sin  &>**  + 


(1.24) 


in  which  is  the  modulation  index  (s  K"*  defined  below),  &r»4^|=  ^/C,. 

If  the  brightness  pattern  is  subject  to  fading  then  amplitude  A will  be  a function  of 
time  (~  Alt))  so  that  the  scattered  wave  will  be  amplitude  modulated. 

If  N particles  participate  in  the  scattering  process  during  an  observation  time  the  time 
history  of  the  scattered  w ive  will  be 

a-r  (i-25) 

in  which  "P  particles  move  with  a signal  (frequency  ) and  fK?  move  randomly.  If 
particles  all  move  with  same  signal,  but  have  random  components  of  motion  as  well,  then 

E.W.  Z<a 

* (1.26) 

This  model  says  that  the  n’th  particle  will  undergo  k.  random  steps  (the  exact  number 
depending  on^),  scattering  a spike  at  each  step,  while  simultaneously  scattering  a 
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Scattered  Field  In  terms  of  the  Simplified  Model  (Edwards  et  al  J.  Appl.  Phys.  42,  837  (1971)). 

Let  the  finite  scattering  volume  A J have  H particles  at  time  tsO  . The  distance 
moved  by  the  /ftttiparticle  is  f0<*  (“0*  (o)  + (i)  + Vi  4 k 

The  scattered  laser  field  may  be  modeled  as  follows: 

E A Z ^ t «i.  <•)+  X*  ) 

~ m & “ 

, „ t [ ( K-K>(/w (.)  + 


(The  significance  of  the  terms  appearing  in  this  math  model  are  detailed  in  Appendix  H) 
Here  we  desire  to  describe  the  scattering  process  in  terms  of  the  simple  model  of  particles 
sampling  fringes.  The  following  points  are  to  be  noted: 

(1)  the  initial  phase  of  £ s is  exp  c (K- 

(2)  each  random  walk  step  samples  a fraction  of  the  fringes  and 

introduces  spikes,  plus  wavering  noise  in  the  scattered  light  associated  with 
phase  changes  exp  1 ( K — ) • h £0*  (t ) . 

(3)  the  steady  velocity  V samples  the  entire  fringe  set  and  generates  a complete 
wave  packet  of  the  sinusoid  cos  (anfrt)  multiplied  by  the  window 
function  defined  by  the  scattering  volume. 

(4)  the  sinusoidal  particle  displacement  ^ t samples  the  fringes 

periodically.  The  number  of  fringes  sampled  depends  on  the  amplitude 
and  number  of  cycles  . If  is  small  then  the  sinusoidal  sampling 
appears  as  a (nearly  constant) velocity  additive  to  V ; 1’  f.  the 

apparent  Doppler  is  (K “ K ) ' ( Y If  many  cycles  occur 


during  the  observation  time  then  the  number  of  fringes  sampled  depend  on  ^ . To  clarify 
the  physical  picture  we  assume  the  steady  velocity  V is  zero.  At  time  "t  and  position 
/v*4>  (o)  inside  the  fringe  set,  the  amplitude  of  particle  displacement  is  [/tc+  (°)J  • 
Let  fringes  ( jl  not  generally  integer).  If  ^ is  less  than  one  fringe 

(namely,  less  than  a half  wavelength  of  light  at  frequency  ) the  scattered  electric 
field  exhibits  little  phase  modulation  (it  is  black,  gray  or  white,  almost  independent  of 
time).  If  4 is  exactly  two  neighboring  fringes,  the  scattered  field  is  nearly  sinusoidal, 
with  equal  plus  and  minus  amplitudes.  If  ^ is  more  than  two  fringes,  plus  a fraction, 
the  scattered  electric  field  is  periodic,  but  not  sinusoidal,  see  Fig.  (9)  ^ 


In  which  the  period  of  sampling  is  T * a%)4 . The  modulation  index  is  (K*  K ) 
•£[$«)]  in  which  The  Jacobi  expansion  of  exp ([ Msin&j  is  found  in 

Watson  "Bessel  Functions.”  page  22), 


oo 


z - sill 


/»*! 


nr*  o 


in  which  is  the  *nlA  order  Bessel  function  (of  the  first  kind),  and  8 1 t,  For 

very  small  M 


% 3©(fl)  +at  Ti( *\)pun.&si  % 


1+  L Y\/tUvi  0>,t 


Thus  for  very  small  modulation  index  the  phase  moduli  tion  appears  as  a (sinusoidal) 
amplitude  modulation.  Considering  only  a partial  product  we  have 

^ (cai^|ct) -t/VA.Xlef  )(l+ t M 

= + M 


The  spectrum  of  the  modulation  is  given  by 

SM*  H j S(uhIU  + U,s)  4 l ((Al&E-tii) 

-tS(0O-IU-  Ct\ 

In  the  positive  half  of  the  frequency  spectrum  the  modulation  appears  as  two  side  bands, 
at  a distance  t Brace  well  "Fourier  Transforms”  p.  184).  Thus  the  spectrum  of 

scatter.-j  the  light  amplitude  is  typically  that  of  a sinusoidal  FM  (or  AM)  modulated  carrier. 
However,  these  spectral  lines  have  no  finite  width,  and  thus  do  not  account  for  noise. 
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Conclusion  to  Part  I 


In  Part  I we  have  presented  a physically  appealing  model  of  the  laser  Doppler  veloci- 
meter  in  terms  of  a set  of  laser  interference  fringes  crossed  by  occulting  particles.  This 
model  is  useful  in  providing  insight  into  the  scattering  process,  thus  allowing  experimen- 
talists to  master  the  details  of  the  phenomena  in  their  laboratory  experiments.  It  is 
however  very  important  to  comprehend  the  exact  nature  of  what  the  experiments  are 
measuring,  and  to  construct  a quantative  model  of  the  complex  nature  of  the  “noises” 
in  the  experiment  which  set  lower  limits  in  detectability  of  the  fluid  motion.  A very 
detailed  math  model  has  been  devised  by  George  and  Lumley  (J.  Fluid  Mech.  60,  p.  321 
(1973)).  We  adopt  this  model  in  all  succeeding  discussions. 
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Part  II 

Mathematical  Modeling 


The  scattering  volume  is  assumed  to  be  three-dimensional  (effective  dimensions 
,0^)  and  to  contain  more  than  one  scattering  particle  in  motion.  A first  step  in 
math  modeling  is  to  project  the  volume  on  to  the  *j|  plane,  Fig.  (lo).  Next  the  random 
velocities  of  the  light  scattering  particles  are  averaged  over  particle  position  on  the  scat- 
tering volume  to  give  the  instantaneous  equivalent  of  an  ensemble  averaged  (random) 
velocity  U»  (-i)  * t V«  ft))  (which  is  the  Eulerian  velocity  reported  by  the  detection 
system  as  producing  the  Doppler  shift),  and  a (random)  deviation  from  this  average, 

A l)  (*,■()  *(dU.  (*,i)  which  the  particle  at£  possesses.  The  displacement  of  a 

particular  particle  (at  X ) in  timet  is  thus  assumed  to  be  the  sum  of  a displacement 
3)o(t)r  ( X(*)  (the  ensemble  average  for  all  particles)  due  to  (<)  and  ( >,t)* 
du®  unique  to  the  particle  at  X. 

' » f ‘ V>^  , 


\ / \N 


" f*io  * 

The  analysis  proceeds  by  restricting  attention  to  a single  dimension,  (say  x-dimension). 
The  random  Eulerian  velocity  i^(<)  is  separated  into  a random  mean  velocity  Us  ft)  and 
a fluctuation  about  the  random  mean  A time  average  of  14ft)  is 

labeled  Uo  The  square  of  the  random  deviation  A U j^X  (•  ,t)  +|  is  ensemble  averaged  over 
particle  position  t and  time  averaged  over  ^ to  give  the  average  mean  - square  velocity 
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deviation  (aur  When  the  motion  is  laminar  and  deterministic  , 

U0-'H  i U.'J  = 0 = ^u)1 

The  time-averaged  random  mean  velocity  (in  general  turbulent  flow)  is  assumed  to 

contribute  the  Dopp'er  spectral  line  in  the  photodetector  (electric)  current,  The  Doppler 

shift  of  this  line  is  U3cz  Km*  where  l\  is  the  x-component  of  the  “vector  wavenumber”. 

The  fluctuating  quantities  Uo^  Alt  are  assumed  to  contribute  spectral  noise  (called 

Doppler  broadening)  centered  on  the  spectral  line  at  <k>Q 

/ 

Assuming  the  fluctuations  ti  , ALL  to  be  Gaussian  distributed  one  can  show  by  math 
modeling  that  the  frequency  bandwidth  which  measures  the  Doppler  broadening 
(«*<•*#)  is, 

(A*Bf  = \C  + K fciu?  + ± 

N ■>  ir-  a 


10? 


Here  the  term  IL  /atf)1  is  the  spectral  broadening  due  to  finite  transit  time  of  the  ‘average’ 
particle  traveling  at  averaged  speed  U,  through  the  effective  (projected)  volume  . In  the 
absence  of  fluctuations  the  broadening  of  the  spectrum  due  to  finite  travel  time  is  called 
the  Doppler  ambiguity  of  RADAR. 

When  other  sources  of  time  fluctuation  of  the  scattering  particles  are  present  they 
can  be  assumed  to  be  Gaussian  distributed.  The  most  important  of  these  contribute  band- 

a 

widths  as  follows:  A A)  due  to  gradients  in  the  (laminar)  mean  velocity  U across  the  volume, 

•>  $ a 

A(ifc  due  to  Brownian  motion,  and  A4)s  due  to  non-monochromaticity  of  the  laser  light 
beam.  Assuming  these  to  be  important  one  rites  the  frequency  bandwidth  of  spectrum 


to  be 


p*)*.  kV  * t?*)2 

(ph)%-  K (4U)*  + t iWj1 


Here  DA  is  called  the  Doppler  ambiguity  of  the  model. 
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Two  important  ratios  can  be  formed  from  these  assumptions  and  definitions.  These 
are:  (1)  Wa  the  ratio  of  the  Doppler  ambiguity  to  the  Doppler  shift,  (2) 
or  the  ratio  of  fluctuation  inU.to  the  laminar  W which  is  equivalent  to  the  ratio  of 
spectral  broadening  due  to  Uto  the  line  spectrum  due  to  These  ratios  are  useful  in 
setting  limits  to  resolution  in  applications  of  the  LDV  to  turbulence  measurements. 

The  basic  model  of  the  LDV  is  the  equivalent  two-beam  interferometer  using  waves 
of  laser  light.  The  Doppler  (electric)  current  t ^reported  by  the  photodetector  of  the 
interference  light  intensity  is 

i(t)~  F (t) c&i KX  (q (t W* !0( = (F+  Q [KX (*) ~ 

in  which  F(-t)^ft^  account  for  fluctuations  in  the  scattering  process  a~d  finite  beam 
widths  of  the  laser  beams,  and  0)  is  a phase  (of  light  intensity)  given  by  tan  ’Wi)  .Thus 
the  photodetector  current  is  ....  presented  by  a phasor  whose  amplitude  and  pnase  are 
random  functions  of  time.  Tne  probability  distribution  of  F and  Gj  are  assumed  Gaussian. 
The  fluctuation  ^4>/cttof  the  phase  is  important  in  the  determination  of  the  Doppler 
spectrum  since  the  equivalent  Doppler  frequency  is  d»(  = - &ft)J  ■ The  corre- 

lation properties  of  <Pft)are  modeled  as  filtered  noise  m an  FM  receiver.  Its  power  spec- 
trum is  labeled  N(o/V  It  has  been  found  (by  experiment)  to  be  white  with  a magnitude 
proportional  to  the  Doppler  ambiguity  DA(i.e.  N ^ N(o)  * 0.  J4  2 3>A  ) « 

Under  the  assumption  that  the  scattering  volume  is  larger  than  the  wavelength  of  the 
laser  light  ( i e K IP* >>  I ^ it  can  be  shown  by  math  modeling  that  the  correlation  of 
F(*)  is  approximately  the  same  as  the  correlation  for  § that  the  cross-conrelation  of 
F and  ^ is  approximately  zero,  and  that  the  time  average  over  fluctuations  (but  not  the  mean)  is, 


K [X(*>  - f 


Here, 


and  X^)i»  the  time  average  (or  mean)  of  the  random  displacement  X^).  Similarly,  under 
the  same  assumption  of  volume  size,  the  math  model  shows  that 
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r«)  go 

F‘(t) 


* j- 


kV 


} 


) 


where 


rS  a (^y.r  (t-t') 


Thus  the  normalized  correlation  of  the  random  photodetector  current  is  modeled  as  the 
correlation  of  a phasor  with  random  amplitude  and  phase, 

. i„b  (_  kV4  X(t)-XftO)  e«,  L M 

^ 4<T|*  J ^ ^ 

• CeiK  [X(t)-X«'j] 

in  which  the  exponentials  account  for  the  fluctuations  in  both  amplitude  and  phase. 

The  basic  math  model  is  first  used  to  calculate  the  power  spectrum  of  the  photo- 
detector current.  In  this  modeling  we  require  an  explicit  form  for  the  average  (time- 
varying)  displacement  X^)  to  be  attributed  to  aH  the  particles.  We  assume  that  all  the 
particles  have  a volume  averaged  velocity  tip  and  an  acoustic  displacement  Xv  sin  «stwith 
this  choice  of  it  is  found  by  Fourier  transformation  of  the  autocorrelation  of  t (t) 
that  the  Doppler  spectrum  consists  of  two  Gaussian  peaks  at  <fc)  - 1 K^t|0with  amplitude 
proportional  to  S,(K«0  each  peak  in  turn  associated  with  two  satellite  Gaussian  peaks 
^(e)r  ±(t^  with  amplitude  proportional  to  T,  making  a total  of  two  groups  of  peaks 
with  three  peaks  in  each  group.  The  "spectral  lines”  in  one  group  are  thus:  Kilo 
KU#  - (0S , The  math  model  also  predicts  that  each  spectral  line  will  be  broadened  by 
fluctuations  and  finite  transit  time.  The  Doppler  broadening  will  be  48  as  defined  above, 

SB 
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This  model  will  be  now  be  used  to  determine  limits  of  detection  of  the  acoustic  signal  in 
the  presence  of  noise,  and  to  estimate  the  limits  of  detection  of  turbulence. 

In  the  case  of  particle  displacements  due  to  a mean  velocity  U,  and  an  acoustic  dis- 
placement X4  the  spectrum  of  the  photodetector  current  centered  at  Kt(«  together 

with  one  satellite  is  shown  in  Fig.  10a.  Here  the  acoustic  signal  is  so  weak  that  the  acoustic 
satellite  is  everywhere  less  than  the  noise  accompanying  the  Doppler  shift  due  to  Kt^. 
Increasing  the  amplitude  of  the  acoustic  signal  raises  the  spectrum  of  the  satellite  to  posi- 
tion (Fig.  10b).  At  the  same  time  the  Gaussian  spectrum  at  jCll«  diminishes  to  , 

4^.  At  the  point  marked  by  the  circle  the  satellite  just  emerges  from  the  noise.  The 
amplitude  of  acoustic  signal  required  to  reach  this  condition  is  the  theoretical  minimum 
detectable  acoustic  signal  for  given  frequency  and  Doppler  broadening  . The  irre- 
ducible minimum  spectral  noise  is  that  of  the  photodetector. 

In  the  case  where  the  particle  displacements  are  due  to  turbulence  one  must  measure 
fluctuations  of  velocity  in  the  fluid  flow  by  observing  the  fluctuations  in  the  Doppler  fre- 
quency • If  the  frequency  broadening  of  these  velocity  fluctuations  are  of  the  same 
order  as  the  broadening  due  to  the  Doppler  ambiguity  (*=DA,  defined  above),  then  one 

t 

cannot  distinguish  between  fluctuations  It  and  fluctuations  giving  rise  to  the  Doppler 

/ _ 

ambiguity.  The  limit  of  resolution  of  U relative  to  XL  (i.e.  of  turbulent  change  velocity 
relative  to  laminar  flow),  is  then  given  by  the  condition 

Thu*  the  measurement  of  turbulent  14  is  limited  by  the  Doppler  ambiguity. 


i 
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Conclusion 

This  summary  has  served  to  note  the  following  points:  the  only  hydrodynamic 
velocity  that  is  measurable  in  a laser  Doppler  velocimeter  (LDV)  is  the  Eulenan  velocity 
which  is  an  average  overall  the  scattering  particles.  The  Lagrangian  velocity  of  in- 
dividual particles  (in  the  case  of  flow  of  many  particles)  is  not  measurable.  The  principal 
sources  of  spectral  noise  associated  with  the  scattering  volume  are  finite  transit  time  of 
particles  crossing  the  scattering  volume,  the  temporal  fluctuations  in  ft}  itself,  and  the 
random  motion  of  individual  particles  about  the  mean  Eulerian  velocity  as  functions 
of  particle  position  in  the  scattering  volume.  The  basic  model  of  the  LDV  is  a two  beam 
interferometer.  The  basic  model  of  the  electric  current  in  the  photodetector  is  that  of 
phasor  with  random  amplitude  and  random  phase.  The  calculation  of  minimum  detectable 
signal  is  based  on  comparison  of  the  size  of  the  acoustic  spectrum  to  the  spectral  noise  in 
the  absence  of  acoustic  signals.  The  irreducible  noise  in  the  system  is  that  of  the  photo- 
detection process.  The  possibility  of  reducing  all  other  noises  remains  open  to  research. 
Model  of  George  and  Lumley 

It  is  known  (see  Appendix  C)  that  the  light  intensity  S falling  on  the  detector  is 
rp  « a a 

S*  S-  J—  CL  . where  CL  is  a unit  vector  normal  to  the  photosensitive  area, 
and  is  the  reference  beam.  Assume  the  coordinates  of  the  photosensitive  area  are  x. 

it  J ''j  ll  A A X 

x and  that  an  element  of  area  so  that  j * , where  is  the  unit  vector 

in  the  direction  of  ^ . Let  M be  a constant  which  converts  light  intensity  into  electric 
current  (dimensions:  ).  The  total  electric  current  detected  by  the  LDV  is. 
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~JJ  {{ 1 5*  <*#*')£■  y * ")l  " 


(ii-i) 


in  which  there  are  an  scattering  particles.  If  the  scattering  particle  is  small  relative  to  a 
wavelength  of  laser  light  the  scatterer  is  a simple  monopole  so  that  in  the  coordinates  of 
the  reference  beam  (double  prime  coordinates)  the  scattered  field  is 


*X<  * 


(H-2) 
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t*  * 

in  which  )L.  is  the  position  of  the  particle  in  the  coordinates  of  the  reference  beam.  X* 

A*  p 

is  the  position  of  the  same  particle  in  the  coordinates  of  the  scattered  beam,  C^is  a scat- 
tering coefficient  constant,  and 


(II-3) 


(II.4) 


(Note  again  that  the  double  primed  coordinates  are  on  the  photodetector).  Since 

j]  ft dt)  tti Mtf!}  it’jis  £ (>.&:,«;) 

it  is  seen  that  the  electric  current  due  to  one  particle  is 

(II.5) 

The  reference  beam  and  scattered  beam  cross  each  other  at  angle  8 . The  photodetector 
is  set  to  detect  the  reference  beam  only. 

As  noted  earlier,  LDV’s  are  essentially  based  on  two-beam  interference.  The  “slit- 
sources"  are  two  laser  beam*  which  are  characterized  by  finite  v idths.  It  is  customary  to 
assume  that  the  distribution  of  light  amplitude  in  each  beam  width  is  Gaussian.  Thus  the 
light  beams  are  modeled  as 

A .-ifl.-t  .. 


> ( a")  ■ 

(*')=  {- 


(II -6) 


Kt  - 


(H.7) 
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in  which  are  amplitudes  at  the  center  of  the  beam  (as  measured  at  the  focus)  and 

<T  is  given  by 

O'- 

ltd 

(118) 

where  £is  the  focal  length  of  the  lens,  is  the  wavelength  of  light,  and  dis  the  dis- 
tance (in  the  cross-section  of  the  light  beam)  between  the  f£  intensity  points  of  the  beam. 

It  is  convenient  next  to  introduce  a local  coordinate  system^  fixed  at  the  intersection 
point  of  the  two  crossed  beams.  A particle  in  X is  located  at  Xjs  - f^J,  where 

one  defines  X^  in  terms  of  X , X , y t ^ by  the  relations, 


© * A I 

*'*  tv***  *&'“?). 

*"*  ®4  >)>  /'*  t <IL9> 


In  terms  of  the  net  current  due  to  a particle 


is 


b 


K - a Ke  /k*  * 


(11.10) 


(11.11) 

This  is  the  analog  of  the  classic  formula  for  a two-beam  interferometer  (see  Eq.  A5),  in 
which 


<r.  = 


jr_ 


T*a=-* 


• ^ = 4- 

1 A ' If  A 


(11.12) 
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The  exponential  term  in  Fig.  11-10  accounts  for  the  finite  cross  section 
of  the  laser  beams. 


Note  that  the  plane  /p,  fa  is  normal  to  the  plane  of  the  two  laser  beams,  and  effectively 
describes  the  observing  screen  ir.  the  conventional  2-beam  Young  interference  experiment. 
The  coordinate  is  normal  to  the  effective  observing  screen. 

Now  are  Eulerian  coordinates  of  particle  position,  and  so  are  functions 

of  Lagrangian  initial  position  a.  1 ($;i)  (see  Eringen,  Mech.  of  Continuies  p.  7). 

Letli/'tf  the  velocity  components  (in  the  local  coordinate  system)  of  identified  par- 
ticles,  i.e.  U-Ub,i)  etc.  Then 

0 

(11.13) 

However,  IX  (<}.-{)  is  not  observable  in  the  case  under  consideration  where  there  are  many 
particles  in  the  scattering  volume.  The  observable  velocity  is  l/ak)  which  is  an  ensemble 
average  over  all  scattering  particles.  To  perform  the  averaging  process  one  introduces  a 
probability  function  j/qjthat  a particle  of  specific  size  is  present  in  the  scattering  volume. 
The  scattering  volume  itself  is  specified  by  (4, inhere  TT(x)is  the  exponential 

in  Eq.  11.10.  Thus  the  observable  velocity  is 


lie  (*  )-' 

(11.14) 

and  is  the  expected  number  of  particles  in  the  volume.  At  any  instant  a particle  ve- 
locity deviates  from  the  observable  average  by  amount  V.  (&  >i)  - Uc  k)  . The  time-averaged 
deviation  is 

0 

(H.15) 

If  the  effective  displacement  over  all  particles  is  m then  by  definition, 
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Uo  (^)  ) Jt  ) 


(11.16) 


yn>.  ^ 


Hence  the  (Eulerian)  coordinate  of  the  particle  at  Q is, 


^ (a^  * **  4X^)  * A (a fit) 


(11.17) 


that  is,  a particle  which  is  initially  at  moves  an  average  distance  Xfy)(in  the  system) 
plus  a random  distance  (along  the  Xp  coordinate)  of  amount  Afapi).  Hence  the  total 
current  due  to  tm  particles  is 


n 


(11.18) 


This  current  is  again  a function  of  particle  positions  (cl')  and  is  unobservable.  The  current 
t (+)  averaged  over  all  £ will  be  observable.  Hence  one  again  introduces  the  probability 
function  £($}and  writes 


i(i).  (i,*)  <*3 


(11.19) 


Two  random  variables  are  in  tfaji^namely  random  particle  position (=4^  and  random  par- 
ticle displacement  (s  ,-t)).  Noting  that  GOS  (o/ 4^  = -/4av>ct  , one  can 


write 


ilt)*  XX  -(*•+<$  ) 


(11.20) 
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in  which 


F(*)--  (J  A ^ K [%  + 

<*(+>  -f  A K [a  -f  4 (*>*)[  fk\k 


(11.21) 


(11.22) 


(11.23) 

Since  the  photodetector  gives  a signal  proportional  to  the  time  derivative  of  the  phase  the 
output  of  the  velocimeter  is  a frequency  shift 


(n.24) 

The  observed  (averaged)  velocity  U»fi)  is  a random  function  of  time,  with  mean 
value  fluctuation  !(•(*) ^(^and  mean-square  fluctuation  = 

[v,(t)-  . The  local  (non-averaged)  velocity  u(v,-t)  is  also  a random  function,  with 

mean  value  U(jf  )and  fluctuation  Uf^t^-iffjrJ.The  covariance  I?,,  (*,*")  is  the  x-component 
time  averaged  product  of  at  £ and  £ \ 

T?„  lx,*')*  (u<M>  u( - 5 (£)] 

(11.25) 

Now  the  time  averaged  mean-square  fluctuation  U can  be  written  in  terms  of  pro- 
vided the  space  average  product  is  known.  It  can  be  shown  (George  & Lumley, 

p.  360)  that  if  the  particles  are  distributed  as  a Poisson  process,  then 
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uljL&JL  ■ 


}(*')*  Ue  + t{  *(a-a') 


(11.26) 


where  ii  the  expected  number  of  particles  in  the  scattering  volume.  The  formula  to/ 
1 is  then 

y»y-  — 

(;n).  <r,f,V, 


Noting  that 


•R„(x,xV  ^ S„ 


(11.27) 


and  defining 


9f(i.)  * •&($) 


where  wfJ^is  the  Fourier  transform  of  Vil),  it  is  seen  that 


(11.28) 


(11.29) 


(11.30) 


under  the  assumptions  that  the  fluctuations  of  !(.(*) conform  to  the  picture  "homogeneous- 
turbulence?  Note  that  It  * is  the  particle-averaged  variance  of  tt(x)(over  the  scattering 
volume)  and  It'  is  the  averaged  square  of  U(ji)  at  )^.  Under  the  further  assumption  that 
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Ne»l  and  '7tfk)|>  N£'  , one  reduces  this  last  formula  to 

U*  a ^ I 


(11.31) 


In  words:  the  time  and  particle  averaged  mean-square  fluctuation  of  the  observable  veloc- 
ity (in  theX  direction)  is  obtained  by  integrating  the  11-component  of  the  power  spectrum 
of  fluctuating  flow  at  wave  number  It  times  the  weighted  volume  of  ^ space  ) cik 
(which  tells  what  weight  to  assign  to  (k  The  weighted  volume  is  essentially  a 

statement  of  the  finite  boundaries  of  intersection  volume  in  real  space.  From  the  defini- 
tion of  it  « seen  that 


)--  Mb-  ( — + All 
‘ > 


in  which  A = (kt,  ^and 


(11.32) 


V - _L_  T - 1 7-1 

' JSC.  > 


/r  r> 


(11.33) 


This  form  of  ^H}^)  has  an  important  meaning.  Noting  that  (T, , IT,  ^ essentially  define 
the  * boundaries  of  the  scattering  volume  one  sees  that  and  ^ are  the  ^cut- 

off wave  numbersarising  from  the  spatial  Fourier  decomposition  of  the  scattering  volume. 
Similarly  fei,  are  the  wave  numbers  associated  with  the  spatial  Fourier  decomposi- 

tion of  the  fluid  velocity  field.  Now  if  *,>e,  , that  is,  if  the  spatial  variation  of  fluid 
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V..OC*,  » mam  .h»  *.  ~ vo.u»,  in  .h.  x-direction  the- 

small.  This  means  that  the  n,  component  (or  its  corresponding  velocity  change)  is  se- 


verely attenuated  by  the  measurement  volume.  This  conclusion  is  of  prime  importance  in 


the  question  of  velocity  resolution  of  the  velocimeter. 

The  subscript&of3||  have  special  significance.  Let  ^(jjand  (fc)  be  two  velocity 
vector  fields  (in  Cartesian  index  notation).  Then  the  velocity  correlation  tensor  is  the 
spatial  average 


Kj  V4  It)  (*-*£) 


(11.34) 


In  particular  1?,  | (a  ) is  the  correlation  of  a velocity  vector  field  at  Jt  and  at  The  spec- 
tral density  (tensor)  is  the  Fourier  transform  of?^f*JP  and  (kj  is  the  spectral 
density  component  of  the  average  product  atjt  and  XfA  of  the  same  velocity  field.  Since 
only  one-dimensional  velocity  spectrums  can  be  measured  one  defines 


F,/  ( k, ) - I„  (k,,  k.MUUK 


(11.35) 


The  notation  Fjj(fei)  means  the  following:  the  11  subscript  signifies  the  same  velocity 
field,  and  superscript  1 indicates  tb  j first  (sag  x)  component.  Note  that 


M'*  = if 


(11.36) 


i.e.  the  mean-square  fluctuation  of  the  x-component  of  velocity  spatially  averaged  is  given 
by  the  integral  shown.  The  measured  one-dimensional  spectrum  ^(kt)  is  quite  different 
since  gives  the  temporally  averaged  mean-square  value  of  the  average  (effective)  velocity  in 
the  scattering  volume, 
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(11.37) 


«.'(*>*  * J Po  (It,)  A 

-O© 


and 


-a* 

,w^-  h'k?\dUi> 


, (H-38) 

Note  again  that  is  the  measured  spectrum  and  I", , is  the  theoretical  spectrum. 

The  ratio  /F]|(^i)  **  the  velocimeter  transfer  function.  Typical  values  based  on 

the  assumption  of  isotropic  turbulence,  and  a special  form  of  the  energy  spectrum 
is  shown  in  Fig.  4,  p.  322,  George  and  Lumley . 

Particle  Displacement 

The  observed  x-component  instantaneous  velocity,  averaged  over  all  particles,  is  Volt), 
This  is  a random  variable  since  the  number  of  particles  varies  from  moment  to  moment  in 
the  scattering  volume.  If  we  assign  a probability  distribution  to  tds  we  can  find  its 

•WMB  — ■ A 

mean  value  (see  Fig.  ||)a  The  effective  displacement  due  to  Z i$(i)  is  )(H),  where 


A ft  — - — , 

x U)=  \ tuft,)#,  -at 


(11.39) 


Here  U.  is  a fictitious  average  constant  velocity  such  that  the  effective  displacement^  is  the 
product  U t . Now  X^)  itself  is  a random  function  of  time.  Thus  at  timet  the  fluctua- 
tion  in  observable  (x-component  ofjdisplaccment  is 


x'M-XM-  X«) 


(11.40) 
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Fig.na>'*b> ,,c 

In  distinction  to  observable  UofO  there  is  the  unobservable  x-component  of  velocity  of  a 
particular  particle  , vi $ . The  difference  between  Wit)  and  averaged  over- 
time, results  in  a random  x -displacement  The  (random)  instantaneous  dis- 

placement is  therefore 


In  this  form  A (v >s  unobservable  because  of  A method  of  averaging  over  all  par- 

ticle position  is  discussed  in  the  next  section. 

Measurable  Correlation  of  the  Photodetector  Current 


We  desire  to  form  (1)  the  total  current  L(t)and  (2)  the  normalized  correlation  of  current, 

0)  id)*  f i (td) $(i)A 

(Z ) id)  ilt')  <■■-«> 

<Jf<) 


(11,43) 

in  which  the  random  variables  are  Mify  The  function  Q(a ) is  already 

assigned  (see  Eq.  11.26).  To  X^)we  assign  the  probability  density 


*1 


(11.44) 

and  toAfo^l  considered  a variable,  we  assign 


?[a  («,♦),  Aft, 01* 


(George  and  Lumley,  Eqs.  4.1.2,  4.1.3)  (11.45) 

in  which  % and  ^ are  standard  deviations.  A full  calculation  leads  to  the  (measureable) 
result  that 
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(George  & Lumley.  Eq.  4.1.6) 

This  is  the  correlation  coefficient  for  the  Doppler  current  between  times  t and  t . 

We  return  now  to  the  quantity  fa/t)]  which  '.is  the  unobservable  particle  velocity 

at  point  X , (where X is  the  Eulerian,  or  fixed,  coordinate),  and  the  averaged  (over  particles) 
particle  velocity  lit  (+)  (independent  of  position  in  the  volume)  and  form  the  mean-square 
fluctuation. 


uUMj-v-hI 


(11.47) 

This  quantity  differs  for  each  X . If  we  average  over  volume  the  result  is  an  average  mean- 
square  velocity  deviation,  labelled  (Alt)  . It  is  the  mean-square  difference  between  the 
averaged  mean  square  fluctuation  at  the  center  of  the  scattering  volume,  and  the  mean- 
square  velocity  fluctuation  averaged  over  the  entire  volume.  Thus  it  is  a measure  of 
velocity  deviations  across  the  scattering  volume.  If  the  fluid  medium  is  incompressible, 
and  statistically  homogeneous,  stationary  and  isotropic,  it  is  directly  seen  that 

or  (11.48) 
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(11.49) 


Hence  accounts  for  the  difference  between  the  volume-averaged  fo(^i)  spectrum 

and  the  true  spectrum  at  the  coordinate  point.  We  next  assume  the  “optical  size”  of  the 
scattering  volume  in  the  x-direction  (namely  Kr,)  is  of  such  magnitude  that  K <T^>^  | 
Taking  Tig  being  the  optic  wavelength,  this  assumption  means  that  the 

scattering  volume  is  many  wavelengths  of  light  in  the  x-direction.  Now  the  autocorrelation 
of  photodetector  current  depends  on  the  quantities  , where  is 

the  standard  deviation  of  the  assumed  Gaussian  form  of  the  joint  probability  of  the  efffec- 
tive  random  displacement  of  the  averaged  particle  at  the  center  of  the  scattering  volume 
at  times  £;t^and  0^  is  the  standard  deviation  of  the  assumed  Gaussian  form  of  the  joint- 
probability  of  the  random  displacement  difference  A (*; between  a particle  at  £and  the 
“averaged”  central  particle.  For  small  time  difference (“t-ty! 


r„*  a tt'"  (<"<’/;  ofatiZ?  (t-t'f 


(11.50) 

(again:  U a is  the  mean  square  fluctuations  of  the  center;  GmJ*  is  the  mean  square  fluc- 
tuations relative  to  the  cente^.  Using  ly^  j , kV;  ■I.  as  a characteristic  number,  it  can 
be  shown  that  when  K C*,  » • , F(t)  F (<■')  3?  5 (t)Q(if)>  F(l)  OrtJ 

CosK(X-X  r Sxp ^ -X')  , valid  to  the  order  of  the  characteristic  number.  Thus, 

in  this  approximation, 

KttV) a 1 Kr/(t,t')j. 
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(11.51) 


A more  convenient  form  is 


Ut ) i fe+t) 
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(11.52) 

in  which  the  superscript  ^ means  time  average,  and 
■i 

F(t)  ffm) . A f_ *V(r),  Xfo-Xftn)? 
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(11.53) 

The  availability  of  the  normalized  correlation  of  the  photodetector  current  leads  directly 
to  a calculation  of  the  power  spectrum . It  will  be  convenient  to  discuss  each  factor  of  the 

correlation  separately.  + » 7 

The  Fourier  transform  of  Q(r)»  (Si  K [/(t)- X 6-^(«“ely  of  j-  l 

Q(»»  { j-  7^  J + { 


The  Fourier  transform 
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(11.54) 
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Here  (as  before)  it  is  desirable  to  define  all  relevant  terms  again;  thus.'U*  is  the  time 


averaged  value  of  ^is  the  time  averaged  value  of  the  time  averaged 

value  of  the  non-random  X (t )/f  > is  the  time  averaged  mean-square  deviation  of  the 

random  velocity  off  center  from  the  velocity  on  center.  The  Fourier  spectrum  <I(k)of  the 


Dhotodetector  current  is  then 


ilM-  pT(»)  Q (<■>- *>,)</«, 
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(11.55) 


In  words:  the  normalized  power  spectrum  of  the  photodetector  current  arising  from  the 
motion  of  particles  which  scatter  light  in  a heterodyne  laser-doppler-velocimeter  consists  of 
two  Gaussian  peaks  centered  at  £ Isl^  where  *14 » is  the  time  average  of  the  velocity  u»k), 
which  itself  is  the  (observable)  average  of  the  velocities  of  all  the  individual  particles  at 
some  instant  t present  in  the  scattering  volume.  The  symbol  K is  9 *4  , where 

eA^  X 

is  the  included  angle  between  the  two  beams 
of  the  interferometer.  The  Gaussian  peaks  are  broadened  by  the  effects  of  several  time- 
varying  quantities.  The  first  of  these  is  ^ /at,1  . To  see  its  full  meaning  we  discuss  it  as 
if  it  were  the  only  quantity  affecting  the  spectrum.  The  entity  0^  is  recognized  as  the 
effective  size  of  the  scattering  volume  (in  the  x-direction),  and  jX  is  the  time  averaged 
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particle  velocity  of  all  the  particles  in  ‘.he  absence  of  turbulence.  Essentially  U is  the  veloc- 
ity of  all  particles  through  the  volume,  acting  in  laminar  flow.  Hence  J is  the  effective 
time  of  transit,  and  is  the  bandwidth  of  frequencies  associated  with  this  time  of 
transit.  In  effect  the  laser  light  is  modulated  by  the  on-off  switch  of  the  finite  scattering 
volume,  resulting  in  an  uncertainty  in  the  Doppler  shift  that  is  ascribed  to  l^. 

When  the  flow  becomes  turbulent  one  sees  two  additional  effects  which  serve  to 
broaden  the  spectrum.  First  there  is  the  effect  of  the  fluctuation  of  the  (particle-averaged) 
central  velocity  Ui«)  about  its  mean  . The  time  scale  cf  this  occurrence  is 


t--  (K  V)  * 


Uc(*)-  Vo 


and 


_ <IL56a> 

the  associated  bandwidth  is  (K  Uq  */  . The  fact  that  Voftjexists  introduces  an 
uncertainty  in  the  Doppler  ascribed  toVt . Second,  the  velocities  of  the  particle  are  no 
longer  laminar.  At  any  time£  there  is  a difference  between  the  velocity  of  an  individual 
particle*^  (j  ,4)  and  the  velocity  averaged  over  all  particles  The  time  average  of 

the  square  of  this  (random)  difference  is  feu)1  l.  The  time  scale  of  this  turbulence-induced 
velocity  variation  across  scattering  volume  is 


-i 


1 — A (IU6b) 

and  the  associated  bandwidth  is  ( AU)  J . 

Power  Spectrum  of  the  Photodetector  Current 

In  the  presence  of  a sinusoidal  phase  modulation  due  to  an  acoustic  signal  we  repre- 
sent the  photodetector  current  by  the  form 

(11.57) 
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Since  this  function  is  periodic  one  can  separate  the  modulation  index  from  the  time 
variation  , viz 

_ 30 


in  which 


F«)  F (<’) 


f’m 


(11.60) 

The  power  spectrum  of  l ft) can  be  obtained  by  taking  the  Fourier  transform  of  this 
function.  We  note  first  that  the  FT  given  by  Eq.  11.55  accounts  for  all  terms  in  the  power 
spectrum  except  the  acoustic  term.  By  the  modulation  theorem  (p.  108  Bracewell)  if 
F the  transform  of  -fftjthen  the  transform  of  4 ft) Cot  43^  tf  ^ ^ ftiMtfo)  4 i F 
Thus  the  power  spectrum  T(w)of  the  photodetector  current  in  the  presence  of  an  acoustic 
signal  is  given  by 
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(11.61) 


,3  (A*)5 

K’iiJ5  -t  K + kW> 

v 1 a r,' 

In  words:  In  the  presence  of  an  acoustic  signal  which  gives  an  average  displacement  of 
to  all  the  particles  of  the  scattering  volume,  and  in  the  presence  of  a steady 
(constant)  average  velocity  Vo  of  all  the  particles,  the  power  spectrum  of  the  photodetector 
current  exhibits  the  following  character: 

(1)  there  are  two  Gaussian  peaks  at  frequencies  £ (=Cose  o$  w>aoV 

(2)  the  standard  deviation  of  the  distributions  in  each  case  is  the 

Doppler  ambiguity  A-B. 
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(3)  there  are  four  additional  satellite  Gaussian  peaks;  two  centered  at  KU0  i 

and  two  centered  at  — A (~  Cause 

(4)  the  relative  amplitude  of  the  satellites /tfz}  to  the  main  peak  is  proportional 
to  the  ratio  of  the  Bessel  functions  To(K*s)/t/K>s)  . 

Instantaneous  Signal 

The  photodetector  current  i(t)  has  been  modeled  as 

i (*).  Ff+)  an  [KXW]  + <5 

in  which  Fft)  ^j(t)are  functions  obtained  by  integrating  over  random  initial  position  and 
random  deviation  from  initial  position  for  all  particles  in  the  illuminated  volume  and  X&) 
is  the  time  averaged  displacement  of  the  random  “center  velocity”  Voty-  This  formula  can 
be  altered  to  give  a different  (physical)  insight  into  the  detection  process.  The  basic  pro- 
cess (as  noted  earlier)  is  a two-beam  interferometer  experiment;  which  calls  for  an  inten- 
sity of  light  (note:  l(t)  is  proportional  to  light  intensity)  with  a sinusoidal  component 

(see  Eq.  11.20).  In  the  above  modelXf*)  is  the  analogy  of  and  the  illuminating 
(i.e.  scattering  volume)  source  has  a random  character  described  by  Rt)  and  Cf(t).  One 
can  say  that  the  equivalent  two-beam  interferometer  exhibits  a fluctuating  amplitude  versus 
time.  Thus  one  can  write  tft)  as  a phasor  with  fluctuating  amplitude,  and  fluctuating 
phase, 

I 

i (t).  Ah) 


(11.63) 


At  time  t there  are  four  random  quantities,  ^*(^t  &»),<?  (*l)- 

Similarly  at  time  ^ there  are  four  additional  random  quantities.  The  autocorrelation  1?Ni 
i (i)  is  that  of  an  eight-dimensional  Gaussian  distribution.  It  is  known  that 
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(11.64) 
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Now  according  to  Part  II.  Sect.  1 (above)  the  functions  F,  G have  Gaussian  correlations. 
Thus 


From  this  one  obtains 


hr)=  O'?') 


The  Fourier  power  spectrum 


is  therefore  given  by 


N(a)*  F^4>(t)}  * ^(r)e:avjr 
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(11.65) 


(U.66) 


(11.67) 


(H.68) 


This  is  approximated  by  the  two  values  depending  on  the  relative  magnitude  of  Si  and  DA^ 


to 
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N (°)  = 

a,  lloi  DA 

(Si  < da). 
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(11.69) 


We  define  the  Doppler  ambiguity  spectrum  as  the  power  spectrum  due  to  temporal  char- 
acter of  the  center  x-  component  of  velocity  Vo(t)  (<•€.  t/.  + lO  This  is  given  by  spec- 
trum of  finite  transit  and  spectrum  of  Uq.  We  define  the  turbulence  power  spectrum  as 
that  associated  with  the  fluctuations  in  phase  as  defined  above.  Now  the  rate  of  change 
the  phase  is  the  important  contributor  to  the  Doppler.  One  labels  this 

^ tio (t)  -<f> 


(11.70) 


Both  Uo  and  C can  be  written  as  a Fourier  expansion 
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(11-71) 
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Since  the  amplitudes  <J$i)and<M(Xl)  are  orthogonal  (in  frequency),  we  write 
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(11.73) 
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Now  the  limit  of  resolution  of  turbulence  spectrum  is  set  at  a frequency!^  such  that  the 
Doppler  ambiguity  spectrum  is  equal  to  the  turbulence  spectrum,  i.e. 

(Doppler  Ambiguity)«(turbulence) 


H((L) 


(11.74) 


(dimensions:  (l)Kl  M";  («)R  i HV*i(»)U':_Ms' ; (4)0:  S^fS)  tf  : S-<  ). 
It  is  convenient  to  define  tie  /n  » ^ej  Ajc  |=  fee  j|<f  ' tj;  / which  f is 


the  Kolmogorov  rate  ofdissipation  of  turbulent  energy  per  unit  mass,  7 is  the  kinematic 
' is  the  Kolmogorov  microscale.  Adopting  the  approximation 


viscosity,  and  f|a  j 
that 


for^  - [K'.V  + K j 


we  are  required  to  find  appropriate  values 
For  smae.1  scattering  angle  one  can  write 
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(11.75) 

(George  and  Lumley,  Eq.  3.2.8). 

Hence 
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(11.76) 

(George  and  Lumley,  Eq.  5.5.2).  Here  ^ is  the  Reynolds  number  based  on  the  smallest 
length  that  can  be  resolved  in  the  mean  flow  direction,  i.e. 
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(11.77) 

The  above  equation  for  p (ko\  **  a 'tey  reiult.  In  words:  the  one  dimensional  velocity 
spectrum  of  turbulence  {-‘{(fc  \ can  be  determined  for  given  Reynolds  number  Tit  of 
flow,  given  angle  6 between  the  two  beams  (scattered  and  reference)  of  the  equivalent 
2-beam  interferometer,  and  given  size  Ofot  scattering  volume  in  the  x direction  as  it  appears 
in  the  parameter  K| 1 • In  the  square-root  brackets  the  first  term  (namely  unity) 

represents  the  (relative)  contribution  of  finite  transit  time.  It  contributes  the  Doppler 
ambiguity  spectrum.  The  second  term  gives  the  (relative)  contribution  of  turbulence  to 
the  broadening  of  the  spectrum  of  velocity.  For  fixed  ^ and  6 the  spectral  broadening 
in  F„  ( f ) is  a function  of  the  size  parameter /Cf)  . Two  regimes  of  variation  of 
Fii  (L)  wth  Ki  rj  can  be  found.  In  the  first,  is  made  very  small.  This  is  equiva- 
lent to  making  the  scattering  volume  very  large  (relatively  speaking).  The  spectrum  broad- 
ening then  is  due  to  turbulence,  and  increases  as  the  size  of  the  volume  of  scattering  in- 
creases. In  the  second  regime  is  made  very  large,  i.e.  the  scattering  volume  is  made 
very  small.  The  spectrum  broadening  is  then  due  to  finite  transit  time  (Doppler  ambiguity 
broadening),  its  magnitude  increasing  with  reduction  in  size  of  the  scattering  volume.  Be- 
tween the  two  regimes  there  is  evidently  a minimum  broadening  of  the  spectrum.  This 

I ^ 

occurs  at  the  value  of  j which  makes  a minimum.  By  inspection  the  con- 

dition for  minimum  is  given  by  making  the  contributions  of  finite  transit  time  and 
turbulence  of  equal  magnitude,  i.e. 
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The  smallest  value  of  the  turbulence  spectrum  that  than  be  measured  with  minimum  ambi- 
guity due  to  finite  transit  time  is  thus 
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(11.79) 

In  words:  when  it  is  desired  to  use  the  laser  doppler  velocimeter  to  measure  turbulence 
it  is  found  that  the  power  spectrum  of  velocity  turbulence  is  broadened  by  “large”  scat- 
tering volume,  which  renders  the  determination  of  the  "spectral  line”  associated  with  a 
given  wavenumber  k,  uncertain.  The  broadening  of  the  spectrum  is  reduced  by  decreasing 
the  size  of  the  scattering  volume.  Hence  one  expects  that  the  smallest  available  volume  of 
scattering  would  generate  the  narrowest  spectral  lines.  However  as  the  volume  of  scattering 
is  reduced  the  broadening  due  to  finite  transit  time  is  increased,  making  the  spectral  line 
of  velocity  again  uncertain.  Evidently  the  narrowest  spectral  line  is  achieved  by  using  a 
compromise  size  of  scattering  volume.  This  compromise  size  is  fixed  by  » The 

smallest  spectral  level  F,|  (kj  that  is  optimized  by  scattering  volume  is  [FJ|  opt. 

Hence  the  highest  value  of  that  can  be  measured  with  minimized  ambiguity  due  to 

I 

finite  transit  time  is  found  by  expressing  Fj(  (k%)  as  a function  of  some  turbulence  model, 
and  solving  for  k,  implicitly.  ^ 

The  implicit  solution  for  k,  can  be  illustrated  by  use  of  the  model  of  isotropic  tur- 
bulence. Here, 
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The  form  of  Eft).  particularized  by  choice.  Let  Ef*)be  that  of  Pao  (Phys.  Fluids  8, 


1063,  (1966)), 
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Substituting  this  form  into  the  above  integral,  and  performing  the  integration  leads  to  a 
model  of  F»fe)  which  can  be  plotted  versus  . On  the  same  graph  one  can  plot 
Tp  VClM  for  various  choices  of"R*  and  fr.  The  intersection  of  the  turbulence  model 
(Pao’s  model)  with  the  ambiguity  model  J opt.  gives  the  value  of  ^ , which  is 

the  wavenumber  for  which  the  ratio  of  the  turbulence  spectrum  to  ambiguity  spectrum  is 
unity.  This  is  the  largest  value  of  wavenumber  of  the  turbulence  power  spectrum  that  can 
be  measured  with  the  minimum  ambiguity  (that  is,  under  conditions  of  the  minimum 
spectral  broadening  due  to  finite  transit  time),  for  fixedly  and  $.  ^ 

Seeing  that  every  Choice  of  Re  and  & is  equivalent  to  a choice  of  k , one  can  plot 

) versus  , Suppose  we  desire  the  maximum  l (•«.*!  4)  to  be  unity. 

Then  the  graph  shows  "R*  Q.lO  . For  an  experiment  in  which  6 is  fixed  such  that 
• ft  ' 

-0.2  and  the  wavelength  "fa  C ■3*/0  Cm  it  is  found  from  the  Reynolds  number 
that  the  smallest  measurable  velocity  (the  optimal  velocity  in  the  above  sense)  is  K.  2f  /-S’ 
cm/sec. 

In  actual  experiments  the  optimum  scattering  volume  can  be  determined  by  a simple 
formula  once  the  Reynolds  number  and  angle  of  scattering  are  known, (see  Appendix  B). 
Concluding  Remarks  to  Part  II 

In  Part  II  the  laser  Doppler  velocimeter  is  mathematically  analyzed  using  the  model 
of  George  and  Lumley.  The  photodetector  (electrical)  current  L (i)  in  this  model  is  a 
phasor  which  is  a random  function  (in  both  amplitude  and  phase)  of  the  particle  displace- 
ment across  the  laser  fringes  (see  Eq.  11.20).  By  assigning  probability  densities  to  the 
fluctuating  component  of  the  volume  averaged  displacement,  and  to  the  deviation  displace- 
ment of  an  individual  particle  from  the  volume  averaged  displacement,  one  can  find  an 
ensemble  average  of  i(t),  and  then  form  the  normalized  autocorrelation  of  i(t)  at  times  t 
and  t (see  Eq.  11.46).  Under  the  assumption  that  the  size  of  scattering  volume  is  several 
wavelengths  of  laser  light  large,  and  that  t't  it  not  too  great,  one  can  calculate  the 


63 


power  spectrum  of  the  random  current  t(t)(see  Eq.  11.55).  The  particle  displacement  is 
then  considered  to  contain  a sinusoidal  component  due  to  the  presence  of  an  acoustic 
signal  in  the  scattering  volume.  The  calculation  of  the  power  spectrum  of  photodetector 
current  is  then  repeated  leading  to  the  result  shown  in  Eq.  11.63.  In  applications  to  mea- 
surement of  turbulent  velocities  one  is  required  to  find  the  true  one-dimension  power 
spectrum  (k)  of  these  random  velocities  for  different  scale  sizes  (i.e.  wavenumbers  k). 
Because  of  finite  size  of  the  laser  beam,  and  the  existence  of  random  motion  from  parti- 
cle to  particle  within  the  scattering  volume  the  measured  spectrum  fb  (^differs  from  the 
true  spectrum  F()  (k) , The  measured  spectrum  Fo  (k)can  be  summed  for  all  values  of  k. 
to  give  the  time-averaged  mean-square  fluctuation  of  the  observable  Eulerian  velocity 
The  true  spectrum  F|(  (^)  can  be  summed  for  all  values  of  k to  give  the  time- 
averaged  mean-square  fluctuation  of  the  unobservable  velocity  tr^which  is  < u*H )* 
whereAU  is  the  deviation  of  across  the  scattering  volume.  Under  certain 
assumptions  fit(k)  can  be  calculated  (see  Pao’s  spectrum  Eqs.  11.80,  11.81).  A theoretical 
limit  of  spectral  resolution  (of  the  turbulent  velocity  spectrum)  can  be  calculated  using  a 
good  model  for  the  noise  due  to  the  phase  fluctuations  of  the  equivalent  phasor  i(i) i The 
model  is  that  of  filtered  noise  in  an  FM  receiver,  called  By  approximation  H(&)  « 

taken  as  a white  spectrum  with  a fixed  value  MM -<5.34?  *(*ee 

Eq.  II.69).  The  limit  of  resolution  of  the  turbulent  velocity  spectrum  F)(  (k)Can  then  be 
obtained  by  setting  HO)  as  the  threshold  (see  Eq.II.74) of  (k^The  limit  thus  found  is 
Eq.  (11.76)  where  F„  ft.)  is  found  to  be  a function  of  the  Reynolds  number  (based  on  the 
laminar  velocity  "u)  and  the  wavenumber  of  the  laser  light  IC-(^  sin^  , where  & is 

the  angle  between  the  two  beams  of  the  laser  interferometer.  When  the  magnitude  of 
F||  [kj  is  found  for  given  Reynolds  number  and  O , the  value  of  ko  is  then  (implicitly) 
determinable.  This  is  the  "cut-off  wavenumber”  (see  Eq.  11.79).  At  this  wavenumber  the 
ambiguity  due  to  finite  transit  time  and  turbulence  is  minimal.  It  is  essentially  the  largest 
wavenumber  (hence  the  smallest  velocity)  that  can  be  measured. 

The  model  of  George  and  Lumley  reviewed  and  extended  above  is  the  most  complete. 
Additional  models  are  briefly  discussed  in  Appendices  F,  G,  H and  I.  These  serve  to  dis- 
play somewhat  different  points  of  view. 
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1 


Part  III 


Communications  Theory  Model 

In  the  model  of  George  and  Lumley  (Part  II)  the  photodetector  current  L (i)i&  a phasor 
with  random  amplitude  and  random  phase, 

i(t)=  [f(t)  < [ KX(t) - *(t)\ 

we  restrict  the  generality  of  this  model  and  consider  only  the  cosine  function,  replacing 
the  random  amplitude  by  a constant.  Furthermore  we  take  the  displacement  )((t)u>  be 
deterministic,  and  write 

X(t)* 

in  which  tij  )lf  and  ^ are  defined  in  Part  II.  The  only  random  quantity  is  a phase  4>(i) 
This  fatt)  represents  noise  in  the  laser  Doppler  velocimeter  and  is  modeled  as  if  shaped  by 
an1?C  filter,  bandwidth  power  spectrum  %l(f) and  covariance^  given  by 


J 


(HI.1) 


It  is  convenient  to  regard  the  frequency  K u a carrier  frequency,  the  quantity 

as  a phase  modulation  of  this  carrier,  and  a noise  modulation.  In 

complex  notation  the  photodetector  current  then  is  assumed  to  have  the  form 

i (t  ) = A 1 | 

Since  this  is  a periodic  function  it  can  be  represented  by  a Fourier  expansion, 

Imteisti  ld)ot-L4(t) 


(IH.2) 


The  covariance  function  of  this  modulated  carrier  is  then  derived  to  be 

K,  (C)  * / I ) 

9 » i TTT 


(III.3) 


(Middleton,  “Intro.  Statis.  Comm.  'Theory.”  p.  612). 
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Here  £ - the  Bessel  function  of  first  kind  of  integer  order,  and  £=  I when  /th*  C 

s i tjl&i  m % 0 . The  spectral  distribution  corresponding  to  this  covariance 
function  is  where 

0* 

Z 


Z * 

/n'o  /f>!  ^ ys*» 


A _ ^sl  - &)fl  ♦ /)M 

-**'  — Z5?r“ 


(III.4) 


(Middleton,  p.  613).  (Note,  only  positive  frequencies  appear  here.) 

This  power  spectruin  has  a dual  character.  It  is  noted  in  the  formula  for 

(kWyi'i  fc*)*D-l<N+|£f...  ] 


'that 


(III.5) 


The  Fourier  transform  (ft)  of  unity  is  a delta  function  at  the  origin  of  frequency.  The 
FT  of  cos  <dp^  is, 

FT  { cei  av*}=  { &(&-&•)  4 £ Kmu3*) 

If  we  multiply  this  time  function  by  cos  m al$tthen 

FT  ( oh  a.t  «**t  } = ^ S(<D  -".♦*>*)»)  X 

Thus  centered  around  (O*  tide  there  are  two  delta  functions  (or  discrete  spectra)  for  every 
value  of  m.  The  magnitudes  of  these  spectral  lines  are  proportional  to  T/m  ( ,^Xj  ) We 
return  now  and  examine  HM"  ht  the  special  value  n * 0 and  (V  «•  ( t)o  V“*fV>tiat  In  this  case 
ft,  so  and  ft  s - . The  value  of  the  first  term  is  known,  vi?, 

I4*  r-» 


po  ^4/C  1 J 


(III.6) 


(see  Morse  & Feshbach,  Meth.  Theo.  Phys.  813).  Thus  the  term  n^^ayields  a discrete 
spectrum  and  the  term  jt  ^ 0 yields  a continuous  spectrum,  L e. 
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KI-).  /VK.WWl4#i^)] 


(III.7) 


Now  let /ins)  , and  take  fa  . ^ * to-  60,  -i  (Ot  .We  then  choose  a 

frequency  (V-tty  -MU.  Then^+)-»  0 and  Hence  for  the  choice  m - I 

the  spectrum  again  has  both  discrete  and  continuous  values, 

, ^ 

. T,  l(KxsM  l * S ( m-teo+tom 

1 A/«i..  V 


Wl 


/Vi 


/ tfi-flOo+tti  \v  4 / \x 

\ AUk  I ' Mu  ) 


)} 


(III.8) 


A sketch  of  the  total  spectrum  is  shown  below: 


FIG.  \Z 


We  next  consider  the  continuous  spectrum  of  M*o  at  the  value  of  ti~  tic*  and  take 
Kx  to  be  small.  Then  the  noise  due  to  the  peak  at  4^  is 


> _ 


K* 


*>  A e To’(Kxs)i 

A tin 


/ms  o 


kW  , 

14  /^L  f 


(III.9) 
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The  continuous  spectrum  of  a i*s  I a+  k>-  <5tfot  tig  («, 


.a  -fe)'  , 

( O).  +% )'  £;  e T,  ( Kxs)  K V } /+ 


(III.10) 


i '♦ean 

There  are  two  signals  in  this  special  case:  (1)  noise  due  to  the  main  peak  at  the  Doppler 
frequency  KtTo  (2)  noise  at  the  frequency  of  the  acoustic  satellite.  We  now  assume  that 
the  criterion  for  detection  of  an  acoustic  signal  is  that  it  should  equal  the  satellite  noise 
at  the  satellite  frequency.  Hence  we  state  that  an  acoustic  signal  of  magnitude  equal  to 
the  continuous  noise  ai  the  satellite  frequency  will  be  detectable  if  it  rises  above  the 
noise  contributed  by  the  Doppler  ambiguity  of  the  central  peak  (t>e,  in  the  absence 

of  an  acoustic  signal.  This  means  that  we  must  find  a value  of  such  that, 


or 

I,1  (iW  „ 

/ fel*  \ 

W * (***,) 

(111.11) 

u \A 

:fl 

These  formulas  will  be  used  in  the  following  sections  to  make  numerical  estimates 
of  the  minimum  detectable  modulation  index. 
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Part  IV 


Calculations 


Summary  of  Measured  Displacements 

Massey  (1968,  “An  Optical  Heterodyne  Ultrasonic  Image  Converter,"  Proc.  IEEE  56. 

2157)  used  a local  oscillator  heterodyne  apparatus  to  sense  the  vibration  amplitude  distri- 
bution on  a reflective  resonant  diaphragm  placed  in  a liquid  acoustic  medium.  The  phase 
modulation  of  the  laser  beam  was  ^(t)  * 4 IT  $ in  which  “Xj  is  the  la- 
ser wavelength,  is  the  acoustic  signal  frequency  and  “ the  amplitude  of  vibra- 

tion of  the  diaphragm.  The  dominant  noise  process  was  shot  noise  in  the  photodetector. 
For  a circuit  bandwidth  B,  and  electronic  charge  e the  rms  noise  current  is 


in  which  currents  in  the  detector  due  to  the  local  oscillator  beam,  or  signal 

beam  acting  alone.  Assuming  a total  laser  power!#  quantum  efficiency  ft . and  setting 

> • 

. one  finds  that 

jlce 

in  which  h is  Planck’s  constant,  is  the  speed  of  light.  Now  the  minimum  detectable 
amplitude  is  found  by  equating  the  noise  ^ to  the  first  side  band  current.  It  is  found 
that 

k 


7 - I 


Thus  the  minimum  detectable  amplitude  depends  on  the  square-root  of  the  circuit  band- 
width and  inversely  as  the  square-root  of  the  laser  power.  In  a typical  case 
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.•cc 


i 


i 


! 


0 


i 


I 


J 


XE  « 4,-  I MHz 

"g*  2o  kHz. 

^ • Sx/fiT2, 

7»*  io'V 

it  is  found  that 

Z„iti  * 2.2  C-^eak) 


Yeh  and  Cummins  (1964,  Appl.  Phys.  Letters.  4,  176)  also  used  a LDV  to  measure 

fluid  flow  based  on  the  local  oscillator  heterodyne  principle.  They  considered  that  they 

_c,  .1  o 

could  detect  constant  velocites  as  low  4<|D  nS  at  a scattering  angle  of  30  . 

Edwards  et  al,  (1971,  J.  Appl.  Phys.  42  837)  analyzed  and  measured  steady  flow  in 
a LDV  based  on  the  local  oscillator  heterodyne  technique.  In  order  to  resolve  the  Dop- 
pler signal  they  estimate  that  the  (constant)  velocity^ must  be  such  that 

v > D K , ^ K 

in  which  D is  the  translational  diffusion  coefficient,  andlKW  4IL  /usn  ^ . In  the 

case  of  Tig  = C'JlS  >/0  ^8  = 4$°,  "Jr  1 0 AlV  . they  found  that  approximately 

|v|  > 10'*  Ms'1 

Wang  (1972;  J.  Phys.  E.  763)  found  the  following  S/N  dependencies; 


local  oscillator  heterodyne:  S/N 


'V 


differential  or  symmetrical  heterodyne:  S/N 


in  which ff)  is  the  partical  number  density,  ^ is  the  Mie  scattering  cross  section  of  the 
partical,  Xg  is  the  laser  wavelength,  Is  is  the  linear  dimension  (say  depth)  of  the  scatter- 
ing volume,  and  j is  the  average  partical  size.  Thus  for  the  case  of  the  local  oscillator 
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heterodyne  a high  particle  concentration  and  small  focal  volume  are  needed.  It  is  insensi- 
tive to  particle  size  and  has  broad  frequency  spreading.  In  contrast  the  differential  (or 
symmetrical)  heterodyne  system  requires  low  particle  concentration,  is  sensitive  to  particle 
size  and  has  low  frequency  spreading. 

Criterion  for  Signal  Detection 

We  assume  to  be  small,  and  consider  only  two  terms  (i.e.  4W  / ) in  the 

equation  for  the  power  spectrum  of  the  photodetector  (Eq.  11.61).  We  then  state  the  fol- 
lowing criterion  for  signal  detection:  an  acoustic  signal  is  considered  detectable  if  it  has  a 
magnitude  given  by  r,q.  11.61  at  M*~{  at  a frequency  dd=  fa*  , for  that  magnitude 
of  & which  makes  it  equal  to  the  noise  of  /fn=D  at  the  same  frequency:  in  symbols, 
r (SKs)*)*  ")  7 

r,'(Kxs)j  i+  e"5^  j 


or 


- £ , />  -£ 


- 2(£l)y 


t.'  (m) 


(Mp  K'l C + 


Discussion:  When  the  noise  band  width  exceeds  the  frequency  of  the  acoustic  signal  to 


the  extent  that  then  and 0<7  For  this 

I,  0,7  To(^s) 


case 


with  a solution. 


Ktfx  - 11*5  a M 

Xe 


or 


_7 

We  choose  /Ie  X.  Zt/Q  Hence  the  threshold  signal  is 

.7 


_ ? 


i_i  .5,/0  = A.4X/0  H 

47T 
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In  contrast  when 


I * I > C 6 . if  we  choose 

3^l), fc>£*<f4=6.©ll  ,V<>  0 ‘/OS',  Thus 

T,  (Kxs)  - o-ioy  Io(Vxs) 


with  the  solution, 

^240,2,  of  xs  aj  Sx/o'^M 

*n 

A table  of  minimum  detectable  signals  is  presented  here,  based  on  the  approximation  that 

j.(Kxs)  a : i.e.  \ £TiXE  z &Jc  //o'7 

1 2.  7 47C  ** 


Table  I 


)(s(meter) 

(4B) 

5 

1.5  x 10’10 

6 

9.8  x 10-12 

7 

3.8  x 10-13 

_ /* 

Thus,  in  order  to  achieve  detection  of  the  particle  displacement  3 14/ /O  assoviated 
with  sea  state  1 the  noise  band  must  be  less  than  j/y  of  the  radian  frequency  of  the  sig- 
nal. At  100 HZ  we  therefore  require  A B 4 ( /rJ  ?fcC7T  * ; tfi  4^14- 

Hz. 


A different  criterion  of  signal  detectability  is  to  assume  that  the  discrete  spectrum  of  Eq.  III.  8 
at  (t)  ■ is  equal  to  the  continuous  spectrum  at  the  same  frequency.  Thus  the 

condition  becomes 

* r i l 


Kx 


s 

S 


< 


I + 


|4*"if 

\4&Vt/ 


= It 
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As  an  example  let  ^5*^  Then 

K2/s  s it  (0,33) 

Xs  = ts-t/o1  * 6.4//o**tf 

47T 

Shot  Noise 

We  assume  the  optical  signal  is  strong  enough  so  that  the  dominant  noise  process  in 
the  electronic  system  of  the  photodetector  is  shot  noise  (see  Oliver,  1961  “Signal-to-Noise 
ratios  in  photolectric  mixing”,  Proc.  I RE  (Correspondence)  49,  p.  1960).  The  rms  noise 
current  in  the  local  oscillator  heterodyne  system  is  then 

in  * J zetiufils)* 

(see  discussion  in  "Summary  of  Measured  Displacements”).  Under  the  assumption  that 
one  has 

Hence 


jl  , _ 34.  „ 

Now  choosing  Xg-S.  Uffv/O  H,  r)=  O.0S^z  4.^,  > /o  joule  sec,  CE » 3 y /Q g 

m/sec,  one  has 

s ft  y /0'7 

4 t //& **  *3>  /0? 

Also,  the  electronic  charge  is 

.If  _ . 

£>«  ).C  * /0 

Hence  j , _l/  * 

(,N  -r^l.  6 >/0  10  "g  /fa's 


^CE 
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If  the  system  of  detection  is  shot-noise  limited  then  this  electric  current  is  the  threshold 
of  detection.  In  the  general  case  of  an  unselected  laser  this  noise  threshold  is 


lh  = 5.09  */0 


in  which  all  quantities  are  in  n< s units. 

Minimum  Detectable  Signal  In  the  Presence  of  Shot  Noise  (Massey  1968,  loc.  cit) 

We  assume  here  that  the  only  noise  in  the  system  of  a local  oscillator  heterodyne  is 
shot  noise.  The  photodetector  electric  current  in  the  presence  of  an  acoustic  signal  is  then 


t(t)«  Tu,  * Ts  ♦ z J IL0 Is  (b >(<+*  + fa 

c 1 T 

in  which  Ti0*  _S  Itli  The  side  band"  structure  has  the  form 
X > “ * * 

in  = J Jo  (K* i)  c*  ^ T*>  (&*)  “*[(*>*+*> 

I 


“I  ; (Kxs)  &[(*-/»*)*]  I 

m-i 


Assuming  1^ « 1,  then  To  (Kx,)*!,  T,  (fa)*  faUjn 
Then  the  amplitude  of  current  in  one  side  band  is 


If  this  is  equated  to  the  shot  noise  current  the  minimum  detetestable  acoustic  amplitude  is 


Choosing  Xg 


(*'  Km 

-7 

S-  IAS  >10  } 


. (ICcl'-i  l*  • 

' OTT/ 

a -5*  5 

*)■  6.05,  *\s  i 6^/0  f Ce  r 3x/0  t it  u seen  that 


* 4X  x/o" 


,U/JL 

/ *p 

y ro 
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Scattering  Equation 

The  scattering  function  is  definded  by  the  relation 


(M  s-fr’1) 


in  which  1(8)  is  the  light  intensity  (watts)  in  the  8 -direction, E is  the  irradiance  (watts/  J 
is  the  volume  of  scattering  (M').  The  total  scattering  coefficient  for  ; 


non- 


polarized light)  is 


b=lTt^  (W“*) 


Now  let  tne  differential  scattering  volume  be  the  simple  geometric  figure  of  a cross  section 
A normal  to  the  beam  and  a depth  1.  The  product  is  the  power  of  the  incident  beam 
T*l . Hence  the  power  scattered  into  the  solid  angle  dO.  is 

A?s(e}«  aI(8)aIU&) 

The  incident  power  Is  related  to  the  laser  output  power 

_ o/T| 

* y*  e 

in  which  V is  the  attenuation  coefficient  and"R,  is  the  distance  travelled.  The  solid  angle 
Q is  that  intercepted  by  a lens  of  aperture  A 4,  at  distance  Kj  ^ i- 1 , 

Thus  the  power  received  over  AA  is  \ 

^2 
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Calculations 


From  Duntley  C963),  J.O.S.A.  53,  (214-233)  the  “average  ocean  water”  can  be  de- 
scribed as  having  an  attenuation  d m 0.05/fl  , and  a scattering  coefficient  ^(/Si  )-  6X/0  / 


meter  steredian. 


Assuming  :t  is  seen  that 


3^ . 6 »,crV0''o1W 


If  we  take  a receiving  lens  of  diameter  0,1  meter,  a depth  of  scattering  volume  at  least  a 
half-wavelength  of  acoustic  signal  at  lOOHl  (•  7.5  meters),  and  a working  range  of  30 


meters,  it  is  seen  that 


— _ x r TT 


* I,  H x/0 


If  “clear  water”  is  considered  then 


lered  then  ^ (&)%  ^ >t /t> 


0 , and  the  backscattered  light  is 


Modulation  Index  and  Doppler  Shift 

The  modulation  index  appears  in  Eq.  II.  57  as  the  symbol  KXs  in  which 

The  threshold  particle  displacement  at  |00HZ  is  3.4  H • Choosing  the 

o 

angle  of  observation  of  the  scattered  light  to  be  180  we  see  that  the  threshold  modula- 


tion index  is  calculated  to  be 


The  Doppler  shift  in  the  presence  of  an  acoustic  signal  is  given  by  Eq.  11.61  . Here 
it  is  seen  that  a constant  velocity  K#  shifts  the  frequency  of  the  carrier  by  an  amount 
■ The  acoustic  signal  shifts  the  frequency  be  the  amount  mfc)s, 
where  m is  the  number  of  the  side-band.  We  see  then  that  an  acoustic  modulation  of 


the  laser  carrier  does  not  appear  as  a true  Doppler,  namely  as  a frequency  shift  dependent 
on  velocity,  but  rather  as  a frequency  shift  due  to  phase  modulation,  whose  magnitude 
(of  shift)  is  independent  of  the  associated  modulation  index,  but  whose  amplitude  of 
power  spectrum  at  the  frequency  of  shift  does  depends  on  this  index.  The  detection 
method  is  not  that  of  an  laser  Doppler  velocimeter  measuring  a constant  velocity,  but 
rather  of  an  optical  heterodyne  measuring  sinusoidal  displacement.  Ultimately  the  de- 
tection process  rests  on  the  ability  to  measure  magnitude  of  spectrum  rather  than  the 
ability  to  measure  frequency  shift. 

Particle  Size,  Particle  Density,  Volume  Scattering  Function,  Attenuation  in  the  Ocean 
Suspended  material  in  the  ocean  which  is  retained  by  a 45  micron  filter  is  called 
particulate  matter.  The  amount  of  particulate  matter  (in  milligrams  per  liter)  ranges  from 
0.04  in  the  surface  water  of  the  North  Atlantic  to  18  in  the  region  of  coastal  waters.  An 
average  oceanic  total  is  0.8  to  2.5  (See  Table  II). 

Particle  size  distribution  in  various  ocean  bases  is  shown  in  Fig.  12(a)  in  which  is  a graph 
of  the  number  of  particles  per  cubic  centimeter  versus  the  diameter  of  the  particle.  A 
statistical  analysis  leads  to  a statement  that  the  mean  squared  particle  radius  is  about  14 
x 10'12  M2,  so  that  the  “mean  particle  diameter”  is  about  7 microns.  According  to  Mie 
theory  the  scattering  cross-section  of  a single  sphere  (^ ) whic^  is  much  larger  than  an 
optical  wavelength  is  twice  the  geometric  cross  section. 

Hence  ^ - 2Jt  x 14  x 10'12  - 8.8  x 10’11  M2 

Table  III  shows  the  vciume  scattering  function  of  pure  water  and  sea  water  as  a func- 
tion of  wavelength.  Fig.  12(b)  shows  the  attenuation  of  sea  water  as  a function  of  wave- 
length. 

Motion  of  the  Scattering  Prrticle 


, X'l\  A 


Table  !•  Amount  of  suspended  particulate  matter 


Depth  Suspension 

(m)  (mg/1) 


Reference 


Oceanic 


deep  water 


North  Atlantic  surface  water 
Oceanic  — 


Pacific,  coast 
Coastal 

(b)  organic  fraction 

Atlantic 
North  Atlantic 
North  Atlantic 
Central  Pacific  s 
Central  Pacific 

(c)  inorganic  fraction 

Atlantic, 

offshore 

Coastal 


deep  water 
surface  water 
deep  water 


0-05  Jacobs  and  Ewing, 

(average)  1969 

0 04-0- 15  Folger  and  Heczcn, 

1967 

0- 8-2-5  Lisitsyn,  1959 

(average) 

1- 6  Goldberg  >.t  al.,  l'W2 

6-0-18-0  Postnia,  1954 


004-0-17  Rfiey  d al.,  1965 

0 05-0-2  Gordon,  D.  C.,  1970a 

0 01-0  02  Gordon,  D.  C.,  1970a 

0-02  Gordon,  D.  C.,  1971 

0005  Gordon,  D.  C.,  1971 


0 05-1-0  Armstrong,  1958,  1905 

0-16-1-20  Armstrong,  1958,  1965 


10  10'  10‘ 
Diameter  (pm) 

Fig. (JftExnmpln  of  particles  size  dis- 
tribution!: 

Upper  icale: 

A Kullenberg,  Pacific  deep,  1953; 

|— | B run -Cot  tan,  Coulter  counter,  500  m 
depth,  Mediterranean,  1971. 

Lower  scale : 

• Gordon  D.C.,  microscope,  organic 
matter,  lurface.Atlantic,  1970; 

x Carder  et  a l..  Coulter  counter,  Pacific 
surface,  1971; 

A Jerlov,  microscope,  fiord,  1955; 

O Ocliakovsky,  microscope,  Mediterra- 
nean, 19GGo. 
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Table  H Volume  scattering  function  at  90°  and  total  scattering  coefficient  for  pure  water  and  sea  water  as  a 

function  of  the  wavelength. 


A(nm) 

350 

375 

400 

425 

450 

475 

500 

525 

550 

575 

600 

/UltHm-1) 

647 

2-80 

218 

1-73 

1-38 

112 

003 

0-78 

Pure  water 

6*(10~*  m-1) 

58- 1 

44-7 

34-9 

276 

22-2 

179 

140 

125 

109 

0M(1O-4  m->) 

841 

624 

4-72 

3-63 

2-84 

225 

1-80 

146 

121 

101 

0-88 

Pure  sea 
water 

6*(10~*  rrr‘) 

1345 

99-8 

755 

581 

454 

359 

28-8 

233 

19-3 

13  2 

141 

(S  = 35- 
39%.) 

* Computed  according  to  oq.  (II)  with  8 * 0*09  which  load*  to  6 =*  lt»  0 / 0 (90). 


fig  | jbAttenuation  curvet  in  the  near  ultraviolet  and  in  the  visible  part  of  the 
spectrum. 

A Lenoble  Saint  Guily  (1055),  path  length:  400  cm; 
x • • Ilulburt  (1034)  (1945),  path  length:  364  em; 

• Sullivan  (1903),  path  length:  132  cm; 

O — Clarke-Jamet  (1030),  path  length:  97  era  (Ceresm  lmed  tube); 

0 • - • James  Birgo  (1038),  path  length:  97  cm  (Silver  lined  tube). 

Total  scattering  coefficient  for  pure  water  and  pure  tea  water  as  a function  of 
wavelength,  according  to  Tablotlf. 

s'! 

<5 
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in  which  • If  the  period  of  oscillation  is  made  infinitely  long,  then  <2,X«. 

<<  ^ and  the  term  in  (i.e.,  the  inertial  increment  due  to  the  gross  mo- 

tion of  the  fluid)  becomes  negligible.  Under  these  conditions  the  magnitude  of  the  force 
exerted  by  the  fluid  on  the  sphere  is  Va lL  This  is  the  stokesian  force  resisting 

the  slow  descent  of  a particle  in  a viscous  fluid  under  the  action  of  gravity.  The  equation 
of  forced  motion  arising  from  this  resistive  force  is  therefore 

% „ 4x^a.U 

in  which  g V is  the  fluid  dynamic  viscosity.  Thus,  for  particles  of  diameterT) , 

ei*t 

ft 

SinceTJ  is  the  differential  motion  of  the  water,  i.e.,  , we  can  assume  Ijj  to 

be  the  acoustic  excitation  and  write  Iff.*  Uafcpjrf  At  the  steady  state  frequency 
it)  the  velocity  (Js  reduces  to 

u - Uo 
4*U/ 

% v 

This  equation  defines  the  properties  of  a low-pass  filter  with  a cut-off  frequency  * <A. 

As  long  as  the  motion  of  the  colloidal  particle  will  faithfully  follow  the  motion 

of  the  acoustic  wave.  To  calculate  the  magnitude  of  the  cut-off  frequency,  we  note  that  ^ 

the  dynamic  viscosity  of  water  at  20  C^»n  centipose'(  *«  0*01  5Cc./c*y,*  = 10  • 

Assuming  p>u  tide  sizes  of  order  7 x 10’®  meter,  one  calculate*  jC  * Oio  *0  /fs  • 

If  the  density  of  the  colloidal  particle  is  the  same  as  that  of  water,  the  cut-off  frequency 

^ Ill'S4*  5 8.S  kHz. 
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By  choosing  (t)<,<(^,  one  sees  that  the  difference  be- 
tween the  motion  of  the  colloidal  particie  and  that  of  the  acoustic  wave  is  negligible. 

It  appears  from  the  above  argument  that  for  (low-frequency)  periodic  motion,  which 
is  long  relative  to  the  optical  process  duration,  the  forces  involved  are  purely  resistive, 
and  the  particle  motion  is  not  sensibly  different  from  the  wave  motion.  However,  when 
the  wave  motion  is  complex  due  to  presence  of  soft  reflective  walls,  diffraction,  etc.,  the 
forces  acting  on  a colloidal  particle  are  along  more  than  one  coordinate.  The  motion  of 
the  particle  is  then  oval,  or  elliptical,  and  the  equations  written  above  no  longer  hold. 

It  could  then  be  said  that  in  complex  sound  fields,  it  can  hardly  be  expected  that  a colloi- 
dal particle  will  faithfully  follow  the  motion  of  the  acousitic  wave  since  it  would  then  be 
required  to  be  nearly  indistinguishable  from  the  medium  itself^ assumed  free. 

Brownian  Motion 

According  to  the  theory  of  Brownian  motion,  the  mean  of  the  square  of  the  distance 
travelled  by  a particle  in  a fluid  during  a time  t is  given  by  the  equation 

hF  = GT>t 

(See  Lundau,  Lifshitz  “Fluid  Mechanics**,  p.  227ff).  HereJ)  is  the  diffusion  coefficient 
(dimensions  M S j.  For  spherical  particles  of  radius  Q,  diffusing  slowly  in  a me- 
dium of  dynamic  viscosity  U , it  is  known  thatD  . Thus,  the  mean  distance 
travelled  in  timet  U 


It  is  important  to  es  mate  the  time  required  for  the  particle  to  mow?  a "decorrelation 
distance’’.  Arbitrarily  (but  reasonably),  this  distance  is  taken  to  be  a quarter  wavelength 
of  the  laser  light  .Let  7^  * 3#60  A t 10  RM  5 , & * 3 1 10  M } Boltzman’s 
constant  ■ 1.38  x 10*23  S K*  1 "f=  300°  K Then  the  decorrelation  time  is 


--(}■  U '.'ll  K ’/o'**  3>'oj_  , 

' A I O O - , -v  .A^ 


I • **  5 t to 


j 
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The  “decorrelation  speed”  is  v A1  Ut  or 


Q . (3-  ) 


a.lA/O'2- 


-4.4  y/0*  Ms' 


Thus,  Q is  of  the  order  of  10  J/tnm  /sec* 

The  possibility  of  decorrelation  of  the  motion  due  to  the  acoustic  wave  by  the  Brown* 
ian  motion  of  the  suspended  particles  must  be  seriously  considered,  if  the  duration  of  the 
process  required  to  sample  the  wave  exceeds  0.02  sec.  At  100  Hz  this  allows  about  2 com- 
plete cycles  to  be  sampled.  However,  at  50  Hz  only  1 cycle  can  be  sampled.  Thus,  there 

if 

is  a threshold  frequency  for  doppler  detection  of  particle  velocity,  estimated  here  at  50  Hz. 
Spectral  Broadening  Due  to  Brownian  Motion 

The  probability  of  finding  a Brownian  particle  in  the  distance  interval  Aand  fL+dft  is 
proportional  to  Hyp (-  /i* /Ot)  ( See  Landau,  Liftshitz,  p.  227).  For  two-way  travel  the 
phase  change  A<j>  due  to  Brownian  motion  is  Hence,  the  probability  of 

finding  the  motion  of  the  particle  in  phase  A$  is  proportional  to  exp CA4>)V<47ll^vt)}  = 
JUbi  I where B = ^5-  . In  accordance  with  the  mathematical  model  sketched 

*' 1 t * 7t 

earlier,  the  spectrum  of  the  first  sideband  due  to  acoustic  modulation  is  1/2  of  this  quantity. 
Thus,  the  spectral  broadening  due  to  Brownian  Motion  is 

A „ li  . ii  ir  4T 
1 3 ty* 

For  a laser  wavelength  in  water  of  b,  84  X iQ  M , It*  (Xs  3 x/0  M 

mij  | J * 

Wr  J r»  300V  the  spectral  broadening  is  155  Hz. 

Multiple  Scattering 

The  transport  equation  can  be  used  to  model  multiple  scattering.  Let  % (h , / -tjcl/c.  d ¥ 


d A.  & V represent  the  rate  at  which  photons  are  introduced  into(£ , j by  the  laser. 

Then  the  most  general  time  dependent  transport  equation  is 

♦This  conclusion  is  restricted  to  the  case  of  s single  scsttering  particle  whose  Brownian  motion  disturbs  the  scoustic 
particle  velocity  over  short  times,  snd  over  an  assumed  chars cteriatic  length  In  the  multiple  particle  case 

the  Brownian  motion  is  random.  Over  long  (enough)  times  it  only  adds  noise  to  the  detection  process  but  does  not 
set  threshold  frequencies.  The  characteristic  length  is  then  not  significant. 
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U + v-  y$  + v<r  (h,v.t)  - f i/V  $ c /V  dfv'dt  +£ 
dt  ~ ~ " J 

This  is  a differential-integral  equation  in  the  unknown  function  £ and  known  function  0*. 

It  is  too  intractable  to  solve.  Simplification  is  usually  accomplished  by  considering 
only  the  axisymmetric  case  in  which  /i  is  replaced  by  ft  ( 2 being  the  direction  of  prop- 
agation, jJ  the  polar  angle  to  this  direction),  and  considering  only  uniform  velocity^  so 
that  $ is  independent  of  velocity.  Then,  the  transport  equation  reduces  to 

+ £)  5^)’  ft  frwS)*  WW 


where  (T,(^  are  the  scattering  and  absorption  cross-sectioiu^respectively,  and  L=  diffusion 
length.  The  solution  is  assumed  to  have  the  form  . If  multiple  scattering  is  ne- 

glected, then,  approximately 


in  which  I,  is  number  of  photons  simply  scattered  from  the  initial  set  j^.  The  first  term  on 
theM.S.  is  the  number  of  singly  scattered  photons  that  are  lost  by  absorption  and  scattering. 
By  changing  subscript  0 to  1 and  subscript  1 to  2,  this  same  formula  can  be  used  to  estimate 
the  number  of  photon  that  are  doubly  scattered.  A repetition  of  this  bookkeeping  process 
leads  to  the  GE  conclusion  that  27%  of  the  received  laser  light  will  be  multiply  scattered 
and  73%  simply  scattered  (See  GE  Report  973-SH-347,  March  1970,  R.  M.  Ameigh,  et  al, 

p.  28). 

This  conclusion  on  the  partition  of  the  incident  photons  into  single  and  multiple  scat- 
tered photons  is  based  on  neglect  of  the  integral  in  the  transport  equation  which  accounts 

I 

for  the  transfer  of  photons  from  direction  to  direction  jl  ; i.e..  the  volume  effect  is 
neglected  leaving  only  the  linear  effect.  A true  account  of  multiple  scattering  must  rely  on 
solving  the  integral  equation  per  se. 

While  single  scattering  accounts  for  observed  scattering  effects  over  short  path  lengths, 
it  is  not  valid  over  long  path  lengths.  In  the  latter  region  multiple  scattering  is  the  dominant 
factor. 
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Part  V 


The  Effect  of  Medium  Distortions  On  the  S/N  Ratio 
of  an  Optical  Heterodyne  System 


Fried  (Proc.  IEEE  55,  p.  57,  1967)  has  developed  a theoretical  formula  for  the  limit 
of  improvement  of  the  S/N  ratio  in  the  presence  of  a distorting  medium.  We  review  here 
the  conclusions  most  relevant  to  our  study. 

Let  the  photodetector  be  modeled  as  a circular  aperture  of  c’iameterD  . By  use  of 
an  appropriate  weighting  function  for  the  contribution  of  each  elementary  area 
of  this  circle  the  total  photocurrent  is  derived  to  be 


i‘  { J<«i  %.)  <|  (&,+ (E^E.) . T Vr)  (Es +E  jfc  ♦&) 

D 

(V.l) 

in  which  is  the  quantum  efficiency  of  the  detector,  E4  is  the  complex  scattered  signal 
field  and  E0  the  complex  local  oscillator  field  of  the  optical  heterodyne.  On  the  assump- 
tion that  Eo  > it  is  seen  that  the  information-carrying  part  of  the  photocurrent  is 


ts=  7L?  |€o\ \ES\  Ort  (2TC^t  +A<t>) 


(V.2) 

in  which  s ■£$-  ^s"  , anci  IE.I.N  are  the  amplitudes  of  the  local 

oscillator  beam,  and  the  scattered  beam  respectively. 

Let^  be  the  gain  associated  with  the  current  amplifier,  and  "Ri  resistance  of  the  load 
in  the  detection  circuit.  Then  the  time  average  of  the  signal  in  the  photodetector  is 


(V.3) 
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C»..  •j.  J&t  sfc. 


For  circuit  noise  fluctuation  current  ^ («  shot  noise),  the  amplified  noise  is  (q  and 
the  noise  power  per  unit  bandwidth  is 

N = et^f)  15  [E^R 

(V.4) 

in  which  6 is  the  electronic  charge.  Hence  the  S/N  ratio  per  unit  bandwidth  is 

" e s 

(V.5) 

x A 

Thus  S/N  is  proportional  to  the  total  signal  power  collected  (i.e.  K "D  \HS|  ) and  pro- 
portional to  the  quantum  efficiency  % of  the  detector. 

When  the  wave  field  scattered  from  the  medium  is  random  one  replaces  the  non- 
random  phase  4^  by  thej-andom  phase  (*  ) . and  the  non-random  amplitude  (Ej  j by 
the  random  amplitude^!  *xp  [ £ ( £)1  , in  wliich  is  the  ms  value  of  E (gj , and 

Hi)  is  the  log  amplitude  i.e.  ji(t  '^lie  •Ufnal  current  I5  is  then  de- 

rived to  be 

■s.  r,  iE.lIiU? 

(V.6) 

The  information  bearing  signal  is 

s • (tj'il’it  %\UQ)R  { 5 h^1*- -4f 

(V.7) 

‘S.  WE.feiQ)R 

Ln  the  last  step  the  time  average  has  been  dropped  because:  (1)  the  average  over  <fae  vanishes; 
(2)  the  spectra  of  (x')  and  \ (&)  are  within  the  information  bandwidth  of  the  detector, 
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(x)  and  X (jl)  therefore  being  virtually  constant  over  the  averaging  period.  Since S u 
random  one  can  form  the  ensemble  average  In  terms  of  a change  of  variable,  ftsjt*)!1, 
this  average  is 

< s >=  l (n  jjtotc'  W<4'*  iil)W(iA-iti)  [-i&ki] 

(V.8) 

where  ia  the  wave-structure  function  defined  as  the  sum  of  the  structure  function 

for  phase  (k)  , and  amplitude  (fi.)  , i.e. 

M - <t+i«-W3> 

3,  (a}  --  <&)-W> 

A = |£\  (V.9) 

By  choosing  W (£  ) properly  it  can  be  shown  that 

<S>-  u (,||e.lp4€i)*  ^>A«h.k'.(A^*|,[-iS61(A)] 

o 

KM-i  [*'-»■!:&)- a 


(V.10) 


The  shot  noise  in  the  detector  system  is  the  same  as  Eq.  V.4.  Thus  the  S/M  ratio  is 

<s>/ti  * 4.  i IfiL*  £ ftdx  K>M [- { 3K*)] 

" = (v.n) 

%’  i>  1 

(V.1*> 
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The  symbol  /i»  is  introduced  for  the  following  reason.  A graph  of  V*  (^h«^vs.  % 


in  log-log  coordinates  has  the  following  form, 


FIG).  115  (V.13) 

This  graph  is  based  on  several  assumptions:  (1)  both  (i  ) and  ) obey  Gaussian 
statistics  (2)  the  turbulence  in  the  fluid  is  isotropic  and  homogeneous  (3)  the  Kolmogoroff 
theory  of  turbulence  applies  (i.e.^S^A  , ) , $(h)x 


theory  of  turbulence  applies  (i.e.^S=i^A  , A» 

l_0  j (^structure  constant  defined  by  Tatarski,  “Propagation  in  Turbulent  Media” 
1961  Eq.  3.50).  From  the  graph  it  is  surmised  that  / N reaches  an  asymptotic 
limit  as  the  aperture”])  is  increased,  and  that  ho  is  the  smallest  D for  which  this  limit  is 
available.  Thus  increasing  the  aperture  beyond  Ae  does  not  increase  the  ^5  >/N  ratio, 
given  the  existence  of  wave  front  distortion.  (See  next  section  for  critical  comment.) 

Fried  has  calculated  A*  for  many  cases  of  distortion  in  atmospheric  turbulence  (see  e 
his  Fig.  8).  We  note  here  that  for  the  case:  (1)  electromagnetic  wavelength  X*  OiSjsrt  ■ SOOflA 
(2)  zenith  angle (3)  altitude  of  receiver  H a 1 kilometers  (4)  daytime J the  magni- 
tude of  ho  is  approximately  0.12  meters  (say  5 inches).  However  in  the  case  of  under- 
water turbulence  the  calculation  of  A6  is  a subject  of  further  investigation. 

Recently,  considerable  progress  has  been  made  in  the  mathematical  theory  of  the 
propagation  of  laser  beams  in  strongly  turbulent  media.  A summery  of  developments  is 
provided  by  Prokhorov  et  al.  Proc.  IEEE  63,  p.  790  (1975),  where  the  subject  of  saturation 
of  intensity  of  scintillations  is  treated  in  great  depth. 
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Limits  on  VAS  Imposed  By  The  Inhomogeneities  of  the  Medium 

The  VAS  is  equivalent  to  an  array  of  remote  sensors  located  in  an  inhomogeneous 
medium.  The  medium  imposes  limits  on  its  performance.  Fried  (see  previous  section)  has 
discussed  the  effect  of  the  medium  on  the  photodetector  performance  for  a single  sensor. 

He  concludes  that  the  S/N  ratio  cannot  be  \nproved  by  increasing  the  receiver  aperture, 
once  a critical  aperture  size  is  reached.  This  conclusion  must  be  reviewed  in  light  of  known 
statistical  antenna  theory.  Since  this  theory  is  very  extensive  we  focus  here  attention  only 
on  one  key  theoretical  point,  namely  on  the  maximum  attainable  directivity  of  a statisti- 
cally perturbed  array. 

Shifrin  (“Statistical  Antenna  Theory’  Golem  Press,  1971)  provides  a convenient  sum- 
mary. Let  T)0(*  iL  be  the  directivity  of  a continuous  line  of  length  i_  located  on 
the  x-coordinate  axis  receiving  a (monochromatic)  wavelength  A . Also  let  4(x)  be  a ran- 
dom phase  of  (normalized)  position  X * A i/L  . It  is  first  assumed  that  4(X)  >s  normally 
distributed  with  zero  mean,  variance  oL  (independent  of  coordinate  .:),  and  correlation 
function  • The  mathematical  form  of  ^ is  assumed  to  be  of  the  form 


( i |- 


,a 

( *-  x'l 


1 


7 


L4  - ^4 


(.i) 


in  which  is  the  normalized  correlation  length  for  phase,  and  (p  is  the  non-normalized 
correlation  length  for  phase.  The  mean  directive  gain  3)  is  given  by 


D 


■ - W 

« if 


dx 


(.2) 


Let  ^ be  the  far  field  pattern  angle,  measured  from  the  normal  to  the  center  of  the  line 
array,  and  let  the  generalized  angle  be  defined  by  Ip  • 

When  discussing  directivity  of  an  array  that  sees  a wavefront  with  (Gaussian)  phase  correla- 
tion A it  is  convenient  to  define  a function  1 ip|')  in  which 
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I (U.V.-'Yi')  * 3lU,  <“*'&>) 

Jr  P (-3) 

In  terms  of  these  quantities,  it  is  seen  by  expanding  the  exponential  that  the  mean 
directive  gain  is 

5-a»oe“  [ i+±  T(  «*0,0\1 

I r' 

- . u '•4) 

A plot  of  .D  vs  -M  /a  for  various  values  of  phase  variance  06  and  correlation  distance 
is  shown  below: 

j Z A d*o 

^ /c(*l 

S^y5  FIQ  14 

/y  / /ot“°  (Shifrin,  Fig.  12.3) 


J.  ^ 4,  * \ ;e ^ 

Here,  SO  H curves),  s |oqK  (^)ro'lCen  curve).  These  curves  can  be  discussed 

in  terms  of  the  parameter  /\Jq[  , i.e.  the  ratio  of  normalized  correlation  length  of 

phase  to  the  standard  deviation  of  phase.  Two  cases  are  of  interest. 

Case  I.  Let  , that  is,  the  antenna  dimension  is  much  smaller  than  the  correl- 
ation length  . Then  choose  Lot»l . It  is  seen  that  . This  means  the  fluctu- 

ations of  the  medium  have  no  effect  on  antenna  gain.  Now  choose  Let  , meaning  that 
the  standard  deviation  tor  phase  is  very  large.  Then 


»■& 


i 

i 
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in  which  "V  is  the  mean  value  of  the  square  of  the  (random)  angle  of  arrival  (see 
Fig.i5). 


L 

It  is  seen  thac  antenna  gain  X)  does  not  depend  on  the  length  of  the  antenna.  This  is  the 
region  of  directive  gain  saturation.  However  as  the  length  of  the  antenna  increases  the 
gain  saturation  disappears  and  the  mean  directive  gain  increases  with  antenna  length  with* 
out  limit. 

Case  II.  Let  1 (that  is,  the  antenna  dimensions  are  much  larger  than  the 

correlation  length^  ).  Once  again,  for  sufficiently  large  phase  fluctuations  there  is  a 
region  whereD  is  independent  of  length  of  the  antenna.  This  is  the  region  of  directive 
gain  saturation.  When  the  antenna  length  is  increased  beyond  this  region  the  gainD  in- 
creases with  length  of  the  array  beyond  all  limit.  ' 

Summarizing:  An  inhomogeneous  medium  limits  the  performance  of  a linear  array 
only  over  a limited  region  (region  of  gain  saturation)  of  array  lengths.  Beyond  this  region 
the  directive  gain  of  the  antenna  increases  with  length  beyond  all  limit.  The  limits  of  S/N 
discussed  by  Fried  in  the  previous  section  are  seen  to  apply  only  over  limited  detector 
apertures.  If  the  aperture  is  made  sufficiently  large  (according  to  Shifrin), there  is  always  a 
directivity  gain  for  a fixed  standard  deviation  of  phcse  fluctuation.  Thus  S/N  can  be  in- 
creased by  increasing  the  aperture  provided  the  aperture  is  large  enough. 

♦ 
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Part  VI 


Signal  Processing 


We  take  a large  two-dimensional  seismic  array  to  be  roughly  analogous  to  the  projected 
VAS  array,  and  discuss  seismic  signal  processing  in  a brief  review. 

Seismic  signals  are  always  mixed  with  seismic  noise.  Several  methods  of  processing  to 
reject  noise  are  known.  In  the  first  method  (Burg,  Geophys.  29,  p.  693-713,  1964)  both 
signal  and  noise  are  assumed  to  be  stationary  multidimensional  random  processes  with 
known  cross-correlation  functions.  A linear  (Wiener)  filter  is  designed  which  produces  an 
output  that  is  a minimum  mean-square-enor  version  of  the  signal.  In  a second  method 
(Claerbout,  Geophys.  29  p.  197,  1964)  the  known  noise  correlation  function  is  used  to 
design  a linear  processor  for  the  seismometer  outputs  which  provides  a minimum-mean- 
squared-error  pred;ction  of  the  noise  over  a short  interval  ahead.  This  prediction  is  sub- 
tracted from  the  actual  seismometer  output,  thus  greatly  reducing  the  noise  level.  How- 
ever the  signal  waveform  is  distorted.  In  the  third  method  (Capon  et  al.  Proc.  IEEE  65, 
p.  192,  1967)  the  seismic  noise  is  assumed  to  be  a zero-mean,  time  stationary  multi- 
dimensional random  process.  The  signal  is  assumed  to  be  a single  plane  wave  propagating 
in  a homogeneous,  linear,  nondisperslve  medium,  but  is  an  otherwise  unknown  time  func- 
tion. The  signal  processing  therefore  takes  on  the  character  of  generating  an  estimate  of 
an  unknown  time  function  and  uses  the  theory  of  maximum  likelihood  (or  minimum- 
variance)  unbiased  estimator  approach.  The  essential  feature  of  this  approach  is  to  design 
a noise-rejection  filter  based  on  a matrix  of  filter  weighting  coefficients  which  are  obtained 
by  use  of  the  calculus  of  variations  on  a system  of  equations.  Briefly,  let  there  be  a total 
of?  sensors  in  the  seismic  array,  and  let  the  noise  in  each  sensor  constitute  one  component 
of  a multidimensional  random  process,  Nj  - t(  M ,(0  * N,(t)  ♦ • • • * 4- 1** 
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Assume  further  that  each  N J is  Gaussian,  so  that  the  total  N has  a multidimensional 
Gaussian  distribution  with  zero  means.  The  covariance  matrix  of  N is  defined  as 

K'h  (*„M  Kp(t,A)  \ 


with 


(t'>M  (MO  j 

K(,  (MO* g C ^ (<*>] 


where?  is  the  expected  value. 

Now  let  the  time  function  of  the  incoming  wave  be  sampled  A V+l  times  from  -V  to 
li  , and  designate  two  distinct  samples  by  integers*, ft.  Then  the  minimum-variance  un- 
biased estimator  has  weighting  coefficients  6^  ,^=1,  2...  "P  , determined  by  solving  the 
matrix  of  equations  ^ y 

Z Z.  (/m'I*)  Kfk  (*"VJ  + \m~° 

in  which  '^^are  1 Lagrangian  multipliers  chosen  to  satisy  the  constraints  that 
■P 

ZL.  6,  (d*  U)  - t ten. 

pi  6 ' 


Thus  knowing  , and  the  constraints,  one  finds  a set  of  "P  coefficients  , 

1,2  ...T.  These  coefficients  form  the  basis  of  constructing  a digital  filter  which 
(by  its  construction)  gives  the  optimum  estimate  of  the  unknown  time  function  con- 
stituting the  signs],  simultaneously  rejecting  noise.  The  synthesis  of  the  digital  (two 
sided)  filter  is  carried  out  in  the  frequency  domain.  It  is  noted  that  the  estimator  thus 
found  is  also  the  maximum-likelihood  estimator. 


Diversity 

Diversity  is  a signal  processing  technique  aimed  at  improving  the  reliability  of 
reception  of  signals  that  are  subject  to  fading  in  the  presence  of  random  noise.  To 
discuss  diversity  one  requires  a brief  summary  of  the  nature  of  a fluctuating  channel. 
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We  assume  that  there  exists  in  a communications  channel  many  “copies”  of  a 
signal,  travelling  different  paths.  Each  copy  has  an  identifiable  amplitude  and  phase 
which  vary  slowly  with  time.  At  a particular  receiving  point  the  instantaneous  time 
signature  of  arrival  of  signals  is  the  sum  of  two  or  more  sinusoids  with  vary  amplitudes 
and  phase.  This  sum  shows  fluctuations  in  signal  strength  called  multipath  interference 
fading.  By  choosing  many  different  channels  one  can  obtain  many  time  signatures  which 
differ  noticeably  from  each  other,  some  fading,  some  flaring  up,  etc.  If  these  signatures  (consisting 
of  signals  plus  noise)  are  widely  enough  spaced  in  location  (or  time,  or  frequency)  as  to 
be  statistically  independent,  one  can  combine  them  in  an  appropriate  manner  so  as  to 
obtain  better  or  more  reliable  reception  of  the  “message”.  Such  combinations  form  the 
basis  of  diversity  methods  of  signal  enhancement.  It  is  to  be  emphasized  that  enhance- 
ment is  possible  only  when  there  is  a random  character  to  the  received  signal  itself,  and  when 
diverse  copies  of  the  signal  can  be  obtained,  or  generated,  from  a single  original. 

A convenient  and  much  used  model  of  a fluctuating  channel  is  the  "Rayleigh- 
fading”  model.  In  this  model  a transmitter  projects  a signal  at  frequency  fcle  which 
arrives  at  the  receiver  in  the  form 

e(t )*  A(i)  CM 

The  random  amplitude  /4  (t)  is  Rayleigh  distributed,  with  probability  density 


2 


in  which  (Lis  the  J A*  in  an  (appropriate)  interval  T . The  random  phase  is 
usually  assumed  to  be  uniformly  distributed  between  0 and  27F  with  a value  IT/a,  . 
Thus  i(t)  is  a narrow  b'1  ,d  gaussian  process  with  zero  mean  and  variance  OJ  = 0-  /i  . 
If  <*CO  i*  expressed  in  quadrature  form 

e(t)*-  e«(Otei  (0ct  + €i(i) 


then  0S  , are  independent  gaussian  random  variables. 
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In  these  (see  D.G.  Brennan,  MIT  Lincoln  Lab.)  probability  densities,  curve  0~  is  «-Rayleigh 
channel,  curve  b shows  mild  fading,  curve  C shows  deep,  frequent  and  persistent  fading. 

We  assume  now  that  a number  of  copies  of  the  signal  (or  “message”)  are  available. 
These  are  to  be  sent  from  source  to  receiver.  Upon  arrival  each  copy  is  assumed  to  show 
fading  in  the  presence  of  noise  of  the  transmission  path  itrelf  or  of  the  receiver.  To  apply 
the  technique  of  diversity  signal  processing  the  fading  patterns  of  the  original  copy  must 
be  statistically  independent.  To  assure  this  several  methods  of  transmission  are  used:  (1) 
copies  of  the  signal  are  sent  in  time  succession  as  a sequence  of  pulses.  The  signal  is  said 
then  to  be  time  diversified.  (2)  copies  of  the  signal  are  transmitted  on  different  carrier 
frequencies.  The  signal  is  said  to  be  frequency  diversified.  (3)  Copies  of  the  signal  are 
received  from  different  spatial  locations  in  the  medium.  The  signal  is  then  said  to  be  spa- 
tially diversified.  In  all  cases  statistically  independent  noise  is  assumed  to  comipt  the  signal. 

When  K copies  of  the  received  signal  are  available  they  can  be  signal  processed  in 
various  ways.  A simple  technique  is  coherent  processing.  In  this  procedure  the  1C  received 
signals  are  added,  and  then  squared  to  form  the  power.  The  signal  power  is  proportional 
to  K and  the  noise  power  is  proportional  to  C Hence  the  improvement  in  S/N  ratio 
is  equal  to  the  number  ^ of  independent  copies  available.  Such  processing  however  assumes 
the  quality  (good  or  poor)  of  all  signals  is  identical.  A more  general  type  of  processing  is 
to  weight  each  available  copy  (■  by  a weighting  factor  C+.  and  then  add  them  together 
to  form  (ty,  where 
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< 

7(0-  Z G+  {k b) 

K-i 

Different  choices  of  Q^lead  to  different  processing  techniques.  If  only  one  ^ is  different 
from  zero  at  any  time,  one  has  switching  techniques.  Trying  each  one  at  a time  is 
called  scanning  diversity.  Trying  a group  of  {^simultaneously  and  selecting  the  best  is 
called  optimal  selection  diversity.  If  all  the  are  used,  and  given  equal  magnitudes  the 
processing  is  called  equal-gain  combining.  If  the  flip’s  are  adjusted  according  to  the  quality 
of  ) so  as  to  yield  the  maximum  S/N  the  result  is  called  maximal-ratio  combining. 

These  different  procedures  can  he  calculated  to  show  the  S/N  gain  of  diversity  over 
signal  channel  reception.  A set  of  charts  have  been  prepared  by  Brennan  (1959,  Proc. 

IRE  47  p.  1075-1102)  for  the  cases  of  two-channel,  four-channel  diversity,  and  theoretical 
gain  for  a large  number  of  channels.  For  example,  a typical  case  shows  that  ten  (fading) 
copies  can  be  made  to  deliver  a 10  dB  improvement  over  a single  copy  when  the  technique 
of  maximal-ratio  combining  is  used. 
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Conclusion  to  the  Report 

We  have  presented  above  a comprehensive  review  of  all  factors  pertinent  to  the  detec- 
tion of  acoustic  phe  lomena  by  use  of  a laser  heterodyne  detector.  The  principal  problem 
in  proposed  Navy  upplitV-ions  of  this  technique  (as  replacement  of  standard  Navy  hydro- 
phones) is  the  very  small  signal  to  noise  ratio  in  a single  beam.  This  limits  detection 
threshold  in  noisy  environments  to  (a  few)  orders  of  magnitude  greater  than  the  thresholds 
currently  available  to  the  Navy.  However,  onsiderable  improvement  in  S/N  ratio  is  fore- 
seen in  the  proposed  use  of  multiple  beams  (i.e.,  “diversity”).  The  gain  is  anticipated  to  be 
proportional  to  the  number  of  beams  employed. 
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Appendix  A 


Analysis  of  Two  Beam  Interference 

Two  point  sources  of  light  separated  a distance  d in  a dark  screen  superimpose  on  an 
observing  screen  at  distance^ . The  total  field  at  point  ?(*)  in  the  observing  screen  is 
Ep*  ~ exp 

If  we  assume:  _ t- . (A-1) 

(t  } F,  * Ea  * tt 
'2.)  4>(^1)  * 4^0  * 1 

(3 ) |Aa I - Ar  > AT « ml,  Ar«|r,l 

^ >2E;CM(St  i^  (A.2) 

Now  let  the  angle  of  intersection  of  the  two  beam  be  6 . Then 

A.  » IfiU.  | (A.3) 

The  intensity  at"P  is  XL,  given  by 

i„  - i fi-.E 

•’  = Ctf*  E*[»4  X ( <a.5) 

" 1 £ <0  Et’  [h 

in  which  Ko  is  the  dielectric  constant,  and  c * £ £ 

i 

If  the  sources  have  finite  width  QQ  the  projection  of  which  on  to  a plane  normal  to 
OP  is  ^ and  if  of  is  the  angle  enclosed  by  the  rays  from  0 and  Q then  a path  difference 
is  introduced,  as  measured  form  origin  0 at  the  center  of  the  width.  The  phase  difference 
is  Kt^.Eq.  A.5  constitutes  the  basic  formulation  appearing  in  George  and  Lumley’s 
model,  Part  II  of  this  report 
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Appendix  B 


Dimensionless  Wavenumbers  and  Optimum  Scattering  Volume  (George  ar.d  Lumley  loc.  cit). 

The  symbol  Y|  is  the  Kolmogorov  microscale  ( (**%> ) ) and  has  the  dimensions 

of  length.  The  symbol  ot  is  a frequency,  appearing  in  the  spectra  for  velocity,  phase  etc. 
The  symbol  VI  is  the  mean  flow  velocity.  The  symbol  I?,  is  the  dimensionless  wavenumber 
for  which  the  ratio  of  turbulence  to  ambiguity  spectrum  is  unity.  By  definition 

k.  ■ <».!> 

Thus  if  R„m|.0  the  metui  velocity  u,  is  equal  to  the  velocity  of  the  fluid  of  displacement 


equal  to  the  Kolmogorov  microscale  multiple  by  frequency  &o. 

H 


The  symbol  ^ » (p'V.)  where  is  the  size  of  the  scattering  volume  in  the 
x-direction.  The  symbol  ^ e ^ is  thus  a ratio  of  the  Kolmogorov  microscale  to  the 
“x -component  of  scattering  volume."  The  statement  that  kf  is  optimum  when 


1.27 


./ JL\ 

(j9  0",  /fr 


means 


o ^ ^ 


(<r.L, . Jli  _ 

v 1 .37  41 


in  which  the  Reynolds  number  is 


zn  TC 
V(S0 


55 t ~ W tl* 


(B.2) 


(B.3) 


(B.4) 


The  symbol  where  is  the  size  of  the  scattering  volume  perpendic- 

ular to  the  flow,  but  in  the  plane  of  0",  For  small  6 , 

r 1,  aM.  AhH.  ‘b-5» 

Thus, 

$ . zli*.  $ » ^ 


f 


(Um  ■ 


y.  i.37  I 


The  optimum  volume  of  scattering  should  be  determined  by  this  formula 


(B.6) 
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Appendix  C 


Optic- Acoustic  Interactions 


Optic-acoustic  interactions  are  reviewed  in  Appendices  C,  D and  E. 

We  consider  first  a nonstationary  medium  in  which  the  time-rate  of  change  of  the 
index  of  refraction  is  much  slower  than  the  frequency  of  the  electric  field.  In  symbols 

T where /n  is  the  index  of  refraction  andT  is  the  period  of  oscillation 

of  the  electric  field.  Now  let  the  change  in  index  (s  m,}  be  caused  by  the  passage  of  a 
soundwave,  i.e.  caused  by  a change  in  mass  density.  Thus  when  the  sound  wave  is  present 
the  instantaneous  index  of  refraction  (-  "*)  is  the  sum  of  the  unperturbed  value  /ho  and 
the  perturbation  m . In  the  absence  of  true  sources  the  propagation  of  the  electric  vector 
E is  governed  (to  quantities  of  first  order  in  ^i]by  the  linear  wave  equation 

)'  £|_  - 0 . or  <7  £ - a g 
- <H*  ^ d at* 

in  which  terms  of  order  /Mohave  been  omitted.  The  perturbation  in  the  index  of  refrac- 
tion /r)t  is  proportional  to  the  perturbation  in  fluid  density  p,  t l-t 

f, 

7Z,  ' £ (C-2) 

thus  we  write  the  wave  equation  for  E in  the  form 

* p _ i d § «J  f i I 2>  f 


- Cc 


The  transit  of  the  electric  vector  in  the  medium  creates  a mechanical  body  force  rl 
(dimensions:  N M ) which  acts  as  a source  of  acoustic  pressure  j,  the  propagation  of 


which  is  governed  by  the  equation 


vV  JL  «-£-  7.  F, 

r c,*  at* 


Now  'die  acoustic  density  is  proportional  to  the  acoustic  pressure,  and  the  mechanical 
stress  is  proportional  to  the  electric  field.  We  write  t f E --  E • 

Thus, 
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(CA] 

^ ai*  - 

Taken  together  we  see  that  the  waves  of  density  and  electric  field  form  a coupled  system 
given  by  Eqs.  (Cl)  and  (C2).  Returning  to  Eq.  (Cl)  we  'et  P,  - p,  ( X ,4)  and  seek  a 
solution  E (y/t)  by  Fourier  transformation  of  the  time  coordinate  with  zero  iniual  condi- 
tions, 

(vW)  E (»,»)  =(|JP,  (>,U)  * (-|)EM=  «?(«,*>)  (C61 

in  which  * signifies  the  convolution  opera-ion.  'Hie  symbol  Q (*,£*>) describes  a 
(fictitious)  source.  The  electric  field  at  any  point  is  obtained  by  the  use  of  the  Green’s 
function  for  the  infinite  domain 


E(x,u>)s  ^ 


This  formula  is  an  integral  equation  in  the  scattered  part  of  the  e’ectric  field,  the  scatter- 
ing itself  being  due  to  inhomogeneity  in  mass  density  of  the  medium.  A few  important 
applications  of  this  formula  to  the  interaction  of  acoustic  and  electric  fields  are  discussed 
by  Morse  and  Ingard  (Chap.  XIII).  They  are  (a)  the  Debye-Sears  effect  (b)  Bragg  scat- 
tering (c)  Brillouin  scattering.  These  topics  are  briefly  reviewed  here. 

A.  Debye  — Sears  Effect 

Let  the  perturbation  of  the  mass  d^  -sity  be  due  to  a periodic  progressive  ultrasonic 

fa) 

wave  travelling  in  the  positive  x-direction.  The  harmonic  frequencies  are  a>s  ' , and  the 
.(*>) 

wave  numbers  are  Rs  J ,y\=  •>  i It  is  assumed  that  the  width  of  the  £<>und  beam 

is  small. 

Now  let  the  incident  light  be  a monochromatic  plane  wave,  E = B*  cos  ( |^e ' - ) 

where  Kj  »(t  C6i  <X  f f (*><  i k Cfiii ) | ( A , jj  , 7 , being  direction 

cosines,  and  j^e  = 0-/cB  . This  light  is  scattered  by  the  "grating  effect"  of  the  periodic 
ultrasonic  wave.  Choosing  for  example  then  n’th  Fourier  component  of  this  wave,  and 
writing  the  mass  density  as  sinusoid;  ^ . 

Y.’  " ( fesx  - t) 

one  then  seeks  a solution  to  Eqs.  (3)  and  (4).  In  general  the  solution  process  is  difficult 
to  carry  through.  Approximate  answers  are  obtainable  by  the  following  procedure.  First 
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it  is  assumed  the  perturbation  is  so  small  that  the  electric  field  in  the  source  term  ^ 

( £,£i»  is  effectively  the  incident  light.  This  is  the  Bom,  or  single  scattering,  approxi- 
mation. Secondly  it  is  assumed  that  the  width  of  the  sound  beam  ( = Z L ) is  so  small  that 
multiple  reflection  (as  from  parallel  layers)  is  negligible  in  establishing  the  scattered  electric 
field.  Thirdly,  both  electncal  and  mechanical  fields  are  taken  as  2-dimensional  (i.e. 
coordinates  x,  z),  so  that  'Ke  * t ^ ^4  and  * fv/i  • Thus  using  the 

2-dimensional  Green’s  function  (74  ) ( |Ks  I v/(#-Xft)*+  (f  -fcf  • takin« 

to  be  parallel  to  the  y-axis,  one  arrives  at  the  following  approximation  to  the  scattered 
field  due  to  the  n’th  harmonic  components  of  the  periodic  ultrasonic  wave, 

E = E*  •£  - (C.7) 


where 


1 [ &L/c 

C (A ] (Si l/c  Sf  ^ 


K±  -jr  -•  ” 


‘^1  . [m  ) — - --  - ■— 

4^°^  * |&  I ***  ± *s , (C10> 


The  symbol  £),  is  the  angle  of  scattered  light. 

In  words:  each  Fourier  component  of  the  periodic  mass  density  perturbation  induced 
by  the  ultrasonic  wave  scatters  the  light  (incident  at  angle  b ) into  two  plane  waves 

a(fr\)  f •*). 

^ and  (7j  respectively  where 


Aun  5 'lilt Jg)-  &t<V)  ,0.11) 

When  « 1/.  (as  is  usually  the  case)  the  -cattering  angles  are  given  by 

' *»•*  . (C-1J 
Let  ^ be  the  fundamental  component  of  the  sound  wave,  and  ^ the  wavelength  of 


(C.12) 


the  incident  light  wave,  then 


!,  then.  , a 

^ 6 * ,,  4s. 


(C.13) 


Thus  on  each  side  of  the  direction  of  the  incident  beam  will  be  found  two  diffracted  plane 
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waves  for  each  harmonic  of  mass  density  appearing  in  the  periodic  ultrasonic  beam.  The 
frequencies  of  these  waves  are  -fl  (“up  Doppler”)  and  J2  ("down  Doppler”). 

Their  amplitudes  are  proportional  to  the  amplitude  of  the  n’th  order  harmonic  of 
the  ultrasonic  wave  (which  diminishes  with  increasing/n. ).  For  an  arbitrary  periodic  dis- 
tribution of  mass  density  a family  of  plane  waves  rill  appear  in  space  distributed  in  angle 
according  to  Eq.  (6C).  These  diffracted  waves  mane  up  the  Debye-Sears  effect.  It  is  a 
notable  feature  of  this  phenomenon  that  the  diffracted  waves  appear  for  any  angle  of 
incidence  $ . 

Bragg  Reflection 

A plane  wave  train  of  monochromatic  light  ( wavelength  ) crosses  a plane  wave 
train  of  monochromatic  sound  (wavelength)^)  at  an  angle  of  incidence^.  The  sound 
of  train  constitutes  parallel  layers  of  density  from  which  the  light  is  reflected.  We  con- 
sider three  wavefronts  of  sound  at  phases  0,  2 n and  4 rt  radians  (that  is  2 X^long),  and 
a single  wavefronts  of  light  at  zero  phase.  The  light  is  first  reflected  by  the  zero  phase 
wavefront  of  sound  at  an  angle  8'.  In  order  for  a second  reflection  to  occur  the  light 
must  reach  the  second  layer,  travelling  a distance  tic.  It  is  seen  that  must  be  equal  to 
or  less  than  XXj.Thus  the  angle  8 J is  given  by 


■ a-  - tie 

^ e.  -- 


(C.14) 


Now  we  choose to  be/^(one  wavelength).  This  sets  the  requirement  that  the  second 
reflected  wavetrain  be  in  phase  with  the  first  reflected  wave.  This  is  equivalent  to  the 
geometrical  requirement  that  the  wavefronts  1 and  2 be  reflected  into  wavefronts  3 and 
4 shown  in  the  Fig.  Cl.  This  is  the  Bragg  angle  requirement  that 


(C.15) 


PIG,.  Cl 
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For  given  and  "X  it  is  seen  that  only  one  angle  will  satisfy  this  equation.  If  we  set 
8^  » 90°  the  light  is  reflected  exactly  backward,  provided  the  wavelength  of  light  is  twice 
the  wavelength  of  sound,  hi  going  from  wavefronts  1 and  2 into  wavefronts  3 and  4 the 
light  wave  is  scattered  through  a total  angle  6,  " 2 0’.  However  the  wavelength  ^ is 
created  by  a wave  moving  at  velocity  0 , and  characterized  by  a frequency  1 ^//tj  . 
The  frequency  of  the  light thus  Doppler  shifted  by  an  amount  — that  is, 


by  an  absolute  amount, 

1 (C.16) 

This  equation  gives  the  spatial  and  temporal  relations  between  all  the  vrriables.  If,  for 
example,  we  are  to  observe  reflection  of  light  from  a sonic  wave  train  at  “scatter”  angle 
9,  then  for  a selection  of  frequency  of  light  J^we  must  employ  a frequency  of  sound  given 
by  the  above  equation.  In  a typical  case  (Morse  and  Ingard,  p.  819)  j-^is  taken  as  10*5 
Hz,  CE  as  3 X 10®  Ns'*  and  ss  10®  Then  the  observation  of  back  reflection 
(at  Bragg  angle  180°)  requires  an  ultrasonic  frequency  of  some  10*®  Hz.  If  ^ is  less 
than  10®  Hz  the  reflection  angle  is  very  small,  that  is,  the  reflection  is  very  nearly  in  the 
same  direction  as  the  incident  wave. 
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Appendix  D 


Thermal  Modes  of  Fluid  Motion  (Ref.  “Theoretical  Acoustics”,  Morse  & Ingard,  Chap.  XIII.) 

The  generation  of  acoustic  pressure  and  the  diffusion  of  acoustic  temperature  7* 
form  a coupled  system,  described  by  the  equations, 

(D.l) 

7T-  /,  f <D-2> 


in  which 


ui  mucu  / 

£.  A./ 

1 Cl  ^ ** 

£ ' £ 

1 ■ jpi* 

a-  y> 

C;(f +p)f* ; ^ % '> 

Here,  C^,  Cv  are  specific  haats  at  constant  pressure  and  constant  temperature,  tj  is  the 
bulk  viscosity  modulus,  jl  is  the  shear  viscosity  modulus,  £ is  the  coefficient  thermal 
expansion,  f(  is  the  thermal  conductivity  and  Kr  is  the  bulk  (elastic)  modulus  at  constant 
temperature.  Eqs.  (1)  and  (2)  are  acoustic  equation*.  A more  general  equation  of  fluid 
dynamics  is  the  Navier-Stokes  equation, 

f *-£•*  -f*  (o.4) 

in  which  and  ‘Mf  are  the  longitudinal  and  transverse  parts  of  the  fluid  velocity,-^  is 
the  acoustic  pressure,  ^ is  the  equilibrium  fluid  density. 


I 


Three  specie:  of  fluid  motion  are  predicted  by  Eqs.  (D.ll  thru  (D.4).  These  are  (1) 
longitudinal  waves  of  acoustic  pressure  £ (■  Eq.  D.l),  (2)  transverse  shear  waves  (=  Eq.  D.4) 
(3)  thermal  diffusion  of  temperature  {*  Eq.  D.2)  which  then  gives  rise  to  acoustic 
waves  (“  Eq.  D.l).  We  will  discuss  thermal  motion  and  the  accompanying  thermal  modes 
of  sound  generation. 

In  the  thermal  mode  «stT  is  very  much  greater  than  Jb.  The  diffusion  of  temperature 
in  this  case  is  given  by  Eq.  (2).  To  apply  this  equation  we  first  assume  the  acoustic  pres- 
sure on  the is  negligible,  and  take  the  temperature  pertubation  to  have  the  space  and 
time  dependence  of  a plane  highly  damped  travelling  wave,  CC  typ  . Then, 

to  first  order,  _ 

Ar-feV 

V ‘ K 1 (D.5a) 

(D.Sb) 


7ft  = ~ik>r 

From  Eq.  (2)  it  is  then  seen  that  the  effective  square  of  the  wave  number  is  complex, 

4.*-  i-V 

Aci  <"*> 

The  acoustic  pressure  corresponding  to  this  diffusion  of  temperature  is  found  from  Eq.  (1) 
in  combination  with  Eq.  (5a),  (5b),  and  (5c),  i.e. 


i&a  (A-Qr 


(D.6a) 


The  effective  perturbation  of  pressure  is  not  this  alone  but  the  combination  p-  <xT . 
The  term  oiT  describes  the  reduction  in  pressure  due  to  heat  conduction  (i.e.  describes  the 
departure  from  the  assumed  adiabatic  condition  of  wave  propagation).  Corresponding  to 
wives  of  acoustic  pressure  are  waves  of  acoustic  density  ft  which  to  first  order  are  related 
to  each  other  by  the  formula  m<yfi  or 

P.  * oft. 

' c;  c/ 

The  perturbation  in  fluid  velocity  associated  with  thermal  motion  is  found  from  Eq.  (3) 
by  neglecting  the  viscosity  term  on  the  /i-i-i  since  this  term  will  yield  only  a second- 
order  effect  of  temperature  on  velocity.  Thus  one  finds 


(D.6b) 
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k 


% = ”7"  ^4*4-)  (D.6c) 

The  set  of  equations  comprised  in  (5)-(6)  form  a convenient  description  of  the  thermal 
mode  of  fluid  motion  in  terms  of  the  parameters  and  . 

Decay  Rate,*  of  Temperature  and  Pressure  Fluctuations 

Assume  the  modes  of  motion  noted  above  arise  from  random  fluctuations  in  pressure  (Af> ) 
and  temperature  ( aT).  Temperature  fluctuations  decay  at  the  rate 

K v’Ut)  »f  C 1(f) 

r at 

(D.8) 

Assuming  the  spatial  dependence  of  aT  has  the  form  exp  iVp,  we  reduce  this  equation 
to  the  form 


3 (at) 


at 


-{3(*T 


\ 

> > 


Li* 

?<;1 


(D.S) 


Hence  the  time  decay  of  (spontaneous)  fluctuations  in  temperature  of  a fluid  follows  the 
law  ATof  exp  (-pt ),  in  which  p is  the  temporal  decay  constant. 

Now  let  there  be  « pressure  fluctuation  (spontaneous)  of  amount  Af  (^>4)whose 
spatial  dependence  is  cos  ^ . Th'j  mode  of  fluid  motion  if  found  by  solving 


V "Jb*  ofy  = 0 


(D.10) 


for  the  complementary  solutions  in  time.  The  result  is 

dj>(h,t)°ce  [c*i (<j\r+Ht)+<wfyr -<*$1  (d.id 

which  describes  .wo  waves  traveling  in  opposite  directions  (at  angle  ^<j»r)  at  frequency 
^ ■ <JCS  where  ^ (<t>)  satisfi»i  the  dispersion  formula  obtainable  from  (D.10).  The  tem- 
poral decay  constant  is 


109 


^V*/, 


\ 39  e / 2 

which  is  also  obtained  from  the  complementary  solution  of  (D.10).  The  time  t=  d will 
be  called  the  effective  life  time  of  the  wave. 

Under  forced  drive  Eq.  (D.10)  has  a forcing  function  the  By  Fourier 

transformation  of  time  (vnder  zero  initial  conditions)  the  acoustic  pressure  is  seen  to  be 
given  by 


/ , " ^ F (h}U)) 

-|>  (A , 6u  )-  7 ~ 


SlJ1-  + (t)ai 


(D.13) 


at  the  lifetime  of  4^(i.e.  its  effective  influence  on  the  spectrum  of  ) is  o L 

order  (dimensions  of  seconds).  The  bandwidth  otja  under  forced  drive  is  therefore 

of  order  , that  is,  of  order  0^ . 
o( 

In  sum:  Spontaneous  fluctuations  of  pressure  are  propagated  as  damped  traveling  waves, 
and  spontaneous  fluctuations  of  temperature  are  diffused  as  critically  damped  waves.  A 
spontaneous  spatial  sinusoidal  distribution  of  pressure  at  wave  number  propagates  away 
as  two  damped  waves  in  opposite  directions  at  frequency  Each  propagating  acous- 

tic mode  is  associated  with  a propagating  mass  density.  The  mathematical  form  of  these 
damped  waves  have  the  same  form  as  the  forced  drive  of  a damped  harmonic  oscillator 
(Eq.  (D.10)). 

We  next  consider  the  interaction  of  light  waves  with  mass  density  fluctuation  associated 
with  thermal  motion  of  a fluid.  Let  Ap(h,-0  be  a spontaneous  fluctuation  in  mass  density. 
Choosing  (as  before)  the  thermodynamic  variables  the  pressure  and  temperature  we  write. 


( D.14) 


We  consider  the  effects  of  temperature  and  pressure  separately.  We  first  take  the  spatial 
Fourier  transform  of  the  pressure  and  temperature  fluctuations,  and  write. 


(A,t)  (D.15) 

J tT(k,t)e~\  ( D.16) 


where  the  volume  of  integration  is  finite  (as  it  must  be  in  any  real  case).  Treating  £p  as 
a random  variable  we  next  calculate  the  temporal  autocorrelation  Ov/Obf  *be  pressure, 

If  (k>  tV  < A (>(*,*«*)  (D.17) 

in  which  the  symbols  ^ represent  a temporal  average  ( * integration)  where  the  integration 
is  over  finite  time.  The  intensity  of  the  pressure- fluctuation  is  obtained  by  setting  T-  0 . 
Similarly  the  autocorrelation  of  the  temperature  fluctuations  is  given  by 


■Rr  (fe,  r ) = <iT(k,  in)  aT* '('*,<)>  ( D.18) 

Assuming  further  that  the  pressure  and  temperature  fluctuations  are  uncorrelated  we  can 
find  the  autocorrelation  of  mast  density  T^(k,-ty  by  addition,  i.e. 

a , A 


K?  (t’f) 


(D.19) 


From  previous  analysis  (see  Eq.  (0.11))  the  temporal  character  of  a pressure  fluctuation 
has  been  i-  volved  as  the  response  of  a damped  simple  harmonic  oscillates  for  each  wave- 
number  i , and  associated  frequency  &. ■ ks  C . Thus  one  can  write 

^ ^ itj 

CM  (D.20) 


Similarly, 


1 Rr(k>t)  r e 


( D.21) 


These  formulas  will  be  used  in  the  dirc-ission  of  Brillouin  scattering  (see  Eqs.  E6,  F.7,  etc). 
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Appendix  E 


Brillouin  Scattering  (Ref.  “Theoretical  Acoustics”,  Morse  & Ingard,  Chap.  XIII.) 

The  interaction  of  light  waves  with  inhomogeneities  of  mass  density  in  a fluid  has 
been  discussed  in  Appendix  C in  the  case  of  Debye-Sears  effect.  In  the  present  applica- 
tion the  perturbations  of  mass  density  are  assumed  to  be  due  to  changes  in  pressure 
and  to  be  spontaneous  (i.e.  random).  Let  be  the  random  fluctuatioi:  in  pressure 

Its  spatial  Fourier  transform  is  given  by 

^ (L-O*  J -Va. 

where 


(E.1) 

(E.2) 


Assume  the  volume  V of  pressure  inhomogeneities  to  be  finite,  and  assume  the  incident 
light  to  be  a plane  wave  - F0  /vp  t (J£ To  Using  the  Green's  function  for  un- 
bounded space  and  approximating  the  local  electric  field  by  E'(  (“Bom  approximation), 
one  formulates  Eq.  (C.6)  for  the  scattered  wave  by 

ik-r*  * 


ci 


iA.f ) = ^ ^<0  pVopt  A |>(t  i)e 

in  which  the  dimensions  of  Q are  meter-1  sec-1,  and 

cv  ±-  ♦!££!) 

4 n | a -pi 

■/r##  , , 
In  the  far  field  * h.  and  r‘  l/c  ~ A - , where  |<f  = j K | 

forming  all  mathematical  operations  on  (E.3)  one  arrives  at  the  formula 


(E.3) 


(E.4) 


Per- 


(E.5) 
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in  which  we  have  set  k ~ . Now  the  intensity  of  light  waves  (J)  is  given  by  the 

Fourier  transform  of  the  time  averaged  autocorrelation  of  ^ , (brackets  indicate  time  aver* 
age,  or  division,  by  time  interval), 


L&r 


I („,K-K  si)  * f*  <£>(*.,&, tit)??/*  >*)>e 

„ . J /ae\1/EiSU'j2  f~^U-V/r) <*•» 
e1  U^-wac*'  J*  e e 

in  which  0 is  the  effective  duration  of  the  pressure  inhomogeneity.  Hie  autocorrelation 
of  intensity  of  pressure  has  the  same  form  as  that  of  a damped  harmonic  oscillator,  as 
noted  earlier,  (see  Gq.  D20), 

-d|tl 


i?A?(  <'-*>?)  asgr 


(E.7) 


Thus,  the  intensity  of  the  scattered  light  is  given  by 


„ -■*  ^ 
1 ^ c»stvr  e 


Hie  integral  is  standard.  Hie  final  result  is, 

r. 


ic  integral  u lumumxu.  inc  unu  result  is, 


of 


* | d7i  [a1-  (ft-»s)T  1 


(E.8) 


(E.9) 


In  words:  the  intensity  of  scattered  light  in  the  direction  at  distance^  due  to 
spontaneous  pressures  in  a Quid  is  the  sum  of  two  contributions,  namely,  the  two 
Brillouin  lines.  The  first  line  is  at  frequency  _fl  = SI*  and  the  second  line  is  at 
Si-. The  bandwidth  of  these  lines  is  cl  , and  the  corresponding  lifetime  ic  of  . 
When  mass  density  fluctuations  are  due  to  fluctuations  in  temperature  a similar 
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mathematical  treatment  leads  to  the  result  that 


This  is  the  scattered  Rayleigh  line  which  appears  in  the  direction  , centered  at  fre- 
quency!}. (-  incident  frequency)  with  bandwidth  ft  (or  effective  lifetime  of  ft  ). 

The  total  intensity  of  scattered  light  is  the  same  of  scattering  due  to  pressure  and 
temperature.  However,  by  using  an  appropriate  equation  of  stato  relating  pressure  per- 
turbations in  the  fluid  to  temperature  perturbations  one  can  reduce  all  Brillouin  scattering 
to  thermal  motions  of  the  fluid. 

The  intensity  of  light  derived  earlier  was  expressed  in  Lorentizian  form  which  is 
e 'valent  to  the  form  of  the  response  of  a simple  harmonic  oscillator  to  a random  exci- 
tation. An  analog  of  Eq.(D13)  is  the  RLC  circuit  with  forced  (volts)  at  frequency 

■f.  For  a random  emf  the  random  current  is  given  by 

<£>  = <*» 

**' (!-«-£)*  (EU» 


Now  if  the  random  emf  consists  of  fluctuations  in  the  resistor  R (■  thermal  noise)  it  was 
shown  by  Nyquist  that  in  the  small  frequency  interval  between  *f;  £ throne  has 


--  4-R 


(E-12) 


in  which  k ■ Boltxman’s  cons  tan  t^T  ■ absolute  temperature.  The  mean  square  current 
accompanying  thermal  noise  is  therefore 


(E.13) 
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? 


4 


\ 


t. 


Here  ot  ” ^/2L  and  01"  * is  the  effective  life  of  the  fluctuation,  kT  is  energy  per  degree 
of  freedom  and  kT«f  is  the  power  in  the  fluctuation. 

Instead  of  an  electrical  problem  one  can  consider  a single  degree  of  freedom  mechani- 
cal system  of  a mass  M,  spring  K,  a-'d  viscous  resistance  R.  The  equivalent  Nyquist  for- 
mulas for  random  (spontaneous)  noise  in  the  system  are 

<FV»  * a'RfcTt 


in  which  p is  the  force  impulse  (Newton  sec)  which  operates  during  time  t sec,  and 
<tl>  - t iTt  (E.14b) 

H 

in  which  x is  the  amplitude  of  free  Brownian  motion  in  time.  Thus, 

<Fte$)>  = (E.14c) 

is  the  (time)  average  of  forces  with  frequencies  between  f and  f + In  words:  A 
system  at  temperature  T exhibits  tandom  impulses  during  any  time  interval  t,  the  meui 
square  average  of  which  is  given  by  Eq.  (E.14a).  These  impulses  ect  at  an  extraneous 
particle  causing  it  to  undergo  Brownian  motion.  The  motion  is  damped  by  viscosity. 


Experimental  Observation  of  Brillouin  Lines 

Brillouin  lines  are  observed  by  high  resolution  spectroscopy.  A modem  technique  is  the 
laser  interferometer.  According  to  Eq.  C.16  the  shift  in  laser  beam  frequency  d$e  /$g  due 
to  the  velocity  Csof  the  acoustic  wave  generated  randomly  by  the  Brillouin  effect  is 


4Ag 

Xe 


in  which  ^ is  the  speed  of  light  in  the  fluid,  and  is  its  wavelength.  Assuming  0*  I®, 

\ I 5 M|  • 

= 1.5  x lu  Ms”  , (^ “ 2.26  x 10  Ms  (in  water),  it  is  seen  that 


t A« 

s 


_s 

».i  * <c 
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The  wavelength  separation  ^depends  on  the  laser  wavelength  to  be  used.  Let  Xg  = 3. 84*  /0 
M (in  water).  Then 

• 1 

aA6=s.ix/o”  M a 0.05/4,  A$t-  ts  " ?• 6 * io  W*  • 

This  is  the  magnitude  of  Brillouin  effect  that  can  be  detected  with  visible  light. 
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Appendix  F 


A Case  of  Scattering  of  Sound  Waves  by  Inhomogeneitiea 

Let  there  be  a region  R (in  a fluid)  in  which  the  change  in  compreasibility  is  de- 
scribed by  the  factor  (h  ,i  ) 

0 ovJsJt'P.  (F.l 


The  acoustic  pressure  satisfies  the  equation 


Assuming  > «nd  using  the  Green’s  function n ^»^for  unbounded  space, 


it  is  seen  that  the  total  field  is 


"H*  ) * fi  (h  )+  ’fi  (&  ) 


in  which 


Q ( h.  | r»  \ s I — 6 

5 ~ 


t^  lc-c°  I 


In  the  far  fieid 


n / | . I ~ i k /t  • A 


(F.4) 


The  scattered  field  is 


fsU)=i U ~ dn 

4nft  ^ 


(F.5) 


Using  the  Bom  approximation  we  set  the  total  pressure  ^ (/la)  to  be  the  incident  pressure 
incident  beam  be  finite,  and  choose  its  form  to  be 

f,  i'f.)  = A(rj)>if  i(kr.C‘) 


Then, 


, i • u /.  , ( *»- 

)•  i-«l  f/4(£.)X  (c-,+)e  s/V. 


4rt/u 


(F.6) 


If  the  inhomogeneitiea  are  in  the  form  of  particles  the  volume  integration  can  be  replaced 
by  a sum, 


. i 2 <r—  i ^ ^ 

e-  y *« 

4 It  S. 


AV 


(F.7) 


Finally  we  allow  to  be  a function  of  time,  We  see  then  that  the  incident 

wave  is  phase  and  amplitude  modulated. 

When  there  are  many  scattering  inhomogeneities  we  write 
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tsU  )=  1 A(c»)1k(c*)e 

4nh,  _ j 


t ik-^)‘  G* .3 

- A 


(F.8) 


The  value  of  the  integral  depends  on  the  size  of  the  scattering  volume  A relative 

to  the  wavelength  of  the  incident  wave,  ^ . If  4 , then  the  phase  func* 

tion  (given  by  the  exponential)  is  constant  during  the  integration  process.  In  this  case 
(Rayleigh  scattering) 


ii  t kh  t ( ki  ) . lo*» 

* M=  2 * “ " ?(:*» ) 

’ 4KU-  *, 


(F.9) 


where, 

~P  ( j ( A (rov)itc  (/»*  ) d C* 

J (F.10) 
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Appendix  G 


Volume  Scattering  Due  to  Perturbations  of  Fluid  Density 


In  volume  scattering  due  to  perturbations  on  mass  density  the  scattered  field  obeys  the 
wave  equation 


V*E  _L  " 2 i bE 

-s ' 5 w' 


Ag  i fa  ) Af  t \ 4T  . \aH 

Ur^ft  MSfJft 


in  which  AM  represents  chemical  (=  mass)  concentration. 
Consider  only  mass  changes,  Then 


v*e,  - 44 


Then  let 


*ri  (r,<)  ■ C iH  ( 


Assume  next  that  the  electric  field  0n  the  h- Cj  is  the  incident  field  E o t [ K*Jf 

and  use  the  free-field  Green's  function  (fan  )(  0 ' |CrC°i/c  ) to  solve 

the  wave  equation.  Thus 


if>0»  j[  j (r 


or 
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Esu,tl*  iAi\ 

~ ‘ ?AmL  <*  JJ 


-*»  n»*  • /pT  *-  P _ . ~«  ~ ~ (G.6) 

As  before,  assume  jf-ro(  v /*  -(j^'/A*  *"*  ) l‘  • Integration  over 


gives 


t k'£ 
» ■*  ^ 
a 


Es  M--  ~ V* 

C»  “o\pnpr ' 


(G.7) 


Now  of  the  volume  V(Q>^is  large  enough  the  integration  over  ifi  is  an  infinite  integration 

over  an  imaginary  exponential  which  yields  a delta  function 

Thus 

iKk-iOt 


.3(4  +K-Z‘)^ 


e<  (M I k1- y ) * - &9-!£  -ikj  a ff  (£-*,<)€ 

4rtA<£  VoyHhj 


(G.8) 


(G.9) 


The  autocorrelation 


ion  "5^  is 


■Re  B (t» 


dimensions:  field  ^ 


or 


oM(K-KJo^,t(v-v^-o>  ,o,o, 

The  value  gives  the  intensity  of  the  scattered  electric  field.  The  power  spectrum 

of  f.  (dimensions:  field ^ x sec^),  is  the  Fourier  transform  of  over  the 
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finite  interval  T, 


5JT_ 

T 


^E  (t) 


or 


1 T/x 

jt)&H*(l£\ K jin)  M 

% 


(G.ll) 


Now  the  spectral  intensity  density  f the  electric  field  is  the  limit 


Jte  (4)  * U F |E(fl')l 

T-*»  T 


(C.12) 


Hence 


/E.&V  \’ 

\ ' 

\47tACf1  / 

U,T. 

i q x- 

*ir  (dimensions:  field 2 x sec) 


(G.13) 


Ip  words:  If  we  observe  the  scattered  electric  field  in  direction  as  a function  of 

frequency  we  can  find  the  autocorrelation  of  the  density  field  by  Fourier  transformation. 
Alternatively,  if  we  know  the  autocorrelation  of  the  mass  density  field  we  can  calcu- 
late the  electrical  spectral  intensity  density  (dimensions:  electric  field^  x sec). 
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Appendix  H 


Scattering  of  Light  from  Particles  Suspended  in  a Fluid 
(Edwards  et  al.  J.  Appl.  Phys.  42,  837,  (1971)). 


Hie  scattered  electric  field  is  given  by 

U ,0  ■ E.  Si.  i&\  f 


i [(*-*')•£.  -A*7 


Sr 


(H.l) 


Let  us  now  consider  individual  scatters  rather  than  a continuum.  To  do  this  we  replace 
the  volume  integration  over  fi.  by  a sum  over  individual  scattering  particles,  i.e. 


§>M=  6 1 


A - 


e -ftc 


A- 
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electric  field 
acoustic  pressure 

(H.2) 


in  which  A ^ is  the  volume  of  a scatterer.  Furthermore  let  us  take  the  location  vectors 
to  be  functions  of  time  (namely  we  take  the  particles  to  be  in  mcdon)  and  write 


Ae*  ("i ) * ha^,  (o)  + M0*  (*)+?* 

This  reads:  the  n’th  particle,  initially  at  moves  in  time  t relative  to 

the  fluid,  which  is  itself  in  motion  with  velocity/. 

The  density  perturbation  M f is  a function  of  space  and  time.  Let  us  take  an 
arbitrary  time  to  be  the  origin,  and  write 


Also,  at  time  t 

— ^ /Tl  I (H-*l 

* t fo*(o)  + -ntf)]  aVA^) 

In  this  analysis  we  require  the  intensity  of  the  scattered  electric  field.  We  can  ob- 
tain this  by  use  of  an  autocorrelation  function  defined  as 


Kh(t)  = « £>  (f,o)  M(rjt)» 


(H.5) 


in  which  the  double  angle  brackets  indicate  ensemble  averages  over  (1)  initial  position 
fq^  (0)  ; (2)  random  displacements  /I  • To  form  an  ensemble  average  over  initial 
position  we  introduce  (see  Edwards  et  al)  a probability^  ( ft  (#))  per  unit  volume  that  a 
scattering  center  will  be  found  in  volume  1/  . To-  form  an  ensemble  average  over  displace- 
ments'^^) we  introduce  a probability  per  volume*^  h^ft)/!*)  that  the  nth  par- 
ticle will  move  units  in  time?.  The  autocovariance  of  the  scattered  field  then 
becomes 

■Re  ft)  * A A <Rt  e \z%(  ^ <t),tW  -t  hfl 

0r\  I (H.6) 

11  $U(Uu))*^-tr^^,)^'r7  tV\(r^/o)to) An(f^(o)^r^fx) 

t VT ) d V** (*)d36rn  fo*  (»]  % (dimensions:  field 2) 

We  note  that  only  one  sum  o*<n  is  used  In  the  product  2"Z  required  by  the  definition 
of  ^ since  the  contributions  all  vanish  in  view  of  the  assumed  statistical  indepen- 
dence of  the  scattering  centers.  This  equation  can  be  written  in  another  way  by  defining 
an  autocorrelation  function  for  the  density, 

urt)  - 

J ^ ? r -A 

'iMlr^(e)  iAt**<t)*vr)]  d UM  [4V{Cn)J 

(dimensions,  acoustic 

Thus,  pressure2) 
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(t)s  /4A  t ^ Z d*.  (* ~L  4ft*  ^ 

' ( ^l^'>’  ^)  * y^)  “ 1 (}£~f  ) mJft  ^ Af*  (H.8 


TTie  power  spectrum  ) I of  the  scattered  field  is  obtained  by  Fourier  transfor- 

mation of  (V^  over  a finite  interval  7 , 

TU  i SLt 


|Mu')|2-  7?e  MeV 

T J 


(H.9) 


-r/v 


The  spectral  intensity  density  *rP%  the  electric  field  is  the  limit 

3h.  (A’)  /E(fl’)f*  f!>ifr)eia$fr 

C Tina  T ...  . . _ . 


(dimensions:  electric  field ^ x sec) 


In  words:  a time  record  of  the  received  signal  observed  in  direction  K*-  K is  used  to 
form  an  autocorrelation  function!^  ft')  where  7 is  the  time  shift.  After  ^gft)is 
formed  we  then  set'f-d  to  obtain  the  intensity  of  the  received  electrical  signal.  This  in- 
tensity is  a function  of  the  autocorrelation  of  the  density  field.  In  the  absence  of  an 
acoustic  signal  the  autocorrelation  ^ (^will  correspond  to  the  autocorrelation  of  noise. 
The  ratio  of  signal  power  to  noise  power  will  have  the  form 

*6  M, 


s. 

N 


'4S_ 


(H.10) 


If  we  set  a threshold  j 
trie  signal  as, 

Mm 


"fa  (°)n 

we  establish  the  minimum  detectable  intensity  of  elec 


'iQ 


(H.ll) 


The  f unction  "i^(f)is  an  amplitude  weighting  of  the  laser  light, which  depends  on 
the  position  of  the  scatterer  relative  to  a characteristic  size  (or  dimension).  For  example, 
if  the  characteristic  dimension  is  Lc  and  the  position  /u>  then  we  would  expectP  to 
be  very  small.  The  vector  position/t  is  a function  of  time.  If  / is  the  velocity  of  the 
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fluid  then  \/t is  a component  of  A . Hence  P(r)  ■ P(r(t)). 
Now  the  autocorrelation  of  the  scatter  field  is 


Re  « MM  £*(£/!■)» 


(H.12) 


in  which  >)  means  averaging  over  (1)  initial  position  of  the  scatteren  (2)  over  dis- 
placements relative  to  fluid  (*  random  deviations).  Let  Q (&\ ft  it’J  be  defined  as 


Q ( Af -tvt)r  (W.  j 1^^r/o))'P[r'o)  +fr+yr] 

-Too 


(H.13) 


«xp 

If  one  selects  a.  (spherical)  finite  volume  as  a model  of  the  collection  of  scatterers  and 
writes 


•*<*)  - 4 


»/,  rL 

* _ 4TV 


(2rt<?yb 


(H.14) 


then 


0 - I«e~ 


"W1 


(H.15) 


The  ratio,{a  gives  the  residence  time  of  the  scattering  centers  in  the  sample  volume. 
The  effect  of  random  displacements  is  expressed  in  the  Lf  dependency,  conjoined  with  a 
probability  function!^  for;  t)which  shows  the  probability  that  a scattering  particle  has 
moved  bt  in  time?  . Thus  the  generalized  autocorrelation  c§  (^C+yTffc  )vi  } 

3 (tnn.k)*  Q(4TVrMe~l*‘ir 

K (H.16) 


(the  minus  sign  in  the  exponential  is  derived  from  the  conjugate  electric  field).  Thus, 
for  a continuum  of  particles,  the  autocorrelation  function  is 


7 h 


(4,t)=  OSiuf. 


(H.17) 
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The  power  spectrum 


m S (*,«),  then  obtained  by  Fourier  transformation 

iktv)  = £ 


(H.18) 


In  words:  The  detection  process  for  measuring  the  scattered  electric  field  is  modeled  as 
a harmonic  oscillator,  frequency  f£‘¥  , damped  by  two  terms^'ha-  the  hf  effect  due  to 
random  motion  of  the  particle  relative  to  the  average  motion  of  the  fluid,  and  ft effect, 
arising  from  finite  transit  time  of  the  particles  in  the  scattering  volume. 


Assume  the  Af  effect  can  be  modeled  on  the  theory  of  the  random  walk,  i.e. 

{gf 

in  which'P  is  the  diffusion  coefficient  of  the  particles  ir  the  fluid.  Choosing  the  ampli- 
tude weighting  function!^  to  be  given  by  Eq.  (H.ff)  (in  which  a scale  size  <T  is  specified) 


it  is  seen  that  the  power  spectrum  or  the  scattered  field  is 


T>rf> 


-ck-jf  ,,  4. 

c ***  tf*  e 'd  (at  ) <k? 


(H.20) 


We  first  allow  Y to  be  zero.  Then,  using  the  convolution  theorem  for  Fourier  transforms 


one  arrives  at 


vt:vci  , ao  A r-  l 

- <3  V 1 ]jk  Jt  (HJ21) 


An  approximate  spectrum  can  be  obtained  by  expanding  the  diffusion  term  in  a Taylor 
Series  about  l(- 9 , and  then  integrating.  The  result  is 

$(*,»)-<***  f A:i  + Vk j 


£c/«  Jo 
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We  next  allow  V to  be  finite,  and  proceed  again  to  expand  the  diffusion  term  in  a Taylor 

series  about  K *0  . For  small  volumes  of  interrogation  we  can  use  only  the  first  term  t 

in  this  expansion,  namely, 


Appendix  I 


Intensity  of  Scattered  Light  From  Acoustic  Signals  (Edwards  et  al.  Model) 


The  intensity  of  scattered  light  in  a medium  of  dielectric  constant  £ , magnetic  per- 


meability if , is  given  by  the  Poynting  vector,  5 r — 1 1 E X H , or 
1 ~ 41 


S * vr  V JL 


8 - 1 

in  which  Cg.  = 3*  (0  MS 


' y - 


* 47T  11 


(dimension:  E,  volt/meter 

£ , coulomb/(meter  x volt) 

-yf . volt  coulomb/meter^ 

, meter  volt/coulomb 

C,  V meter/sec 

S,  coulomb  volt/(sec  x meter^)) 


Thus  in  MKS  units, 


cv~  /T  E * 


Now  the  scattered  electric  field  is  a random  function  of  particle  motion.  We  therefore 
will  obtain  the  quantity  E by  finding  the  autovariance  of  the  scattered  field  (*  Kg  It)) 
between  two  moments  in  time,  ti,i,  (note  'S'  ),  using  an  appropriate  spatial  aver- 
aging to  average  out  the  initial  position  /^,(»)  , and  random  motion  ) over  all 

space.  Thus,  according  to  Eq.  (H6), 

O0 

1?£hV  A A*<fce'lQ*T  \ Z p*t  )**[>--  [(<-*)•  ^ b)] 

r a-3) 

* K-K  )-  (yr-t 
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in  which 


*«„  Urm(r)+  vrAf^iosr): 

f ^ ft) +tv 

The  magnitude  of  intensity,  averaged  over  random  time  history  of  particle  motion  (that  is, 
ensemble  averaged  through  definitions  of  probability  functions  7^  ( 4^^  and  Tk(r'J°0 
and  average!  over  time  shift  T , is, 


|$  |--  £a  ll  %(*) 

'•  4it  > A* 


(1.4) 


47C  / /* 

The  spectral  intensity  density  (Q,  j of  the  scattered  light  is 

% fe.7)  * j ft  (r)e^jr  as) 

°°  (dimensions:  ) 

A? 

(see  Eq.  H.9).  The  actual  power  spectrum  of  the  scattered  electric  field,  calculated  over 
a frnit;  interval  T,  is 

i*  I es  (a*)!**  i ■** 

, ”T/i  a , (I'6) 

in  which  the  dimensions  of  J )l  V S M . 
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